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Abstract 

The Lagrangian of Quantum Chrome-dynamics is invariant under conformal transfor- 
mations. Although this symmetry is broken by quantum corrections, it has important 
consequences for strong interactions at short distances and provides one with powerful 
tools in practical calculations. In this review we give a short exposition of the relevant 
ideas, techniques and applications of conformal symmetry to various problems of interest. 
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1 Introduction 



Conformal symmetry has a long history. It has been known for more than a century that 
equations of electromagnetism are invariant with respect to the inversion of coordinates x M — > 
x^/x 2 . This symmetry is useful for solution of many nontrivial problems and has numerous 
applications in electrostatics. Later it was realized that the inversion symmetry is related 
to the scale invariance of the underlying U(l) gauge theory, at least at the classical level. 
The scale transformations and the inversion together with the usual translations and Lorentz 
rotations form the conformal group, which proves to be the maximal extension of the Poincare 
group that leaves the light-cone invariant. The conformal symmetry of a quantum field theory 
was intensively studied in the end of sixties and in the seventies, with the motivation coming 
independently from the Bjorken scaling hypothesis (see |T]) and the theory of second order phase 
transitions. The theoretical description of critical behavior in many cases can be obtained within 
the framework of conformal invariant theories. It turns out that conformal symmetry of a two- 
dimensional quantum field theory implies existence of infinitely many conservation laws and 
strong constraints on the correlation functions. In particular, spectacular results were obtained 
in statistical mechanical models in two dimensions. The subject of two-dimensional conformal 
field theories (CFT) has become very fashionable because of connections with string theory and 
quantum gravity (see, e.g., [2]), and has grown into a separate branch of mathematical physics 
with its own methods and own language. In some 4-dimensional theories with supersymmetry 
the /3-function vanishes (to lowest order or to all orders) and some featurs reminicent of two- 
dimensional CFT emerge. This subject attracted a lot of interest, enhanced in recent years in 
connection with the conjecture about AdS/CFT correspondence, see j3j for a review. 

On the other hand, the impact of conformal symmetry on the development of Quantum 
Chromodynamics (QCD) has been rather modest. In the QCD community, a widespread 
scepticism exists to the use of methods based on conformal invariance, for several reasons. As 
we will see below, the conformal symmetry of a quantum theory implies that its /3-function 
must vanish. In QCD, this is only the case for the free theory, a s — > 0, and in this limit it 
essentially reduces to the parton model as far as "hard" processes at large energies and large 
momentum transfers are concerned. The scaling laws of the parton model (e.g., quark counting 
rules etc.) can indeed be derived directly from the conformal symmetry of the classical QCD 
Lagrangian, but this formal connection is usually not what a QCD theorist is interested in. 
The proximity to the conformal limit is determined by the value of the strong coupling and can 
be controlled in hard processes. It is usually thought, however, that the study of this behavior, 
QCD at small coupling, is not really interesting for a "hep-th" theorist, while using specific 
techniques based on conformal symmetry in the QCD phenomenology is complicated and is not 
worth the effort. 

This review is an attempt to overcome this prejudice. We will restrict ourselves to the per- 
turbation theory and demonstrate that the structure of perturbative predictions for light-cone 
dominated processes reveals the underlying conformal symmetry of the QCD Lagrangian. Ex- 
posing the symmetry allows one to simplify calculations, obtain new results and new insight. It 
is worthwhile to mention that although the subject of conformal symmetry is old, the applica- 
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tions considered in this review are rather recent. Loosely speaking, the subject of perturbative 
QCD is calculation of the scale dependence of physical observables, which is governed by evolu- 
tion equations. The known evolution equations for structure functions, distribution amplitides 
and generalized parton distributions can be understood as the renormalization group equations 
for the light-cone operators and have much in common. Using conformal symmetry in order 
to understand the structure of these equations presents our main topic. In addition, we will 
consider the so-called BFKL equation for the perturbative resummation of large logarithms of 
energy, which is relevant in Regge limit. 

The outline of the review is as follows. Sect. 2 is introductory and is meant to provide basic 
information about the conformal group, conformal invariance in a field theory and its violation. 
In particular, the construction of conformal operators, the conformal operator product expan- 
sion and conformal Ward identities are important for QCD applications. Our presentation is 
necessarily heuristic and incomplete. More details can be found in numerous reviews and in 
particular we recommend j3J E] for the first reading. 

Sect. 3 is devoted to the conformal partial wave expansion of hadron distribution ampli- 
tudes. The approach is very similar to the partial wave expansion of scattering amplitudes 
in nonrelativistic quantum mechanics and aims at the separation of variables. In quantum 
mechanics, 0(3) symmetry of the scattering potential allows one to separate the dependence 
of the wave function on angular coordinates in terms of spherical harmonics which realize an 
irreducible representation of the symmetry group 0(3). In QCD the corresponding symmetry 
group is SX(2,IR) and it allows one to separate the dependence of hadron wave functions on 
transverse and longitudinal distances with respect to the scattering plane. The dependence on 
longitudinal degrees of freedom (momentum fractions) is included in "spherical harmonics" of 
SL(2, M.) while the dependence on transverse coordinates is traded for the scale dependence and 
is governed by simple renormalization group equations. Historically this expansion was first 
suggested around 1980 [UJ [7] and was instrumental for the proof of the QCD factorization (see 
8J for a review) for the elastic and transition form factors. This, essentially algebraic, program 
has been then extended to general multi-particle distributions jHj and is still in progress. For 
higher twists, it is important that conformal expansion is fully consistent with QCD equations 
of motion. 

In Sect. 4 we consider a general approach to QCD evolution equations which is based on 
the representation of evolution kernels as formal Hamilton operators depending on Casimir 
operators of the SL(2,M.) group [TOj- In this way the SL(2,"R) symmetry of the leading-order 
evolution equations is made manifest. Using this technique it has recently become possible to 
solve evolution equations for three-parton operators which govern, e.g., the scale dependence of 
baryon distribution amplitudes and the structure function g2(x,Q 2 ) in polarized deep inelastic 
scattering. It turns out that some three-particle evolution equations are completely integrable 
[TTj . since they possess a "hidden" symmetry alias a new quantum number. Complete integra- 
bility is a new and unexpected feature which comes on the top of conformal symmetry and is 
not seen at the classical level from the QCD Lagrangian. The same methods turn out to be 
applicable in a different, Regge kinematics in which case conformal symmetry allows one to 
solve the so-called BFKL equation for arbitrary momentum transfers [T2*] . Miraculously, the 



4 



complete integrability appears here as well [131 E] and allows us to find the Regge intercepts 
of the system of interacting reggeized gluons in the large N c limit. 

Sect. 5 is devoted to applications of conformal symmetry to QCD calculations beyond the 
leading order, in which case change of scaling dimensions for elementary fields and compos- 
ite operators and the running of the QCD coupling have to be taken into account. It was 
conjectured long ago (T^j that the symmetry breaking corrections must be proportional to the 
QCD /5-function and have a rather simple form. The concrete predictions ^H] based on con- 
formal symmetry seemed to be in contradiction with explicit calculations J7j, however, and 
the resolution of this paradox was only found much later ^Hl E3- About the same time a 
pre-QCD conformal prediction, the so-called Crewther relation [20], was verified within per- 
turbative QCD calculations at next-to-next-to-next-to leading order [21]. We explain these 
developments and present the calculation of conformal anomalies that are used to reconstruct 
anomalous dimensions and evolution kernels at next-to-leading order (NLO) [221 122]- Last but 
not least, we demonstrate the predictive power of the conformal operator product expansion 
for the perturbative expansion of light-cone dominated two-photon processes flQ As a 
concrete application, the photon-to-pion transition form factor at next-to-next-to-leading order 
(NNLO) is considered [20]. 

Finally, in Sect. 6 we summarize. The review includes two Appendices: Appendix A is 
devoted to the construction of the orthonormal basis of three-particle conformal operators and 
in Appendix B we sketch the derivation of conformal Ward identities, including ghosts and 
equation of motion terms. 



2 Preliminaries 

2.1 Conformal Group and its Collinear Subgroup 

Among the general coordinate transformations of the 4-dimensional Minkowski space that con- 
serve the interval ds 2 = g liv {x)dx^dx v there are transformations that change only the scale of 
the metric: 

g'^(.x) = u{x)g»v{x) (i) 

and, consequently, preserve the angles and leave the light-cone invariant. All transformations 
belonging to this subclass form, by definition, the conformal group. It is obvious that conformal 
transformations correspond to a generalization of the usual Poincare group, since the Minkowski 
metric is not changed by translations and Lorentz rotations. Examples of specific conformal 
transformations are the dilatation (global scale transformation) and inversion 

x n x 'n = Ax m x n ^ x >» = ?L_ (2) 

X 2 

with real A. Another important example is the so-called special conformal transformation 

a ,a x^ + a^x 2 
x » _> x '» _ (3) 

1 + 2a ■ x + a 2 x 2 ' V ; 
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which corresponds to the sequential inversion, translation by an arbitrary constant vector 
and one more inversion with the same center. 

The full conformal algebra in 4 dimensions includes fifteen generators [27J H] : 



P M (4 translations) 

M MI/ (6 Lorentz rotations) 

D (dilatation) 

K M (4 special conformal transformations) 

and is a generalization of the familiar 10-parameter Lie algebra of the Poincare group generated 
by P M and M^: 

i\P^ P v ] = , i[M af3 , = g afJ Pf3 - g^Pa , 

i[M a p, M^] = g ai Mp u - gp^M au - g^M^ + g^M a ^ . (4) 
The remaining commutation relations that specify the conformal algebra are: 
z[D,P M ] = P M , z[D,K M ] = -K M , 

i[M a/3l K„] = g af J£p - g Pfl K a , ifP^, K„] = -2^D + 2M^ , 

z[D,M^] =i|K AtJ K i ,] = 0. (5) 

The generators act on a generic fundamental field <&(x) (primary field, in the language of 
conformal field theories) with an arbitrary spin as 

z'[P M ,$(x)] = d^(x), 
i[M" v , $(x)] = (x^d v - x v & 1 - TT) $(x) , 
z[D,$(x)] = (x-<9 + £)$(x), 

i[K M , $(x)] = (2x M a; ■ 9 - x 2 d fl + llx^ - 2x v lT v ) $(x) . (6) 

For example, if we consider infinitesimal translations x M — > = x M + e M , the field 
transforms as $'(x) = [1 + e M 5p]$(x), etc. is the generator of spin rotations of the field <3>. 
For scalar, Dirac spinor (quark) and vector (gluon) fields 

TT<i>(x) = 0, E^V = "^V*, = - V mAV > ( 7 ) 

respectively, where <r Ml/ = i[7 M , 7^/2. The parameter £ is called the scaling dimension and it 
specifies the field transformation under the dilatations. In a free theory (i.e., at the classical 
level) the scaling dimension coincides with the canonical dimension £ can which is fixed by the 
requirement that the action of the theory is dimensionless. In the quantum theory I ^ £ can , in 
general, and the difference is called the anomalous dimension. 

An ultra-relativistic particle (quark or gluon) propagates close to the light-cone so that 
one may envisage that separation of transverse and longitudinal coordinates will prove to be 



= 

5 D <$>(x) = 
5^(x) = 
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essential. For further use, we introduce shorthand notations for the projections on the two 
independent light-like vectors. For arbitrary four-vector we define 

A + = , A_ = A„n" , n 2 = n 2 = , n-n = l (8) 

and the metric tensor in the directions orthogonal to the light-cone 

fli = 9iw ~ n i^v - riyfi^ . (9) 

We will also use the notation A± for a generic transverse projection and A± for the vector that 
only has transverse components. For example 

x» = x_n» + x + n» + x? ± , < = g^x v (10) 

and, therefore, x 2 = 2x + s_ — x\. 

Consider a special case of the conformal transformation in Eq. (j^J) with a M being a light-like 
vector = an^. One finds 

x_ — x 7 = — (11) 

so that these transformations map the light-ray in the X- -direction into itself. Together with 
the translations and dilatations along the same direction, x_ — > x_ + c and x_ — > Ax_, they 
form a so-called collinear subgroup of the full conformal group. This subgroup is nothing else 
than the familiar SX(2, R) group [28J. It will play a central role in our analysis. 

In the parton model, hadron states are replaced by a bunch of partons that are collinear, 
that is all move in the same direction, say n M . Whenever this picture applies, one typically only 
needs to consider quantum fields "living" on the light-ray 

$(x) -> $(on) , (12) 

where a is a real number. In what follows we will often use a shorthand notation $(a) = $(an) 
in order not to overload the formulas. We will further assume that the field $ is chosen to be 
an eigenstate of the spin operator £ + _ and so it has fixed projection s on the "+" direction: 

£+_$(«) = s $(a). (13) 

This property is automatically satisfied for leading twist operators [2H! • In the general case 
one should use suitable projection operators to separate different spin components. Several 
examples will be considered below. 

Within these specifications, four-dimensional conformal transformations are reduced to the 
collinear subgroup that generates projective transformations on a line: 

/ aa + b 

a — > a = : , ad — be = 1 , 

ca + a 



$(a) -> &(a) = (ca + d)- 2j <$> ( ^±i^ 

\ca + a J 



(14) 
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where a,b,c,d are real numbers and 



j = {l + a)/2. (15) 

These transformations are generated by the four generators P + , M |_, D and K_ which form 

a collinear subalgebra of the conformal algebra (j3J) and In order to bring the commutation 
relations to the standard form, it is convenient to introduce the following linear combinations 

L+ = U + iL 2 = -iP+ , L_ = Li - iL 2 = (i/2)K_ , 

L = (i/2)(D + M_ + ) , E = (z/2)(D - M_ + ) . (16) 

We obtain 

[L ,L T ]= T L T , [L_,L + ] = -2L , (17) 

which is the algebra of SL{2,M) ~ 0(2,1). Action of the generators defined in Eq. ()16|) on 
quantum fields can easily be derived from Eq. ©, and it can be traded for the algebra of 
differential operators acting on the field coordinates 

[L+,$(a)] = -d a ${a) = L+$(a) , 
[L_,$(a)] = (a 2 d a + 2ja)$(a) = L_$(a) , 

[L ,$(a)] = {ad a + j)${a)=L <S>{a), (18) 

where <9 Q = d/da. They satisfy the similar SL(2) commutation relations: 

[L ,L T ] =±L T , [L.,L+] =2L . (19) 

Note that we use boldface letters for the generators acting on quantum fields to distinguish them 
from the corresponding differential operators acting on the field coordinates. The remaining 
generator E counts the twist t = £ — s of the field $: 

[E,$(a)] = ^- S )$(a). (20) 

It commutes with all Lj and is not relevant for most of our discussions. The definition in (J2*U|) 
corresponds to the so-called collinear twist ( "dimension - spin projection on the plus-direction"), 
which differs from geometric twist ("dimension - spin") [SI] that refers to the full conformal 
group instead of the collinear subgroup. Relations between operators, distribution functions 
etc. of definite collinear twist are well understood by the classification in terms of geometric 
twist since the latter respects full Lorentz symmetry, see e.g. [121 HHH Ell- 

A local field operator <&(a) that has fixed spin projection |T3|) on the light-cone is an 
eigenstate of the quadratic Casimir operator 

£ [Li, [L i? $(«)]] = j(j - l)$(a) = L 2 $(a) (21) 

i=0,l,2 
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with the operator L 2 defined as 



L 2 = Ll + L\ + L\ = Ll - L + L.L+ , [L 2 ,Li] = 0. (22) 

Eqs. (|14|) and (}2~T]) imply that the field $(a) is transformed under the projective transformations 
according to a representation of the SL(2, M) group specified by the parameter j, to which we 
shall refer as the conformal spin of the field. As follows from its definition, Eq. f|15jl . j is positive 
and takes (half)integer values. We will see below that we are dealing with infinite-dimensional 
representations of the collinear conformal group. 

Besides the collinear subgroup just described, one can consider another subgroup corre- 
sponding to transformations of the two-dimensional transverse plane x± = (0, X\, x 2 , 0) intro- 
duced in (11 (Jj) . This "transverse" subgroup involves six generators P^, M^, D and and is 
isomorphic to SL(2, C). It is convenient to introduce complex coordinates 

z = X\ + ix 2 , z = x\ — ix 2 = z* (23) 

in terms of which the group transformations become 

az + b az + b ... 

cz + a cz + a 

where a, b, c, d are complex numbers and ad — be = 1. For the fundamental fields living on the 
transverse plane, $ = $(z, z), the corresponding transformation laws are 

z) - $'(,, z) = (cz + d)- 2h (cz + d)~ 2 H |±1) , (25) 

where h = (£ + A)/2 and h — (£ — A)/2, and A is the helicity of the field defined as S z2 $ = A$. 
This subgroup turns out to be relevant for the QCD description of the high energy scattering 
in the Regge kinematics, see Sect. 4.5. Note that the collinear and the transverse subgroups 
share a common generator of dilations D and, therefore, are not independent. 



2.2 Conformal Towers 

The algebra of the collinear conformal group coincides with the algebra 0(2, 1) of hyperbolic 
rotations which are simply Lorentz transformations in 2+1 dimensional space-time. Similar as 
for usual spin rotations, it is convenient to introduce a complete set of states that in addition to 
the conformal spin j also have fixed conformal spin projection on a certain axis in this internal 
space. The standard way to construct such a conformal basis is in terms of local composite 
operators built of fundamental (primary) fields and their derivatives. 

As easily seen from Eq. ()18|) . the field at the origin of the light-cone $(0) is an eigenstate 
of L with the eigenvalue jo = j an d is annihilated by L_ : 

[L_,$(0)] = 0, [L ,$(0)]=j$(0). (26) 
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It defines, therefore, the so-called highest weight vector on the 5L(2,R) representation space. 
A complete basis of operators on this space can be obtained from $(0) by applying /c-times the 
"raising" operator L + : 



o=0 



O k = [L + , . . . , [L + , [L + , $(0)]]] = [(-d + ) k Ha) 

From the commutation relations it follows that 

[L , O k ) = (k + j)O k , [L + , O k ) = O k+1 , [L_, O k ) 



O = $(0) . (27) 



The primary field operator at an arbitrary position on the light-cone $(a) can formally be 
expanded in a Taylor series over local conformal operators 



E 

k=0 



-a] 



k\ 



-O k , 



(29) 



which can be interpreted as the relation between two different complete sets of states, $(x) for 
arbitrary x and O k for arbitrary k. 

This construction presents the simplest example of what is called the conformal tower. The 
lowest operator in a conformal tower is the highest weight vector in the space of representations. 
Higher operators are obtained by adding total derivatives, each of which adds one unit to the 
conformal spin projection on the "zero" axis. Aim of this Section is to construct conformal 
towers for a general situation of local composite operators built of several fundamental fields 
and (covariant) derivatives. Before doing this, we have to give some more definitions. 

Any local composite operator can be specified by a polynomial that details the structure 
and the number of derivatives acting on the fields. For example, the tower of local operators 
introduced in Eq. (j2*Tj) can be presented as 



O k = [V k {d a )${a)] 



V k (u) 



-u 



o=0 



(30) 



where d a = d/da. The algebra of generators acting on the composite operators (|28|) can 
equivalently be rewritten as algebra of differential operators acting on the space of characteristic 
polynomials. Requiring 



[L Q V(d a )$(a)] 
[L ± V{d a )<S>{a)\ 



o=0 



Q = 



V{d a )[L^{a) 
T(d a )[L T <S>(a) 



o=0 



o=0 



(31) 



one finds the following 'adjoint' representation of the generators acting on this space: 

L V(u) = (ud u + j)V(u), 

L_V(u) = -uV(u), 

L + V{u) = (ud 2 u + 2jd u )V(u), 



(32) 
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where d\ = d 2 /du 2 and in order to maintain the same commutation relations (|19|). we have 
defined L_ as the adjoint to L + , and vice versa. Note that the generators in the adjoint 
representation are more complicated compared to the generators (|18|h in particular L + contains 
the second derivative. The following trick fZ§\ allows us to avoid this problem and proves 
to be very convenient in applications. 

Let us introduce a new variable k instead of u by the following rule 

K n , n = 0,l,2,... (33) 



T(n + 2j) 

In this way a characteristic polynomial V n {u) is mapped onto a polynomial in k: 

V n {u)^V n {K). (34) 

In what follows we refer to Vu{k) as the characteristic polynomial in the adjoint representation. 
The rationale for this name is that the generators of projective transformations become in this 
space 

LoVn(n) = (Kd K + j) V n {K) , 

L_V h (k) = -d K V n (K) , 

L+V n (K) = (K 2 d K + 2jK)V n ( K ), (35) 

i.e., they are given by the same differential operators as the original generators in Eq. (|18|) acting 
on the field coordinate. At this point, the construction of the characteristic polynomial in the 
adjoint representation can look very formal. We will later see, however, that the transformation 
(J2HJ) follows naturally from the conformal transformation properties of correlation functions. 
Let us now consider the product of two fields "living" on the light-cone, 

0(a h a 2 ) = $n(ai) $j 2 (a 2 ) (36) 

with oi\ a 2 - Here we inserted subscripts to indicate that the fields may have different 
conformal spins. Expanding this product at short distances \a\ — a 2 \ — > one encounters the 
set of local composite operators 1 

O n {0)=V n {d 1 ,d 2 )<f> jl {a 1 )<S> h {a 2 ) , (37) 

ai=Q?2=0 

where d\ = d/dai etc., and with V n (u\,U2) being a homogeneous polynomial of degree n. 
Generally speaking the operators O n do not have any simple properties under conformal trans- 
formations. Our aim is, first, to construct a complete basis of local operators that form a 



throughout this review we adopt the following notations: j (or J) always stands for the conformal spin; 
n (or N) denotes the degree of the characteristic polynomial (the number of derivatives of the highest weight 
(conformal) operator), I is the total number of derivatives and and k = I — n (the number of total derivatives) 
denotes the fc-th term climbing up the conformal tower, so that |jo| = j + fe- 
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conformal tower, and, second, to work out the operator product expansion (OPE) in this basis 
— the conformal operator product expansion (COPE). 

The collinear conformal transformations of two-particle operators defined in Eq. (JSTIj) or 
(13 7|) correspond to independent transformations of the fields; the group generators are given, 
therefore, by the sum of one-particle generators 

L a = Li,« + U,a , a = 0,1, 2 (38) 

and the two-particle Casimir operator is 

L 2 = (Lm + L 2 , a ) 2 . (39) 

a=0,l,2 

Our task is to work out the conformal tower of local operators with the given conformal spin and 
the conformal spin projection on a given axis, which is the classical problem of spin summation 
in nonrelativistic quantum mechanics, albeit for a different group, SX(2,R). 

We define a local conformal operator O n by requiring that its transformation properties 
under the collinear conformal subgroup are the same as for the fundamental field ()14j). This is 
equivalent to imposing the following three requirements 

[L 2 ,O n ] = j(j -l)O n , 
[Lo,O n ] = Ui+h+n)O n , 

[L_,Gy = 0. (40) 

The first equation in (}4Tij) states that O n is an eigenstate of the quadratic Casimir operator 
i.e., it has conformal spin j. The second equation requires that O n has a given spin projection 
ji + 32 + n. This condition is satisfied trivially once O n is a homogeneous polynomial in 
derivatives of degree n, cf. Eq. (j3*7|) . Finally, the third equation in (j4*Uj) is that O n is annihilated 
by the "lowering" operator L_, i.e., it defines the highest weight on the SL(2, M.) representation 
space, corresponding to the minimum value of the spin projection j for a given spin j. In other 
words, the conformal operator is defined as the lowest state of an irreducible conformal tower. 
Since L 2 = Lq — L + L+L-, the highest weight condition together with the second equation 
in PU]) make the first equation redundant and imply 

]=]i+]2 + n (41) 

which is the reason that we do not need the second label j for the conformal operator. 

Once the conformal operator (highest weight) is known, there are two possibilities to con- 
struct a complete operator basis. The first possibility is to introduce conformal operators for 
each space-time point O n (0) — > O n (x) with the transformation laws © and (jl4j) . the same as 
for the fundamental fields. The second possibility is to build the conformal tower of operators 
with higher spin projections that can be constructed by a repeated application of the "raising" 
operator L + : 

O n , n+k = [L + , . . . , [L + , [L+, $(0)]]] = (-d + ) k O n , O n ,„ = O n , (42) 
.. ' 
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in full analogy with the one-particle case. Here d + is the total derivative. 

In order to construct two-particle conformal operators explicitly one has to find the corre- 
sponding characteristic polynomials (j3*7|) . This task can be simplified going over to the adjoint 
representation V n {u\,u 2 ) — > V(kx,k 2 ) at the intermediate step [2H]. For two variables, the 
corresponding rule is 

► K? 1 ^ 2 (43) 



r(m + 2ji)r(n 2 + 2j 2 ) 

and in this space 

LoV n (Kl,K2) = {^l9l + K 2 d 2 + ji+ j 2 )'P n {K 1 ,K2) , 

L-V n {Ki,K 2 ) = - (<9i + d 2 ) V n {Kx, K 2 ) , 

L + ^ n («!,«2) = («^ + k\Q 2 + 2j lKl + 2j 2 K 2 ) V n ( Kl , «a) , (44) 

where <9i = d/dnx and <9 2 = d/dn 2 . In particular, if P n («i, ft 2 ) defines the characteristic 
polynomial of a conformal operator (highest weight), it has to satisfy 

L_P„(«i, « 2 ) = - (d x + d 2 ) P„(«i, k 2 ) = , (45) 

which is solved by 

P„(«i, k 2 ) oc (/ca - (46) 

The coefficient of proportionality is irrelevant. Going over to the original variables Ui,u 2 we 
obtain 

„,f^-„ V ™i / r(m + 2ji)r(n 2 + 2j 2 ) \ Ul + u 2 j 



ni+n 2 =n 



(47) 

where Pn ,b \x) are the Jacobi polynomials [35] and we have added the superscripts to indicate 
the conformal spins of the constituent fields. Stated differently, we have proven that the coeffi- 
cients of the Jacobi polynomials are the Clebsch-Gordan coefficients of the collinear conformal 
group. The corresponding local conformal operator is 



2J2-1) 





- d+ 




\d+ 







(48 



where we have restored x as the coordinate. The expression in Eq. (|48j) was first obtained by 
Makeenko using a different method. 

Note that the characteristic polynomials corresponding to the conformal operators with 
different conformal spin are mutually orthogonal with an appropriate weight function: 

(P n |P m > = / dm du 2 6(1 - Ul - u 2 ) ul^ul^F^iut, u 2 ) P£' i2 («i, u 2 ) = U6 mn , (49) 
Jo 
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where M is a normalization constant. This property is related to the fact that the Casimir 
operator L 2 is self-adjoint with respect to the so-called conformal scalar product, see Sec. 2.4.2. 
For the operators (polynomials) that satisfy the highest weight condition, the conformal scalar 
product reduces to Eq. ([49)1 . The orthogonality of characteristic polynomials for the conformal 
operators with different conformal spin follows from hermiticity of L 2 in the same way as the or- 
thogonality of wave functions with different energy follows from hermiticity of the Hamiltonian 
in usual quantum mechanics. 

The above construction can easily be generalized for the products of three and more fields 
in which case characteristic polynomials becomes functions of three and more independent 
variables. For instance, for three fields Eq. (|45|) will be modified to 

L_P„(ki, K2, k 3 ) = - (ft + d 2 + d 3 ) P„(/d, k 2 , k 3 ) = , (50) 
which is solved by any linear combination of terms 

F^ n >( Kl ,K 2 ,K 3 ) OC (« 3 - K 2 )^{k 2 - Kl ) n \ (51) 

where n\ +n 2 = n. Existence of multiple solutions means that the representation is degenerate: 
there exist several ways to sum up the three conformal spins j±, j 2 , j 3 to an overall spin 
j = ji + j 2 + j 3 + n. Accepting Eq. (joTj) as it stands and going over to original variables 
itj, one obtains a particular basis of solutions that are known as Appell's polynomials [35j. A 
disadvantage of using this basis is that the Appell's polynomials corresponding to the different 
choice of ni,n 2 are not mutually orthogonal. A better, orthonormal basis, can be constructed 
if in addition to the total conformal spin j one requires a definite value of the conformal spin 
in a given two-particle channel, see [HOI IHZj arid Appendix A. 

The conformal towers formed by local composite operators are the standard, but not the 
only possible choice of the operator basis. In fact, going over to local operators in coordinate 
space is not always convenient. Another option is to integrate the field $(«) over the light-ray 
with a certain coefficient function 

S = J da<p{a)$(a), (52) 

where it is assumed that the field decreases sufficiently fast at infinity so that the integral is 
well defined, and fix <p(a) by requiring that the nonlocal operator S has desired properties 
under conformal transformations, for arbitrary $. In particular, one can define the conformal 
coherent state as an eigenstate of the "raising" operator 

[L+,5 P ] =ipS p . (53) 

Using Eq. (fTSj) and integrating by parts one finds <f> p (oe) = exp(zpa), so that Eq. becomes 
the usual Fourier transformation and the eigenvalue p can be identified with the momentum 
p + along the light-cone: 

S p = $(p + ) = I dae ip+a $(a) . (54) 
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One can reexpress $(a) in terms of the momentum components (alias coherent states) 

*(«) = / ^e-^°$(p + ). (55) 

Note similarity with Eq. (|29)1 : In the both cases we expand a certain vector in the Hilbert space 
over the complete basis of functions. 

Any given momentum component (eigenstate of L + ) can be decomposed into the 

eigenstates of Lo- Similar as in Eq. (fTHj) the action of generators on the field in the momentum 
representation is reduced to differential operators acting on p + . For the discussion of hadron 
distribution amplitudes in Sect. 3 it is convenient to use polynomial solutions in momentum 
space, which can be constructed as follows. Consider the complete set of functions P 

{P+lj > k) = T(2j + k) {iP+)k+2j ~~ 1 - (56) 
As it is easy to see, the action of the generators (fTHj) in this space yields 

L + (p + \j,k) = {k + 2j){p+\j,k + l), 
L-(p+\j,k) = k(p + \j, k - 1) , 

L {p+\j,k) = -(k+j)ip+\j,k). (57) 

Thus L + acts as the step up and L_ as the step down operator; the state with k = corresponds 
to the lowest \jo\ = j and it is annihilated by L_ so that it defines the highest weight in the 
momentum representation. Any field can formally be expanded in this basis 

$ iG°+) = J2v>k- (P+\j,k), (58) 

k 

where (p^ are operator coefficients that do not depend on momenta. 
2.2.1 Conformal Operators in QCD: Examples 

Our discussion so far was completely general, let us now consider some examples. In many 
QCD applications to hard processes one encounters a nonlocal operator built of the quark and 
the antiquark at light-like separation 

Q jU (a 1 , a 2 ) = ^("i)7m[«i, a 2 }^(a 2 ) , (59) 

where 



[at, a 2 ] = Pexp 



ig / dtA+(t) 



a ■> 



(60) 



is the Wilson line. Expanding Qn(ai,a 2 ) at short distances gives rise to local operators built 

of the quark and the antiquark field and covariant derivatives, of the type ijj(D+) ni/ y fJ ,(D+) n2 ilJ ■ 
Our task is to find the corresponding conformal operators. 
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The starting observation is that we cannot identify the quark field operator ijj directly with 
the fundamental field $, but first have to separate contributions with different spin projections. 
To this end we introduce the spin projection operators 

1 „ 1 



n+ = -7_ 7+ , n_ = - 7+7 _ , n + + n_ 

and define "plus" and "minus" components of the quark field as 



(61) 



(62) 



which is the same decomposition as is made in light-cone quantization, although the motivation 
is in our case different. Using Eq. (0) one easily finds that ip + and ■?/>_ correspond to the spin 
projections (fT3|) s = +1/2 and s = —1/2, respectively. Remembering that the canonical 
dimension of a spinor field equals 1$ = 3/2, one finds for ip + (ip-) the values of the conformal 
spin (JT5Jl j — 1 (J — 1/2) and twist (}2T)j) t — 1 (t — 2), respectively. Now we see that different 
Lorentz projections of the operator Q M («i, a^) correspond to taking into account different quark 
field components and have different properties under conformal transformations: 



twist-2: 


Q+ 


= ^+7+^4- 


twist-3: 


Qx 




twist-4: 


Q- 





+ Q 



1/2,1 



Ql/2,1/2 



(63) 



where the superscripts stand for the conformal spins of the quark and the antiquark. The 
corresponding local conformal operators are 



jl/2,1/2, 



X) 



X) 



x) 



(id + ) 
(id + ) 



^{x) 1+ CT [D+ /d + )i>(x) 
^(x) 7+ 7±7-^ 1 ' 0) (B+ /d+)il>{x) 
^(*)7-Cy 2 (d+ /d + )l>(x) 



(64) 



where D+=D- 



D+, d + 



-D- 



D+ and we replaced Jacobi polynomials with two equal 



superscripts by the Gegenbauer polynomials P, 



(i,i) 



a 



3/2 p (0,0) 



Cl' 2 35 . We also added 



J/2,1 



are 



an overall factor i n which is convenient for applications. The conformal operators 
given by the same expression as Qn 1 ^ 2 with an obvious substitution Pn — > Pn ' 1 ^ ■ 

The similar analysis can be carried out for gluons. In this case one finds that the G + ± 
component of the gluon field strength tensor corresponds to s = 1 and therefore a "legitimate" 
fundamental field with the conformal spin j = |(2 + l) = 3/2 and twist t = 2 — 1 = 1, G±± 
and G + - both have s = 0, j = 1 and t = 2, and finally has s = —1, j = 1/2 and twist 
t = 3. Local conformal operators built of two gluon fields of leading twist are, therefore 



G 3/2,3/2 ( 



(id. 



G a + Ax)C 5 J 2 (d+ /d + ) G°, 



(65) 



Depending on the particular combination of transverse indices there exist three towers of gluon 
operators corresponding to the two-gluon states with different helicity. 
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2.3 The Conformal Operator Product Expansion 

It has been known for a long time that the conformal symmetry provides powerful constraints 
for the Wilson coefficients that appear in the operator product expansion (OPE) of currents 
at light-like distances. Using conformal operators as a basis suggests itself and was studied 
in a number of papers, starting in the seventies with the work of Ferrara, Grillo and Gatto 

HOI HU H21 EHj • It has been employed in QCD to find the eigenf unctions of the flavor non- 
singlet evolution kernel beyond the leading logarithms (THj and to relate the Wilson coefficients 
and hard-scattering amplitudes that appear in a number of exclusive two-photon processes 
jlH IIHI EZj , to those known from deeply inelastic scattering [2U 123 ■ 

Consider the operator product A(x)B(0) of two local conformal operators with twist tj\ 
(is) and spin projection on the "+" direction 5,4 (sb), respectively. We want to expand this 
product in the light-cone limit x+,x± — > and x_ fixed, i.e., x 2 — > 0, over a given tower of 
local conformal operators and their derivatives 0^i fe which, as we assume, forms a complete 
operator basis. For simplicity, we will restrict ourselves to the leading twist contributions and 
also neglect (the most singular) contribution of the unit operator which is most often irrelevant. 
In a free field theory, for dimensional reasons, this expansion looks as follows: 



A{x)B(0) = J21Z C ^[^) ' x - +k+A KnlM + higher twist . (66) 

71=0 fc=0 ^ ' 

Here A = si + S2 — sa — sb, where s\ and S2 are the spin projections of the constituent 
fields in the local operators ©^+ fe and C n ^ are the corresponding Wilson coefficients. The 
singular behavior for x 2 — > is governed by the twist of the operators: tA + ts — t, where 
t n = £ n — n — si — S2 = £1 + £2 — si — S2 is the twist of the conformal operators O^ 1 '- 72 . We 
remind that the twist of all operators in the conformal tower and thus the 1/x 2 behavior of all 
contributions are the same. In the interacting theory the scaling dimensions of the operators 
£ n will be modified by the anomalous dimensions and this universal behavior no longer holds. 

The power behavior in x_ is determined by the angular momentum conservation. This can 
be established in the following way. We have assumed that both A and B have definite spin 
projection 

£+-A(0) = s A A(0) , £+_B(0) = s B B(0) (67) 
on the "+" direction. On the other hand 

»* [M_h Oji*U(0)] = (si + s 2 +n + k)& n f +k (0), 

i[M-+,A(x)B(0)] = [x-n-d-x + n-d + s A + s B ]A(x)B(0). (68) 

Insert the OPE (jHEJ) in the second equation. On the l.h.s., evaluate the commutators 
i [M_ + , ®^ifc(0)l , which gives s\ + $2 + n + k for a generic term, see Eq. (jHEJ)- On the r.h.s., 
take the derivatives which gives a factor SA + ss + A + n + k. Comparing the both expressions 
we conclude that A = si + S2 — sa — sb, as anticipated in Eq. (fBT)]). 
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The conformal symmetry allows us to make one further step and obtain Wilson coefficients 
C nt k with k = 1, 2, . . . corresponding to operators with total derivatives in terms of the coefficient 
of the highest weight operator C n ^=o = C n . One way to see this is to act with the generator 
L + on both sides of the OPE and compare the two results in the leading twist sector, see 
e.g. [10]. This is completely analogous to the above derivation of the constraint arising from 
Lorentz invariance. Consider an infinitesimal special conformal transformation 

[L_ , A{x)B(0)] = |ar_ (2j A + x ■ d x ) A{x) - X -x 2 n- d x A{x) } B(0) + ■■■ , (69) 

Here j A = (s^ + ^U)/2 is the conformal spin of the operator A and ellipsis stand for higher 
twist contributions arising from the spin rotation, i.e., proportional to x±Y>_±A(x). Inserting 
the OPE (ffiljj) in Eq. (JfiTJj) and taking into account that 

[L-> = ~k(k + 2j„ - 1)0^^(0) with j n = j x +j 2 + n, (70) 
cf. Eq. (|2*%|) . results in a first-order recurrence relation for the Wilson coefficients 

n _ U-jB+Jn + k 

Cn ' k+1 --(k + l)(k + 2 3n ) Cn > k > (71) 



which is solved by 



k\{2] n ) k 



where (a)k = T(a + k)/T(a) denotes the Pochhammer symbol. Inserting this solution in the 
OPE (jnni) allows us to sum contributions of operators containing total derivatives and yields 
the representation: 

00 / 1 \ tA+t z- tn n+Sl+S2 - SA -s B 

^ / BUA-jB+3n,3B-jA+3n) 

x [ duu SjrjB *' r - 1 (l -u) j ^ A ^ n - l O%' j2 (ux„) , (73) 
Jo 

where B(a, b) = T(a)T(b)/T(a+b) is the Euler Beta function. The total derivatives get absorbed 
in a translation of the conformal operator along the light-cone with an appropriate weight 
function. The case of an interacting theory will be considered in detail in Sect. 5. We will 
find that the scaling dimensions of operators are modified by quantum corrections so that the 
twist and the conformal spin have to be modified accordingly: j n = (£ n + s\ + S2 + n)/2 and 
t n = (4 -S1-S2- n). 

To give an example, consider the operator product expansion of two electromagnetic currents 
J^ 1 = where is the electric charge and the sum runs over (light) quark flavors 

i = {u,d, s}. At tree level only the transverse components of the currents are of interest, 
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which have the spin projections sj = 0. Since the current has dimension three, it possesses 
twist tj = 3. The relevant operator basis is provided by the quark operators Q*' 1 (|64|) with 
conformal spin j n = (£ n + 1 + n)/2 = 2 + n and twist t n = (£ n — 1 — n)/2 = 2. Moreover, A = 1 
and, therefore, we find from Eq. ()73|) 



6-t 



J ± (x)J ± (0)^y2C n (-) — / duu^il-uY^q^iuxJ). (74) 

^ \x 2 J B(j n ,j„) y 

Note that the remaining coefficients C n are those that appear in the forward kinematics: 

{P\j\x)j\Q)\P)^C n (-^\ 2 (-.x„)" +1 (P|Q^(0)|P). (75) 

n=0 VX ' 

Up to a trivial difference in normalization that arises from the definition of conformal operators 
they coincide with the Wilson coefficients known from deep inelastic scattering. 

As a further example, consider the short- distance expansion of a bilocal light-ray operator 
<2+ in Eq. (JH9*|) . In this case we have to set sa = sb = «i = s 2 = 1/2 so that A = and 
£ A = £ B = £ 1 = £ 2 = 3/2. We find 

Q + {a 1 ,a 2 )~y2C n {-i) n {a 1 -a 2 ) n / duu n+ \l - u) n ^^t\ua x + (1 - u)a 2 ) , (76) 

n=0 ^ 

where we set C n = C n /B(n + 2,n + 2). The simplest way to evaluate the coefficients C n is to 
form forward matrix elements. The result reads C n = 2(2n + 3)/(n + 1)!, in agreement with 

Ref. [SHi cni- 

Before closing this section, we want to mention an alternative approach suggested in |3*5] . 
which is based on the introduction of the complete basis of multi-particle coherent states. For 
instance, for two particles, one can define the coherent state 

S hP = JdadPhjfaflQfi Q(a + 0)) Qfc-J 9 )) ( 77 ) 
by imposing the two conditions: 

[L 2 > = 3 U - 1 ) <%p > ( L + , £j,p] = *P , ( 78 ) 

c.f. Eqs. (fH2"|) - (f5*4"J) . These conditions can be translated to a second-order differential equation 
for the function (f>j )P (a, /3), which can be solved by separation of variables and turns out to 
be related to the Bessel equation, see [3S] for the details. An important difference with the 
expansion in local operators is that the eigenvalue j is not quantized in this case, so that we 
are tacitly going over to a different series of representations of the SL(2, R) group |28j . For the 
leading twist quark-antiquark operator (JH9*j) one obtains, in particular 



hp = J^dae ipa J™^ 2 dpH]_ 1/2 (Pp) $ Q(« + 0)J Q(« - , (79) 
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where f/J^^z) = Jj-i/2(z) —iNj_i/ 2 (z) is the Bessel function. The bilocal field product on the 
light-cone can be expanded in contributions of coherent states making use of the Kantorovich- 
Lebedev transformation 

-r rioo 

(35{(3 -/?)= / jdj Jtf) H]{(3 ! ) . (80) 

^ J -ioo 

One obtains, for quarks: 

$ Q(« + /?)) 1^ (~(« -®)=-f% t^ dj U ~ 1/2)(3 ' 3/2 J i-V 2 (P/5) (2+)i,P • 

1/2 " 100 (81) 

Similar representations exist for gluon operators 

Note that a bilocal quark operator on the light cone involves conformal coherent states 
with an imaginary conformal spin for which the weight function if?^^ in Eq. (|79p decreases 
exponentially at large quark- ant iquark separations. Such states can be interpreted as quasi- 
particles (wave packets) moving along the "plus" direction on the light-cone with momentum 
p and having "mass" j in some internal space. We will argue in what follows that operators 
with different conformal spin do not mix under renormalization in one loop. This implies that 
the shape of the conformal wave packets remains the same at all scales, to this accuracy. 

2.4 Conformal Invariance of a Field Theory and its Violation 

Up to this point our discussion has been entirely about the group-theoretical aspects of the 
conformal symmetry and its collinear subgroup: we discussed the field transformations and 
learned how to construct composite objects with simple properties, which is the task of the 
representation theory, in mathematical sense. We are now going to address the field-theoretical 
aspects of conformal symmetry and the first question to ask is: Under which conditions the 
field theory is conformally invariant at the classical level? We demonstrate that the necessary 
and the sufficient condition for conformal invariance is the possibility to introduce the trace- 
less, symmetric and divergenceless energy-momentum tensor. As well known, symmetries of a 
classical field theory can be broken by quantum corrections, in which case one speaks about 
quantum anomalies. This indeed turns out to be the case for QCD. We introduce the confor- 
mal anomaly and show why it is related to the Gell-Mann-Low /3-function. A powerful tool 
to obtain constraints on the correlation functions that follow from the conformal invariance is 
provided by conformal Ward identities (CWIs) jlHl EH EDI EE] • We will see that the dilatation 
Ward identity is nothing else but the familiar Callan-Symanzik equation [22] while the Ward 
identity corresponding to special conformal transformation is the genuine new feature of a con- 
formal theory and has important consequences. Discussion of the conformal invariance in a 
gauge theory like QCD involves subtleties because the symmetry can be broken spuriously by 
the gauge fixing procedure. These difficulties are of technical nature and are discussed in detail 
in Appendix B, where we present a complete set of QCD conformal Ward identities including 
the contributions of BRST exact and equation of motion (EOM) operators. 
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2.4.1 Energy-Momentum Tensor 

According to the Noether theorem, the invariance of a four-dimensional theory under confor- 
mal transformations implies the existence of fifteen conserved currents. The corresponding 
conserved charges are the fifteen generators of the conformal group: P M , M M!/ , D, and K M . We 
will show that these generators can be expressed in terms of the so-called improved energy- 
momentum tensor Q til/ (x) which is divergenceless, symmetric and traceless: 

e^(z) = e^(x) , ^e^(x) = o , g^e^x) = o . (82) 

We will also see that the possibility to define the energy-momentum tensor with such properties 
is closely related to the conformal symmetry of the Lagrangian. 

For simplicity let us consider a generic Lagrangian £($(x), d^(x)) that depends on a single 
fundamental field <&(x) and its first derivative d^(x). Invariance of the action under a global 
continuous transformation of the field, $ — > $ + implies that its variation under the similar 
local infinitesimal transformation characterized by the parameter e a (x) takes the form 



5 F S — I d x 5 F C — I d x 



d(d^(x)) v v " <9$ 



d 4 xJZ(*)^r- (83) 



On the other hand, any variation of the action must vanish at the stationary point, i.e., as soon 
as the fields satisfy the Euler-Lagrange equations of motion (EOM) 



9C TT 9C 



= ^ . n ^) = ET^TZnn • (8 4 ) 



As a consequence, the current J£(x) must be conserved, d^J£(x) = and the corresponding 
conserved charge is given by Q a = J d?x J&(x), dQ a /dt = 0. 

In this way, the invariance of the action under space and time translations leads to conserva- 
tion of momentum and energy, respectively. Substituting 5 e $(x) = s a (x)5p$(x) = s a (x)d a $(x) 
(see Eq. (0)) in Eq. (jHHJ) one identifies the corresponding conserved current Jp a as the canonical 
energy-momentum tensor 



J£ a = T" a = n»d a $-gi* a C, d^ a E ^ M 0, (85) 



where the superscript eom indicates that the fields satisfy the equations of motion (|84|). The 
corresponding conserved charge is the operator of the four-momentum 



P a = J d d xT\(x). (86) 

The invariance of the action under Lorentz rotations, 5 £ §(x) = e a pS'^^(x), (see Eq. (J0J)) implies 
further six conserved currents which are contained in the angular momentum tensor: 

Jm,oP = BoP^ix) - x p T\(x) - IF(x)E Q/3 $(x) . (87) 
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It is given in terms of the canonical energy- momentum tensor ([850 . plus the additional term 
IPX a ^<3> that arises from the transformation of the spin degrees of freedom. The tensor T^(x) 
has no definite symmetry with respect to its Lorentz indices, in general. However, conservation 
of the angular momentum (jH7[) implies that the antisymmetric part T a p(x) — T a p(x) is given 
by a total divergence, since 

•Vm,^ = T ^{x) - Tp a (x) - <9 M (IP(:r)£ a/3 $(x)) E = M . (88) 

This suggests to define the symmetric, so-called Belinfante, form of the energy-momentum 
tensor [5*3] : 

t% = + V (n A s Ml/ $ - n M s Ai ,$ - rL,£ v $) . (89) 

The symmetry T® u = follows readily from Eq. (JSHJ). Since the expression in the parenthe- 
sis in the r.h.s. of Eq. ()89[) is antisymmetric with respect to A and /i, the redefined energy- 
momentum tensor is conserved, d^T® = 0, so that we can equally well use it instead of the 
original Noether current T^ v as the conserved current associated with the translation invariance. 
Note that the operator of four-momentum ([86)1 is not affected by this redefinition. We can also 
use the Belinfante form of the energy-momentum tensor to redefine the angular momentum 
tensor density as 

_ r T B,^ _ rpB,» a Jfl EOM , . 

J M,af3 ~ ^a 1 p • L P ± a ' u ^ J M,a(3 ~ u ' W U J 

The corresponding conserved charges define the operator of angular momentum 

M a/3 = j d 3 xJ^ a/3 (x). (91) 

We stress that the possibility to define the symmetric energy-momentum tensor reflects the 
Lorentz-invariance of the theory. 

Next, consider scale transformations x^ — > Ax M and $(x) — > A * $(Ax) with being 
the scaling dimension of the field. The action of a four- dimensional theory is scale invariant 
provided that C — > A _4 £, i.e., the scaling dimension of the Lagrangian is equal to its canonical 
dimension (four). For the usual Lagrangians that are given by polynomials in the fields and 
their derivatives, this condition translates to the requirement that all couplings have to be 
dimensionless and the scaling dimensions of the fields are uniquely fixed and are equal to their 
canonical values. The conserved current corresponding to the scale transformations can be 
constructed by calculating the variation of the action ([%3~j) under dilatations <5 £ $ — > e(x)5d&(x) 
(see Eq. ©) 



5 D S = J d 4 xe(x) [-d^D" + A D (x)] 



(92) 



The Noether current = is called the dilatation current. It is written in terms of the 
canonical energy-momentum tensor 

D" = x a T^ a + 4IP$ , (93) 



22 



and A D (x) is denned as 



dC dC 

Ad = iu + 1} W^j d "* + u M®- d£ " (94) 

where d — 4 is the space-time dimension. Invariance of the action under dilatations implies 
that Ad(x) is a total divergence. If the Lagrangian £($, <9$) is a polynomial in <9 M $ and $, the 
only possibility is A D (x) = by itself, or, equivalently, £($/A £ *, <9<I>/A £ * +1 ) = £(<3>, <9$)/A 4 , as 
anticipated. 

Last but not least, examining the transformation of the action under the special conformal 
transformations <5 e $ = e a (x)5'^^(x), see Eq. one gets 

S K S = J d A xe a {x) [-d^K afl + 2x a A D (x) + A a K {x)} , (95) 

where 

K afM = {2x a x u - g a v x 2 )T^ + 2x v W{Ug va - S ra )$ , (96) 
A D (x) is defined Eq. (JM|) and 

A£ (x) = 2U V (4 a " + S a ") $ . (97) 

Invariance of the action (J95j) under the special conformal transformations requires that Ad(x) = 
and, in addition, A^(x) has to be a total divergence 

n„{l 9£ r + E° a ')<l> = Pv^x) , (98) 

where a pv is some tensor. Under these conditions the conserved current corresponding to 
special conformal transformations is given by J^ a = ~ 2<7 a M . Explicit calculation shows 
that for physical spin— 1/2 and spin— 1 fields with the canonical dimension i = 3/2 and i = 
1, respectively, the field a ap {x) is identically zero. In other words, in this case invariance 
under special conformal transformations proves to be a consequence of the dilatation (scale) 
invariance. 

Using the condition (J98|) . one can again redefine the energy- momentum tensor to make it 
traceless, and so the expressions for the currents and J^- a are significantly simplified. One 
defines 

= T l + \d x d»X Xpilv , (99) 

where T^ v is the Belinfante tensor and X\ ptiv is chosen in such a way that Q pv remains sym- 
metric and divergenceless, and, in addition, g^ u d x d p X\ PIJiV = 2d a d x a a \ = —2g^ u T^ v . The latter 
condition ensures that ® pi/ is traceless, as it is easy to check using Eqs. (f%KJ) and (jHEJ)- Since 
the difference Q pi/ — T^ v is given by the total divergence, the angular momentum tensor density 
()90|) can be redefined in terms of Q pu with no effect on the conserved charges P p and M M1/ . 
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Moreover, the currents associated with the dilatation and special conformal transformations 
take a remarkably simple form: 

J£ = x v &" , J^ a = (2x u x a - x 2 g ua ) 9^ , (100) 

leading to the conserved charges 



D = / d 6 x r D (x) , K a = J d 6 xJ£ >a . (101) 

From the representation in Eq. (jlOOj) it is obvious that the dilatation and special conformal cur- 
rents are conserved if and only if the improved energy-momentum tensor is divergenceless 
and traceless. 

In a quantum field theory with a running coupling constant these two conditions cannot 
be preserved simultaneously. Due to the appearance of ultraviolet (UV) divergencies, one has 
to define the theory with an UV cutoff. Any regularization that preserves Poincare covariance 
breaks inevitably the dilatation and consequently the special conformal symmetry. This im- 
plies that the renormalized improved energy momentum tensor remains divergenceless and 
symmetric but its trace does not vanish J^l- Breaking of a classical symmetry due to quantum 
effects is called an anomaly, the conformal anomaly for the case at hand. 

The loss of the scale invariance in the quantum theory is rather obvious, because a dimen- 
sionfull parameter appears that plays the role of the ultraviolet cutoff. Consider QCD with 
massless quarks 

£ QC B = -\G%G a ^ + ^P^ } (102) 

where = d^-ig ( -°H a A a (x) is the covariant derivative, and divide the quark and gluons fields 
in "fast" and "slow", with frequencies above and below a certain scale fi, respectively. The 
renormalized QCD action can be thought of as an effective action, integrating out contributions 
of fields with large frequencies (see e.g., [53]) in the background "slow" fields. In the one- loop 
approximation the effective action can be written, schematically, as 

= -\j ^ X ~^ lnMVfl2 ) K^tlow (*)+-> ( 103 ) 

where M 2 is the UV cutoff, /3 = (U/3)N e - (2/3)N f is the first coefficient in the Gell-Mann- 
Low function and the ellipses stand for the quark terms. For the purpose of the argument we 
rescaled temporarily the gluon field g^A^ — > A^ so that the coupling only enters as an overall 
factor in front of G a ^ v G aflv . Making the scale transformation x M — > Xx 11 , Ap(x) — > XA^Xx), 
ip(x) — > A 3 / 2 ^(Ax) and [i — > n/X with the fixed cutoff results in the variation of the effective 
action 

SS = -^L/3 lnA J d 4 xGl u G a ^{x) , (104) 
which implies that (cf. Eq. ([83))) 



Wx) = -^G%G^{x) . (105) 
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Restoring the standard normalization of the gluon field and introducing a s = g 2 / (4?r) we obtain 



X) 



(3 a s 



in 



This is the desired result. The full answer reads jTo] 



where 



d 



A D (x] 



EOM 



2g 



G«G^{x) 



P(g)/g 



(106) 



(107) 



(108) 



and this is correct in operator sense, for an insertion in arbitrary Green functions. Note that 
Ad(x), which defines the variation of the action (|92j) . is the anomalous term in the trace of © M M 
and both differ by extra EOM terms. 

In practical calculations one usually employs the dimensional regularization extending the 
number of space-time dimensions to d = 4 — 2e with e > and subtracting poles in e within the 
minimal subtraction (MS) scheme. This regularization respects gauge invariance and also it is 
explicitly Lorentz covariant so that both the energy-momentum and the angular momentum 
tensor remain conserved. The modified energy momentum tensor in QCD is equal to 



g^£qcD - Gl x G a uX + - Dv ip — ip Dv j„if) + ^ v) 



+ 



(109) 



where we do not show the gauge-fixing terms as well as the terms depending on Faddeev- Popov 
ghost fields [T5|. One can verify that this expression satisfies the defining relations (J82j) . so far 
as the EOM are employed. In particular, taking the trace one finds 



g^e^(x) 



(d-^G^G a ^(x) + 



(110) 



which would vanish for d — > 4 if the operator on the r.h.s. were finite. The renormalization adds 
Z— factors which have poles in (d — 4) and hence a non vanishing contribution to Eq. (jllOj) in 
the limit d — > 4, reproducing the result in Eq. ()107J) . 

There are several subtleties, however. First, since the regularized action (in dimension d) 
has to remain dimensionless, the quark and the gluon fields must have the canonical (mass) 
dimensions £^ an = 3/2 — e, £ c j± n = 1 — e, respectively. These dimensions are fixed by the kinetic 
terms in the QCD Lagrangian and this implies, e.g., from the g^tpfitp term, that the bare 
coupling has to carry the dimension [g^] = — e, g^ oc fi~ e where fi is a certain mass scale. 

Note that the invariance with respect to special conformal transformations is also necessar- 
ily broken and the corresponding anomaly potentially contains a new contribution, A^(x), in 
addition to Ad (re), see Eq. (|95j) . To avoid extra terms and have a simple pattern of conformal 
symmetry breaking governed by A^ alone, the condition (|9"%j) should be satisfied in the regu- 
larized theory. This forces us to set the scaling dimensions of physical fields in d dimensions 
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to be the same as in four dimensions and, therefore, different from their canonical dimensions: 
£</, = 3/2 and I a = 1- In other words, the anomalous dimensions 7 = I — £ can acquire a term in 
e, similar as it happens in the /9-function. 

Second, there are subtleties related to the gauge fixing procedure. Fixing the gauge may 
lead to spurious symmetry breaking already in the classical theory. For example, the Schwinger 
gauge x^A^ = respects classical scale invariance and inversion, but breaks the translational 
symmetry which is then restored in gauge invariant quantities. In turn, the Lorentz invariant 
gauge condition d^A^ = breaks special conformal invariance. Hence, one may require an 
additional gauge transformation in order to have invariance under the action of group generators 
in the simple form © [SH]. Alternatively, without this additional gauge transformation one 
finds that the gauge fixing and ghost terms imply the non vanishing fields a au , defined by Eq. 
()98|) . Since the ghost field u plays the role of the gauge parameter, it is convenient to choose 
£uj = £™ n = and in this case it does not contribute to a av . The scaling dimension of the 
antighost field Co can be chosen as — ^f n = d — 2. Within these settings, the complete result 
for the trace anomaly of the energy-momentum tensor reads 

A(x) = A D (x) - (d - 2)d x O BX (x) , (111) 

where Ad(x) = {(3 /2g)[G a ^ u G a,J,v }{x) + . . . is defined in terms of renormalized operator insertions 
in Eq. (|B.18|) . The additional term here is a total divergence of the operator Ob\(x), given in 
Eq. ()B.15[) . It stems from the special conformal variation of the gauge fixing and ghost terms, 
= 2(d — 2)0%, and is BRST exact. Thus, it does not contribute to the physical QCD 
sector. The expression (jlll|) simultaneously determines the (global) dilatation and the special 
conformal variations of the action 

5 D S = eJ d d xA{x), S K S = e a J d d x2x a A{x) . (112) 
2.4.2 Conformal Ward identities 

Symmetries in quantum field theory imply constraints for Green functions, which are called 
Ward identities. A generic Green function is given as the vacuum expectation of the time 
ordered product of Heisenberg operators and can be evaluated from the following functional 
integral: 

(X N ) = (0|T$(xi)...$(xjv)|0) = N~ l j V^^{x 1 )...^{x N )e iS ^ ) . (113) 

The Ward identities follow from the invariance of the functional integral under the change of 
variables $(x) —> $'(x) = $(x) + 6$(x). One obtains [qTH IqT] 2 

= (0|T<y$(zi) . . . $(x N )\0) + ... + (0|T$(xi) . . . 5<S>(x N )\0) + (0|mS$(xi) . . . $(xjv)|0) . 

(114) 



2 For the situation that we are considering the Jacobi determinant det[2?<I>' is a C-number and, thus, 
irrelevant. 
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For example, the common wisdom that Green functions in a Lorentz-invariant theory only 
depend on the set of invariants x 2 k = (x, — x^) 2 is in fact the consequence of the Ward identities 
for translations and Lorentz rotations, which leave the action invariant, of course. For the 
dilatations one obtains the Ward identity 

- r 

(4 an + xt ■ d t ) (0|T$(an) . . . $(x N )\0) = -i d 4 x{Q\TA D (x) $( Xl ) . . . $(x N )\0), (115) 

where d{ = d/dxi, which tells that scale transformations of the fields have to be compensated 
by the variation of the action. The latter is given by the zero-momentum insertion of the trace 
anomaly of the energy-momentum tensor including EOM terms that are proportional to the 
canonical dimension and count the number of the fields. The Ward identity (|115|) is nothing else 
as the Callan-Symanzik renormalization group equation [221, but in order to see this we have 
to rewrite it in a somewhat different form. First, the Lorentz invariance and simple dimension 
counting suggest that the renormalized Green function in four dimensions can only be of the 
form 

(X N )(x 1 ,...,x N ;fi,g(fi)) = ii m * a G(x 2 k n 2 ;g(ii)) . (116) 

As a consequence, the differential operator on the l.h.s. of Eq. ()115|) can be equivalently written 

as 

N „ 

V (£| an + x t ■ di) (X N ) = ^(Xn). (117) 

Next, consider the operator insertion of the dilatation anomaly <9 M Jp — Ad on the r.h.s. of 
Eq. ()115j) . From our discussion of the one-loop QCD example in Eq. (|1U3|) it is clear that the 
same result (|106J) can be obtained by taking a derivative of the effective Lagrangian over the 
UV cutoff, 

A D (x) = -M-^£ cS (x), (118) 

and it is easy to understand that this relation is quite general, provided the theory does not 
include explicit dimensionfull parameters like masses. It follows that 

i J d 4 x A D {x) e iS ^ ] = -M-^e iS ^ m (119) 

and the derivative can be taken out of the functional integral since the fields $ are just 

integration variables and do not have any M-dependence. We obtain, therefore 

% j 'd 4 x (Q\TA D (x)$( Xl ) . . .$(x N )\0) = -M^(X N ). (120) 

On the other hand, renormalizability of the theory implies that the Green function (X^) de- 
pends on the cutoff through Z$-factors that take into account the field renormalization 

$ unr (x) = y%$(x) , 7 * = ^ In y/Zi = ~M-^ In Jzl (121) 
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and through the coupling constant. In other words 




Combining this with Eqs. (jll5j) and (|117|1 . we indeed obtain the Callan-Symanzik equation: 

(X N ) = 0. (123) 

Let us assume for a moment that the (3 function has a fixed point g* = g(fi*), such that 
f3(g*) =0. In this case by the shift of the scaling dimension 

£ = £ c ™ £(g*) = £ c " n + j(g*) (124) 

one can bring the dilatation Ward identity (J115)) to the simple form that it has in a free theory: 

N 

( W) + Xi ■ di) (0|T$(xi) . . . $(xjv) |0) = . (125) 

t=i 

In this sense the dilation symmetry is only broken by the coupling renormalization. Eq. ()125j) 
is a strong constraint which fixes the structure of two-point Green functions at the fixed point 
f3(g*) = up to an overall normalization factor N(g*), e.g., for two scalar fields 

(010(^)0(^)10) = N 2 (g*)(^)- W) (i ' ' (126) 

~ x 2 ) J 

where the additional explicit scale dependence has been introduced for dimensional reasons. 
For particles with spin and fixed spin projection s on the light-cone ()13j) this is generalized to 

\2{xi -x 2 )-(x 1 -x 2 )+J \{x 1 -x 2 )-J 

(127) 

where it is assumed that (x\ — x 2 )± = 0. It follows that the dependence of the correlation 
function on the light-cone separation {x\ — x 2 )_ (in the light-cone limit (x\ — x 2 ) + — > 0, 
(xi — x 2 )- = const) is determined by the conformal spin j = (£(g*) + s)/2, see Eq. (|15|). 

As explained in the previous Section, a field theory that is invariant under dilatations is 
usually also invariant under the inversion of coordinates, x M — > x^/x 2 . This implies the existence 
of more Ward identities that are generated by special conformal transformations: 

N 

J2 (2<(4 an + Xi ■ di) - 2£*X - x 2 d») (X N ) = -i d i x2x» (0|TA D (x) $(x x ) . . . $(x N )\0) . 
i=i 

(128) 



9 of s 9 
V—+f3(g) 



dfj, 



dg 



N 
i=l 
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The expression on the r.h.s. has, up to factor 2x M and subtleties with the gauge fixing procedure, 
the same form as in Eq. (J 122)) and can be written as 



N 



i / d 4 x2x^(0|TA D (x)$(x 1 )...$(x JV )|0) = ^ 2xf7 0i (O|T$(x 1 ) . . . $(x N )\Q) (129) 
"* i=i 

+ —J d 4 x2x tM (0\Tg-^-C(x)^(x 1 ) . . . $(x N )\0) . 

To be precise, we mention that the renormalization of dC(x)/dg might require total divergence 
counterterms that do not contribute for dilatation. At the fixed point /3(g*) = one obtains 
the same form of the special conformal Ward identities as in a free theory 

N 

^ (2x^U(g*) + Xl ■ di) - 2W v x V i ~ x 2 M) (0|T$(x!) . . . $(a^)|0) = , (130) 
i=i 

but with the shifted scaling dimension (|124j) . The power of conformal symmetry is to a large 
extent due to this identity If combined with the dilatation Ward identity ()125|) and Lorentz 
invariance, Eq. ()130|) is solved by the requirement that the Green functions only depend on 
the so-called anharmonic ratios x 2 k x 2 m / (x^x^) which are invariant under the inversion. Most 
remarkably, conformal invariance fixes the three-point Green functions up to an overall nor- 
malization factor [H7j, e.g., for three scalar fields 



(o|0(z a )0(z2)0(*3)|o) = W)0O~ W) 



1 X 3[« can + 7 (9*)]/2 



(xi - x 2 ) 2 (x 1 - x 3 ) 2 (x 2 - X 3 ) 2 

(131) 

Everything said in this Section remains true if instead of the Green functions written in terms 
of fundamental fields $(x) (|113|) one considers correlation functions of arbitrary local conformal 
operators provided that their transformation properties under the conformal group, albeit with 
different spin and scaling dimensions, are preserved on quantum level. We will discuss this point 
in detail later in Sect. 5.2. An important consequence of the conformal symmetry, encoded in 
the special conformal Ward identity ()130|) . is that the correlation function of two conformal 
operators vanishes in the conformal invariant theory, unless they have the same conformal spin: 

(0|TOi(x)O 2 (y)|0> = 8 jlda F(x - y) , (132) 

where we indicated that this can only be a function of x — y thanks to the translational 
invariance. Using this and choosing y = —x the Ward identity (|13U|) for the correlation function 
f!132|) becomes 

2 [x^{h - £ 2 ) - (£ (1) - Y,^y v x v ] (0|TO 1 (x)O 2 (-a:)|0) = . (133) 

This is a vector equation. Multiplying it with x M we conclude that necessarily l\ — l 2 . On 
the other hand, choosing the coordinates such that x± = and taking the "minus" projection 
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x M — > x_ we are led to £1 + si = £ 2 + s 2 , or jx = j 2 , as promised. In particular, for the 
two-particle conformal operators defined in Eq. (J48|) we have 

(0|T^ a (o;)^ 3 (0)|0> ~ 8 nm (134) 

since the conformal spins of two operators are different for n 7^ m. This equation looks like 
(and is indeed) the orthogonality relation, but its actual evaluation is a nontrivial task in an 
interacting conformal field theory. In the free theory, however, we can evaluate the l.h.s. rather 
easily, writing the conformal operators in terms of their characteristic polynomials (J37)) and 
replacing the four-point function that arises by the product of two propagators ()127j) . Going 
over to the light-cone limit x 2 = 2x + X-, x± = 0, x + — > with X- fixed, one obtains from Eq. 



Pn(0 ai , d a2 ) W m {d hi dfr) (ai - /3i)- 2jl (« 2 - PiY 2n a =a =0 ~ Km ■ (135) 



01 =02=° 



This can be simplified using the representation 



foo 2j'fc— 1 

(a k - f3 k y 2 ^ = I ds k ^--^ex P (-s k (a k - f3 k )) (136) 







and changing the integration variables to S\ = su\ and s 2 = su 2 with 0<s<oo, < u k < 1 
and u\ + U2 = 1. The integration over s gets factored out, and as the result one finds the 
familiar orthogonality relation ()49)1 for the characteristic polynomials of conformal operators 

(P„|P m ) = / [du] ul^ul^Fnim, u 2 )P m (ui, u 2 ) ~ 5 mn , (137) 



where [du] = duidu 2 S(ui + u 2 — 1). We see now that its existence is a direct consequence of 
conformal symmetry and in particular of the Ward identity for the special conformal transfor- 
mation. Notice that this orthogonality property alone allows one to identify the two-particle 
characteristic polynomials as the Jacobi orthogonal polynomials, Eq. ()47|) . so that their explicit 
calculation is not necessary. 

We recall that the derivation of Eq. (J47|) was based on the transformation (J43|) which led to 
particular simple form of the polynomials F n (K,i, k 2 ) in the adjoint representation, Eqs. ([45)1- 
(14 7j) . Using the above construction we can elucidate the meaning of these polynomials. Consider 
the three-point correlation function x^^x^ 2 {0\Qj 1 (xi)&j 2 (x 2 )®-£ ,: > 2 (0)\0) of the conformal oper- 
ator O^' 32 and two fundamental fields. Taking the light-cone limit in both coordinates x± and 
x 2 and denoting (xi)_ = 1/ki, (x 2 )- = 1/^2, one finds that the dependence of the correlation 
function on K\, k 2 is factorized into 



ai=c«2=0 



P„(ki,« 2 ) = P„(«9 Ql ,5 Q2 )(l-/€ 1 a 1 )- 2jl (l-«:2a2)" 2j2 

J ° kit T ^ 
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It is easy to see that this transformation maps the characteristic polynomial P n into another 
homogeneous polynomial P n (/€ 1 ,K 2 ) of the same degree and with the expansion coefficients 
related to each other according to Eq. ([43)1 . Eq. ()138|) defines, therefore, the characteristic 
polynomial in the adjoint representation. 

3 Conformal Partial Wave Expansion 

3.1 Meson Distribution Amplitudes of Leading Twist 

Hadron distribution amplitudes have been introduced for the QCD description of hard exclusive 
processes jHHl El 13 EH] , e.g., form factors with large momentum transfer, where they play a role 
similar to that of parton distributions in inclusive reactions. Much of the discussion in the 
literature has been concentrated on the pion distribution amplitude of leading twist which is 
the simplest case. In fact, this is the only distribution for which sufficient direct experimental 
evidence is available. It is, therefore, natural to consider this example. 

In an attempt to calculate hard exclusive reactions including a large momentum transfer to 
a pion, one is led naturally to the matrix element of a bilocal quark- ant iquark operator between 
the vacuum and the pion state, of the type 

(0|J(0)[0,oon]7 + 7 5 [oon + x,x]M(x)|vr + (p)) = if nP+ [ due~ iu{p - x) f{u, lnx 2 ) + 0(x 2 ) , (139) 

Jo 

where f w ~ 132 MeV is the pion decay constant. We will neglect the mass of the pion and 
choose the frame of reference such that p^ = p+n^ and = x^n fl + Xj_, x + = 0, so that 
x 2 = —x\, cf. Eq. ()1U|). In the light-cone gauge A + = the Wilson lines disappear and the 
matrix element in Eq. ()139j) can be interpreted as the probability amplitude to find the pion in 
the valence state consisting of the w-quark and (i-antiquark carrying the momentum fractions 
u and u = 1 — u, respectively, and have the transverse separation \x±\. It has been proven 
[HI Ej that this amplitude enters the QCD calculations of the pion electromagnetic form factor 
for large momentum transfers Q 2 and relevant transverse separations are small, of the order of 
x\ ~ 1/Q 2 . 

In order to study this limit, we set x± = and consider a new object, the pion distribution 
amplitude, defined by the bilocal matrix element similar to Eq. ()139j) but taken at exactly 
light-like separations 

(0|rf(0)[0,a] 7+ 75M(a)|7r + (p)) = if vP+ [ du e~ luap +Mu, , (140) 

where we have replaced x_ — > a in order to recover the notations used in Sect. 2. As always 
in a field theory, extraction of the asymptotic behavior creates divergences and the nonlocal 
operator appearing on the l.h.s. of Eq. (|14(Jj) has to be renormalized. In turn, the distribution 
amplitude ^(u, n) is renormalization scale- and scheme-dependent. The original amplitude 
f(u,\iax 2 ) is given by the convolution of (f> n (u, //) with a certain coefficient function c(w,x 2 /i 2 ) 
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which can be calculated order by order in perturbation theory and contains no large logarithms 
provided the scale /i 2 is chosen to be of the order of the transverse separation fi 2 ~ 1/x 2 . In this 
way the small transverse distance behavior of the valence component of the pion wave function 
can be traded for the scale dependence of the distribution amplitude, which is the usual trick of 
the renormalization group. Calculation of this scale dependence is the goal of our analysis and 
we will see that conformal symmetry provides one with a very elegant solution to the leading 
logarithmic accuracy, with almost no effort. 

The renormalization group equation for <p n (u,fi) is known as the Efremov-Radyushkin- 
Brodsky-Lepage (ER-BL) evolution equation jEJ [7] 

d 



ft -jjpfaiu, fi) = J dvVfaviasQififafafi). (141) 



The integral kernel V(u, v) is to leading order in a s given by 

a 



V (u,v) = C F — 

Z1T 



" 1 + — *— J 9(u -v) + -[l+ — ) 9(v - u) 



(142) 



1 — V \ u — V J V \ v — u . 

where the "+" subtraction is defined as 

[V(u,v)]+ = V(u,v)-5(u-v) [ dtV(t,v). (143) 

Jo 

Solving Eq. (j!41j) . i.e., finding the eigenfunctions of the kernel ()142j) . does not look simple at 
first sight. 

Another way to address the same problem is to expand the both sides of Eq. f|14U|) in 
powers of a, so that moments of the distribution amplitude are expressed as matrix elements 
of renormalized local operators: 

(0|<0) 7+ 75 (iD+) n u^)W + (p)) = iU{p + ) n+1 I du (2u - l^K/i) . (144) 



Apart from the flavour content and the extra 75 these are the same leading twist operators 
that enter the OPE for the unpolarized deep inelastic scattering. The difference to the latter 
case is, however, that in the present situation we are not restricted to forward matrix elements 
and hence cannot neglect the mixing of these operators with the operators containing total 
derivatives 

CV M = (td + ) k d(0) 1+l5 (tD + ) n ~ k u(0) . (145) 

Note that the operators with less total derivatives can only mix with the operators with more 
total derivatives, but not the other way around. This implies that the mixing matrix is tri- 
angular and, therefore, the entries on the diagonal are true anomalous dimensions. Moreover, 
they are the same anomalous dimensions that appear in deep inelastic scattering (an extra 75 
has no effect in the flavor non-singlet sector) 

^ ) =Cf{i--—^-—^ + ^-\ , (146) 
V (n+l)(n + 2) t^2 m ) 
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so that 



(°) I a 

fa (a)\ 7n /0 ° 

(P\O r M\P) = (P\O n M\P) (^7^ J , A> = (ll/3)tf c -(2/3)ty. (147) 

Knowing the eigenvalues is helpful but we also need to find the eigenvectors of the mixing matrix, 
which correspond to the multiplicatively renormalizable operators. Finding such eigenvectors 
is of course equivalent to finding the eigenfunctions of the ER-BL kernel. 

The problem is solved by the observation that conformal operators in QCD with different 
values of the conformal spin cannot mix under renormalization to leading order. The reason 
why conformal symmetry "works" in this case is that leading-order operator renormalization is 
driven by counterterms that are tree-level to this accuracy, and therefore retain all the symmetry 
properties of the classical theory. As explained in Sect. 2 the Callan-Symanzik equation that 
governs the dependence on the renormalization scale is nothing else but the Ward identity for 
the dilatation operator. 

In our case the relevant conformal operators are essentially Q 1 ' 1 (x), given in Eq. (|64p. 
correcting for the flavor and Dirac structure: 3 



Qi 1 ^) = (id + ) n [d(x) 1+l5 C 3 J 2 (5+ /d+) u(x)] . (148) 

Because of the flavor structure there can be no mixing with operators built of two gluon fields, 
cf. Eq. (jfi5j) . and all operators built of three and more fields have higher twist so that there 
will be no mixing with them as well. We conclude that the operators in Eq. (|148jl are the only 
relevant ones, and, therefore, they must be multiplicatively renormalized. The beauty of the 
argument is that to arrive to this conclusion we do not need to know the explicit form of the 
mixing matrix, alias the ER-BL kernel. Comparing Eqs. ()144|) and ()148|) . we conclude that the 
Gegenbauer moments of the pion distribution amplitude are given by reduced matrix elements 
of conformal operators 

l duC^ 2 {2u - l)^(u,/i) = ((Qn 1 )), (0|Q^(0)|7r + (p)) = iUp^dQl; 1 )) (149) 



and they must be renormalized multiplicatively, with the anomalous dimensions given in Eq. 
ffTO. 

Last but not least, note that the Gegenbauer polynomials appearing on the l.h.s. of Eq. 
f!148|) are ortho normal with the weight function u(l — u) and form a complete set of functions 
on the interval < u < 1 [3*5] : 

f\uu{l - u) CT{2u - 1) Cl{\2u - 1) = 5 mn (n + 2 1 j ( ^ 3 + ) 2) . (150) 



3 We hope that using the same notation for conformal operators built of quark fields with different flavor and 
Dirac structure will not cause confusion; the choice is always clear from the context. 
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As a consequence, the distribution amplitude 4>^(u, //) can be expanded in a series over Gegen- 
bauer polynomials with multiplicatively renormalizable coefficients 

oo 

<P n (u,fi) = Qu{l-u)Y,K^)C z J 2 {2u-l), 

n=0 

(°) I a 

2(2r> +X) fa ( u) \ 7n /Po 

Mti = ~, \ U J j^r ((Qn'V))), ^) = ^o) Ml • ( 151 ) 
3(n + l)(n + 2) \a s (ii ) J 

This is the final result [SJ Ej- The normalization in Eq. (jl40j) is chosen such that J du(f> n (u) = 1 
which implies <po = 1. The other coefficients are genuine nonperturbative parameters (at a 
certain reference scale) that encode the structure of the pion at small transverse separations 
between the constituents. 

We note in passing that the orthogonality of Gegenbauer polynomials allows one to write 
the ER-BL kernel in terms of the anomalous dimension: 

V (u, v) = - u) ± {n f?x n 3 l 2) if Cy&u - !)£$/> - 1) . (152) 

This representation is valid for Vq acting on arbitrary polynomials on the interval < u < 1. In 
Sect. 4 we will develop a powerful approach to the solution of three-particle evolution equations 
which is based on this kind of representations for the evolution kernels. 
At this place, several comments are in order. 

• To put it simple, the conformal symmetry allows one to resolve the mixing with operators 
involving total derivatives. 

• Whether this is enough to achieve multiplicative renormalization, depends on the problem. 
We have profited from working with flavor non-singlet operators. In the flavor singlet 
sector (e.g., for the distribution amplitude of //-meson) one has to take into account two- 
gluon operators. This will result in a 2 x 2 mixing matrix for each conformal spin. For 
the higher-twist conformal operators given in Eq. ()64|) one also has to take into account 
the mixing with quark-antiquark gluon operators of the same twist and conformal spin so 
that the situation gets complicated. We will return to this problem in the next sections. 

• In order to make use of the completeness condition in Eq. (|15()jl it is important that 
the region of integration over the momentum fraction u in Eq. fll4()j) is precisely < 
u < 1. This is always the case for distribution amplitudes, but can be different in 
other applications, e.g., for generalized parton distributions |HJ HE] • Whenever the 
integration regions do not match, it is not possible to invert the moments and write the 
physical amplitude as a sum over multiplicatively renormalizable contributions, although 
the conformal operators are the same. 

• The expansion in Eq. (jl51j) is formal, in the sense that its convergence properties cannot be 
derived from the symmetry considerations. In QCD, however, the anomalous dimensions 
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of conformal operators are ordered and rise logarithmically with the conformal spin. This 
implies that contributions of operators with high spin are suppressed at large scales and 
asymptotically <f> n (u,fi) — > 6w(l — u) at \i — > oo. This expression is called asymptotic 
pion distribution amplitude. Whether the pion distribution amplitude at hadronic scales 
is close to its asymptotic form, has to be decided by experiment. There exist strong 
indications [UU] that contributions of higher-spin operators are indeed small compared to 
the asymptotic term so that the conformal expansion is well-behaved at already low scales. 
It is worth while to note that convergence of the conformal expansion of distribution 
amplitudes is used implicitly in the proofs of QCD factorization in exclusive reactions. 

• The validity of Eq. (|151|) is not spoiled by quark masses, although the mass terms cer- 
tainly break the conformal symmetry already at the classical level. The reason for this 
is that although finite parts of the matrix elements of the operators (jl48j) and therefore 
the coefficients </> n (/x) depend on the quark masses (e.g., they are different for n and K 
mesons), the infinite, UV divergent contributions to the matrix elements are removed by 
the mass-independent counterterms, so that the renormalization group behavior is not 
affected. 

The expansion (jl51|) is generic and is valid for all leading twist quark-antiquark meson distri- 
bution amplitudes, albeit with different coefficients and anomalous dimensions. As a further 
example, consider distribution amplitudes of a vector particle, e.g., the p-meson (SUES- There 
exist two independent distribution amplitudes corresponding to p-mesons with longitudinal 
A = and transverse A = ±1 polarizations: 

(0|J(0) 7+ u(a)|p + (P,A)) = f p m p ef f due~ map + ^.{u^) , 



o 



(0|J(0) ( x ±+ n(a)|p + (P,A)) = iflefp+ [ due^<j> T Au,v) , (153) 



where is the polarization vector and f p and fj are the corresponding decay constants. 

Since 7 j_ commutes with the quark spin projection operators (joTj) . both 7+ and a± + project 
onto twist-two operators built of two "plus" quark fields, so that the conformal operators only 
differ by this substitution. Repeating the above derivation, we end up the same expansions 



Mu^) = 6«(i-«)E^W(-r^) cf 2 (2 M -i), 

0>,/i) = euCl-uJ^^W^^^J C 3 J\2u-l), (154) 



n=0 

T(0) 



where in the second case 7 is subtracted in order to compensate for the scale dependence of 
fj. The anomalous dimensions 7^ are the same as for the pion while are equal to 

7?°» = G F (l + 4f iV (155) 



m=2 
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From all these examples we see that the task accomplished by the conformal symmetry is 
the separation of variables. A very close analogy is provided by nonrelativistic quantum me- 
chanics in a spherically symmetric potential. As well known, the 0(3) symmetry allows for the 
separation of angular and radial degrees of freedom by going over to the basis of states with 
given orbital angular momentum I. For fixed I one is left with a one-dimensional Schrodinger 
equation for the radial coordinate, whereas the angular dependence is included in spherical 
harmonics which are the irreducible representations of the group of rotations. In our case, the 
conformal symmetry invites to go over to the basis of states with given conformal spin, confor- 
mal partial waves. For each spin the dependence of the wave function on transverse coordinates 
is given by a simple renormalization group equation and the dependence on longitudinal coor- 
dinates (momentum fractions) is included in orthogonal polynomials which form an irreducible 
representation of the collinear conformal group, SX(2,R). 

The group-theoretical content of this problem becomes especially transparent using the 
basis of one-particle conformal states defined in Eq. (|56p. A generic multi-particle distribution 
amplitude can be defined as the matrix element of a multi-local operator 

rco pea 

(0|$i(«i) . . .$ m (a m )\p) = / dpi... dp m S(S2pi-p) exp[-i(aipi + . . . + a n p m )} 

Jo Jo 

oo 

x (°\ ( Pk 1 --- ( Pk m \p) (Pi\juh) ■■■{Pm\jm,k m }, (156) 

k \ , .... Aj ffi — 

where we used the momentum conservation, the fact that all states have positive energy, hence 
Pk = {p+)k > 0, and inserted the expansion (jSHJ). The matrix elements (0|</9fc 1 . . . ipk m \p) are just 
constants; call them Wfci,...fc m - The basis vectors (pi\jiki) (JoT)|) are homogeneous functions; hence 
one can redefine them in terms of momentum fractions Ui = Pi/p and the whole expression 
becomes, up to the overall normalization 

(0|$i(ai)...$ m (a m )|p) 
<f)(ui, . . . ,u m ) 

where J[du] = Ldu\ . . .^du m Uj — 1). What we would like to see, on the other hand, is the 
expansion of the distribution amplitude in terms of the irreducible representations of SX(2,R) 
with given conformal spin j 

oo 

<f>( u i)= J2 J2"'ti,R)(Ui\3,0)R (158) 

j=jl+—+3m R 

where the summation over R goes over all (degenerate) representations with the same spin. 
Going over from Eq. ()157j) to Eq. ()158|) is the classical problem of spin summation, with 
the only difference that the spins "live" on the hyperboloid since the algebra of the collinear 



[du] j p{aiUl+ "- +anUm) (j)(u h . . .u m ) , 

oo 

2j ^fci,.-fem ( u i\ji, h) ... (u m \j m , k m ). (157) 



fcl . . . ..km — 
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conformal group SL(2, R) (|18|) coincides with the algebra of the Minkowski group in 2+1 
dimensions, 0(2, 1). Restricting oneself to infinitesimal group transformations one can say that 
the conformal spin j corresponds to the "mass" and |jo| to the "energy" of a particle in this 
internal space. Using this analogy it becomes intuitively obvious that the lowest conformal 
spin ("mass") for a multi-particle state is just the sum of spins ("masses") of the constituents; 
this state is nondegenerate and is given by the product of one particle states "at rest" , that is 
with the lowest values of the "energies" \j \ = j, alias k = 0, in our notation. This lowest state 
defines what is called the asymptotic distribution amplitude in the general case [5] 

Mm) = ( Ul \j — ji + ■ ■ ■ + j m , o) = T^^^ ^ 1 " 1 ^ 2 " 1 • • • ^T' 1 > ( 159 ) 

where the normalization is chosen such that j[dui] 4>g, s (ui) = 1. The constructions of states 
with higher spin reduces to finding the corresponding Clebsch-Gordan coefficients 

\j, k) = ^2 Cjit 1 ,...j m fe m ( n lbl^l) • • • ( U m\jmk m ) (160) 
h\ -)-... -)- k m = k 

In particular, for two particles 

(u lU2 \j,k) = {u l + u 2 ) k ("l) fcl ( ^^O^lljl,^)^^,^) 

ki+k 2 =j+k~j 1 -j 2 

~ (m + u 2 ) "Muu u. 2) if) 'f 1 (^l) , (161) 

J J J \Ui +u 2 J 

cf. Eq. (|47jl . Note that the "raising" operator L + = L± t + + . . . + L mt+ acts as a multiplication 
operator on the distribution amplitudes, L + (f>(ui) ~ (u\ + . . . + u m )(f)(ui), and the prefactor 
is equal to unity because of the phase-space condition J^Mj = 1. Hence only the lowest 
"energy" state with k = is physically relevant. In the general situation with more than 
two particles it can be shown that higher spin states are obtained by multiplication of the 
asymptotic distribution amplitude by a polynomial of degree n = j —j\ —j 2 and the polynomials 
corresponding to different conformal spins are mutually orthogonal with the weight function 
4>as(ui) ()159|) . A convenient complete set of orthogonal polynomials for three particles is defined 
in Appendix A. Note that the whole construction is equivalent to finding of the corresponding 
local conformal operators. 

3.2 Higher Twists 

Contributions to physical cross sections that are suppressed by a power of the large momentum 
are generically referred to as higher twists. The most important question to be addressed in 
connection with higher twist contributions is always to prove that QCD factorization can be 
extended to include power-suppressed corrections. If this is the case, the second question is to 
obtain a general classification of the relevant hadronic quantities in terms of the independent 
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nonperturbative parameters. Since the number of parameters usually prolifiterates compared 
to the leading twist, a practical task is to establish their hierarchy and build selfconsistent 
approximations with minimum nonperturbative input. 

Conformal partial wave expansion proves to be indispensable for the construction of higher- 
twist distribution amplitudes. At present, this program is completed for pseudoscalar [""J Ell 
and vector mesons jHHJEHL for the real photon [HBl E7J and partially also for the nucleon [68 . 
The summary of results goes beyond the tasks of this review. Instead, we will try to explain 
basic ideas and techniques beyond such calculations, and also indicate existing problems. Our 
first example will be the construction of selfconsistent approximations to pion distribution 
amplitudes of twist-three. This is the simplest case which is, on the other hand, sufficiently 
general; applications to twist-four and to other hadrons follow the same scheme. Second, we 
discuss a recent work fo"*"] 170] concerning convergence properties of the conformal expansion for 
higher twist distributions using the concept of renormalons |71| 172*] . Finally, we consider meson 
mass corrections and in particular the so-called Wandzura-Wilczek contributions to distribution 
amplitudes, and how they can be included in the conformal expansion. To this end twist-three 
distributions of the p meson are considered. 

3.2.1 Pion Distribution Amplitudes: Twist-three 

Pion distribution amplitudes of twist-three are defined as matrix elements of nonlocal twist- 
three operators. One finds two two-particle distributions 

f m 2 r 1 

(0\d(0)i l5 u(a)\7c + (p)) = J * * / due- lua ^<P P (u,fi), 

m u + m d J 



— ICY f 771 ^ i 

(0|d(0)<v7 5 u(a)k + (p)> = a( J \ * A p^-p^) / due-^+^u,//), (162) 



6(m u + m d ) Jo 
and one three-particle distribution 

(0\d(a 2 )a IMU 'j 5 gG a p(a 3 )u(ai)\7r + (p)) = i p a (Pf,g uf 3 - p v 9^p) ~ (« «-> P) 



x [d Ui \e- ip+{aiUl+a2U2+a3U3) (f) 37T (u z ,fx) . (163) 



The normalization is chosen such that fdu<j) p (u) = jdu(f) a ( u ) = !• We will use a shorthand 
notation 

R = /X/K + m d) - -(uu) - (dd) , (164) 

where (uu) ~ (dd) is the quark condensate. We will tacitly imply the chiral limit and neglect 
quark mass corrections in what follows. 

The three functions <fi p (u), <f> a (u) and (p^iui) are related to each other by equations of 
motion "73*] [5] . One way to see this is to observe that relevant matrix elements of chiral-odd 
operators (with even number of 7-matrices) correspond to one "plus" component of the quark 
field and one "minus" component. This also explains the twist assignment, since going over 
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from a "plus" to a "minus" component adds one unit of twist. As well known in context 
of the light-cone quantization [73], "minus" components of the spinor fields do not present 
independent degrees of freedom, but can be eliminated using the Dirac equation. A convenient 
way to derive these relations is to make use of the operator identities 



d 



d p \ d(-x)a fiv j 5 u(x) }> = i-^;\d(-x)>y 5 u(x)^ + J vdv d(-x)x p gG pu (vx)j 5 u(x) 



i 



_d_ 

dx l 



+ i I dv d(— x)x p gG pp (vx)a fll/ / y 5 u(x) 



d(—x)a flv j 5 u(x)^ = id„^d(-x)>y 5 u(x)^ + J dv d{-x)x p gG pv (v x)^u(x 



+ i I v dv d(— x)x p gG PPj (yx)a pv ^u{x) (165) 



i 



for the total translation and dilatation of the relevant nonlocal operators; d v stands for the 
derivative with respect to the total translation, for example 

d - 

d v [d(-x)^y 5 u(x)] = — [d(-x + y)^u{x + y)} . (166) 

Taking in Eq. ()165j) the light-cone limit x 2 — > and the matrix elements between vacuum and 
the one-pion state, one obtains the system of two differential equations connecting the three 
functions (fi p , (p a and (p 3n , see jHUHll for details. 

Our aim is to build a selfconsistent description of the distribution amplitudes in terms of a 
minimal number of nonperturbative parameters. To be precise, by a "selfconsistent" description 
we mean the following: 

• The constraints that follow from equations of motion have to be fulfilled identically. That 
is, the relations between p , <p a and (p^ that follow from the exact QCD operator identities 
(|165|) have to be satisfied. 

• The structure of the distribution amplitudes has to be preserved by the evolution. For 
example, if we choose a model for the distribution amplitude at a certain scale as a 
polynomial, the degree of this polynomial remains the same at all scales. We also require 
that the renormalization group equations for the relevant parameters take the simplest 
possible form. 

Both tasks are achieved if the distribution amplitudes are expanded in contributions of increas- 
ing conformal spin. In the previous section we have seen already that conformal expansion 
simplifies the evolution equations: To one-loop accuracy, the renormalization of relevant oper- 
ators is driven by counterterms that are tree-level and retain the symmetries of the classical 
Lagrangian; it follows that operators with different conformal spin cannot mix with each other. 
For the present case there is a complication: The three-particle representations of the collinear 
conformal group are degenerate, there exist multiple operators of the same conformal spin and 
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the symmetry does not prevent their mixing. The mixing matrix is, therefore, not completely 
diagonalized but rather separated in smaller blocks. More importantly, QCD equations of mo- 
tion do not receive any quantum corrections whatsoever and, therefore, respect all symmetries 
of the classical theory. As a consequence, if two operators are related by the equations of 
motion, their matrix elements must have the same transformation properties under conformal 
transformations and, in particular, the same conformal spin. This implies that the relations 
between <p p , <p a and 4> 3n that follow from the operator identities in Eqs. (|165j) can be solved 
order by order in the conformal expansion. 

After these general remarks, let us proceed with the explicit construction. We start with 
the three-particle distribution ()163|) . It is easy to see that the twist-3 part of interest can be 
separated by the Lorentz projection da + ±'y^gG + ±u and corresponds to the maximum possible 
conformal spins of the three constituent fields: jj = 1, j u = 1, and j g = 3/2. The nonlocal 
composite operator in Eq. (j!63|) can, therefore, be expanded in contributions with conformal 
spin j = jd+ju+jg + n = 7/2 + n, n = 0, 1, 2, . . .. The contribution with the lowest spin j = 7/2 
can be referred to as the asymptotic distribution amplitude and higher spin contributions are 
given by the corresponding orthogonal polynomials. Using the master-formula in Eq. (j!59|) and 
taking into account that the distribution <p 3n («,-) has to be symmetric under the interchange of 
ui and u 2 because of G-parity, we obtain 

<M«i) = 360^2^3 c/ /2 + u 9/2 ^(7u 3 - 3) + uof /2 (2 - A Ul u 2 - 8u 3 + 8u 2 3 ) 

+ w" /2 (3uiU2 - 2w 3 + 3u 2 3 ) + ... . (167) 

Here io\ are scale-dependent nonperturbative coefficients; the superscript refers to the conformal 
spin. Conformal symmetry implies that to 7 ^ 2 and to 9 ^ 2 are renormalized multiplicatively (to 
leading order), while u\^ 2 and lo^ 1 ^ 2 can mix with each other. These properties are indeed 
confirmed by the explicit calculation, the anomalous dimensions and the mixing matrix for 
lo 1 ^ 2 and uj^ 1 ^ 2 can be found in [Oj. The expansion can easily be extended to include higher 
conformal spins using the conformal basis defined in Appendix A. 

Going over to the two-particle distribution amplitudes ()162j) we have, first of all, to separate 
contributions with different (anti)quark spin projections. To this end we define two auxiliary 
distribution amplitudes |2] 

(0|J(0)* 7 -7 + 75«(a)|7r + (p)> = R I due- luap +<P n (u, fj) , 

Jo 

(0|J(0)z 7+ 7-75M(a)|7r+(p)> = R [ due~ lua ^ cj> n (u, fi) . (168) 

Jo 

It is easy to see that, e.g., 4>^{u) corresponds to the spin projections sj = —1/2, s u = +1/2 so 
that jj = 1/2 and j u = 1. It follows that the conformal expansion of (t)^{u) and 4>n(u) goes 
over Jacobi polynomials (}48|) 

R M u) = 2u\nf + .f P r\0 + ^^tt) + -T^ 1] (0 + -T^it) + ■■ 
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R<P n (u)=2u 



K 



3/2 

u 



.5/2p(l>»), 



7/2 p (l,0) 



(0 + « 



9/2 p (l,0) 
U ^3 



(0 



11/2 p(l,0) 

U 



p 4 ^(0 + ... ;;69) 



with nonperturbative coefficients k 3 . Here and below in this section we use a shorthand notation 



;i?o) 



The superscripts j = 3/2, 5/2, . . . refer to the conformal spin. On the other hand, obviously 



M u ) = n[<f>n( u ) + M u )l 



d 

du 



(j) a {u) = 6[0 u (m) - (f>n(u)) 



(171) 



so that the knowledge of 4>^{u) and 4>n( u ) is sufficient to determine both distributions <p p (u) 
and (p a {u) A . 

Finally, we can use equations of motion or, equivalently, the operator identities ()165j) in 
order to relate the expansion coefficients k j and After simple algebra (see jUJ) one obtains 
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tt 



- K 9 { 2 = 30o// 2 



— K 



5/2 
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11/2 _ 

IT = 







K 



7/2 

IT 



_7/2 



30 o; 7 / 2 , 



11/2 /o/o\r/1 n /2 H/2i 

/c n ' = (3/2) [4^ - uj 2 \ 



(172) 



etc. We see that everything fits together, equations of motion only relate contributions with the 
same conformal spin, as expected. In particular, the conformal spin of the quark condensate 
contribution R = —2(^>^)) = —{ip^-^ + ip) — {ip^+^-ip) is indeed j = 1/2 + 1 = 3/2. Note that 
k 11 / 2 is not multiplicatively renormalized, but rather is given by a certain combination of the 
two existing operators with spin 11/2. The reason for this is that the quark- ant iquark twist-3 
operators entering ()168|) mix with the quark-antiquark-gluon operator (j!63J) of the same twist. 

The conventional distribution amplitudes defined in Eq. ()162j) can easily be restored using 
Eq. (IT7TD . One obtains 



R(j) p {u) 
R(j)<j(u) 



P + 30c 7 / 2 C 2 1/2 (O + |[4^ 11/2 



6u(l-u)<R + 



50,7/2 _ 1^9/2 

2 



u?' 2 - 2c 9 / 2 



c 4 1/2 (0 



11/2 



10. 



11/2 



C 4 3/2 (£) U173) 



1 /2 3 /2 

where C n (£) and C n (q) are Gegenbauer polynomials. Note that also in this case the ex- 
pansion is organized in a natural way in terms of orthogonal polynomials. This is no more a 
conformal expansion, however. The coefficient in front of each polynomial is equal to the sum 
of two contributions with neighboring conformal spins. In particular, the coefficient in front of 
C^ 2 in the second equation in (|173|) will be modified once contributions of higher spin 13/2 are 
taken into account. For this reason, and also to avoid superficial cancellations between 4> p (u) 
and 4>a{u), the distribution amplitudes 4>^{u) and 4>n( u ) are preferable for the description of 
twist-three effects in exclusive reactions with pions. 



4 The integration constant in the second of equations in Eq. (|171f) is fixed by the normalization condition 
J du4> a {u) = 1, so there is no ambiguity. 
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3.2.2 Renormalon Model of Higher- Twist Distribution Amplitudes 



The next problem to address is the convergence property of the conformal expansion for higher- 
twist hadron distribution amplitudes. Assuming that the anomalous dimensions of higher- 
twist operators increase logarithmically with the conformal spin, similar to the leading twist, 
and repeating the argumentation of Sec. 3.1, one is led to the conclusion that the conformal 
expansion must converge starting at a certain high renormalization scale. From the practical 
point of view this argument is not very assuring, however, since the suppression of higher spin 
operators by the QCD anomalous dimensions is rather weak, at least in perturbation theory. 
For leading twist, the better argument for the convergence of the conformal expansion at already 
the scales of order several GeV comes from experiment: The CLEO data on the 77*71 transition 
form factor [fJD] strongly favor the pion distribution amplitude that is not very different from 
the asymptotic form. For higher twist, a similar direct experimental verification is impossible. 
However, if convergence of the conformal expansion for leading twist distributions (at a given 
scale) is taken for granted, then one can derive a reasonable upper bound for the contributions 
of higher conformal spin operators of twist 4 using the technique of renormalons [7T1 I72| . 

The basic idea is that twist-four contributions to physical amplitudes are in one-to-one cor- 
respondence to the contributions of low momenta in perturbative Feynman diagrams. Consider 
the following example: 

(0|d(0) ^ l5 d(x)\7f + (p)) (174) 



if 7r (px) / due 



-iupx 







1 + ^ Cfc ( M > V 2x2 ) a s (A* 2 ) ® <&r(tt, « 2 ) + x 2 9i(u, « 2 ) + 0(x 4 ) 



k=l 



where x 2 Aq CD <C 1, <8> stands for the convolution, and the gauge factor is assumed along the 
straight line connecting the quark and the antiquark. The expansion on the r.h.s. of Eq. (|174|) 
is nothing else but the operator product expansion and the contributions with the increasing 
powers of x 2 correspond to the increasing twist. The relevant matrix element of the leading 
twist operator defines the usual pion distribution amplitude 4> n (u) and it is multiplied by the 
coefficient function that can be calculated order by order in perturbation theory. gi(u, /J, 2 ) 
is usually referred to, somewhat inaccurately, as the two-particle pion distribution amplitude 
of twist-four. Its physical interpretation is to describe the k\ distribution of the quark and 
the antiquark in the pion, and it can be expressed in terms of three-particle quark-antiquark- 
gluon twist-four distribution amplitudes using equations of motion, similar as for the twist-three 
distributions considered above. Using these relations, and expanding the relevant three-particle 
distributions in contributions of conformal operators one can work out the conformal expansion 
for gi(u), see jH] for the derivation. Taking into account the contributions with the lowest and 
next-to-lowest conformal spin, the result reads 5 



9i{u,n) = g j=3 (n)u 2 u 2 (175) 



u 1 1( 2+ ilrin) - 10m 3 Inu ( 2-3u + ^u 2 J + 10m 3 Inu (/l-Zu+^u 2 



5 In notation of Ref. H g^ 3 = 5/6S 2 , g^>= 4 = l/2eS 2 . 
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where g J=3 (/i) and gi =A (fi) are the reduced matrix elements of conformal twist-four operators. 
They are renormalized multiplicatively, with known anomalous dimensions. Note that the 
j — 4 and higher spin contributions are not given by polynomials in the momentum fraction 
in this case, which implies that twist-four two-particle conformal operators are not given by 
polynomials in the covariant derivatives (i.e., they are non-local), cf. [TSU^]. The question that 
we want to address now is whether the approximation in Eq. (|175|) is sufficient, or higher spin 
operators can play an important role. 

Let us assume for a moment that the scale separation in Eq. (j!74|) is done with a momentum 
cutoff, that is contributions of high loop momenta \k\ > fi are included in the coefficient 
functions Ck{u, fi 2 x 2 ) and contributions of low momenta \k\ < fi in the matrix elements. The 
coefficient functions can be expanded in the Laurent series in /i 2 x 2 <C 1 and apart from the 
usual logarithmic terms ~ ln fc x 2 fi 2 will also contain the power corrections, schematically 

oo \ / oo \ 

1 + c k( u > A^KV) ) — 1 + Yl Cfc ( M ' ln (Z^) - d{u)fi 2 x 2 + , (176) 

k=i J \ k=i J 

where the terms ~ fi 2 x 2 are also accompanied by a series in ln fj, 2 x 2 that we do not show for 
brevity. As well known, the logarithmic scale dependence of the coefficients is compensated 
by the scale dependence of matrix elements, in our case the leading twist pion distribution 
amplitude, so that the product of the leading twist coefficient function and the relevant operator 
matrix element does not depend on the factorization scale to logarithmic accuracy . The crucial 
observation is that independence of the result on the factorization scale to power accuracy 
involves cancellations between contributions of different twists. Technically, this cancellation 
proceeds as follows. As mentioned above, gi(w,/i 2 ) is determined by reduced matrix elements 
of twist-four operators that have not only the usual logarithmic, but also the quadratic UV 
divergence. If, as we have assumed, the matrix elements are calculated with an explicit UV 
cutoff, the results for the matrix elements will be, schematically 

gx(u, /iV) -> fi 2 d(u) ® <f>n(u) + Aq CD gi(u, ln/iV) , (177) 

with the same function d(u) as in Eq. (j!76J) . The result is that power-suppressed infrared 
contributions to the coefficient functions cancel against ultraviolet quadratically divergent con- 
tributions to the matrix elements of higher twist operators. If the factorization is done using 
dimensional regularization (and MS-subtraction) instead of the momentum cutoff, then power- 
like scale dependence does not arise, but the perturbative series becomes factorially divergent 
in high orders. Divergences of the perturbative series in high orders are called renormalons and 
they can be due to both the infrared and the ultraviolet regions. In this language, one observes 
that infrared renormalons in the coefficient functions cancel against the ultraviolet renormalons 
in the matrix elements of higher twist operators. 

Without going into details, we note that the function d{u) can be calculated in perturbation 
theory and fi 2 d(u) ®(j) w {u) can be taken as an estimate ("renormalon model") for the twist-four 
distribution amplitude. The leading-order result reads [03 [70] 

' x dv 



g[ en (u) = ,i 2 <u 
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Figure 1: The twist-four two-particle pion distribution amplitude g\(u) in the renormalon model 
f!178|) (solid line) compared with contributions of the first two orders in the conformal partial 
wave expansion (j!75|) (dashed line) with the coefficients estimated using QCD sum rules p. 
The asymptotic distribution amplitude corresponding to the first term in Eq. (J 175)) is shown 
by the dotted curve. The normalization is adjusted so that f Q dugiiu) = 1. 



+ «r$[ 1 + (s- 1 ) to ( 1 -j)]«4 <™ 



where the scale /i 2 has to be of the order of Aq CD but otherwise is arbitrary. This expression 
corresponds to taking into account contributions of all conformal operators with arbitrary spin, 
assuming "ultraviolet dominance" [721 EO] of the corresponding matrix elements and neglect- 
ing the (logarithmic) effects of the anomalous dimensions altogether. Because the anomalous 
dimensions suppress contributions with higher conformal spins, the renormalon model is likely 
to overestimate them, and thus can be viewed as an upper bound for higher spin contributions. 

The comparison of the renormalon model (|178|) with the contributions of the lowest two 
conformal partial waves ()175|) is shown in Fig.Q where we used the asymptotic pion distribution 
amplitude 4> n (u) = Quit for the evaluation of Eq. (|178|) and adjusted the overall normalization. 
We see that the global shape of gi(u) in both models is very similar, and the only essential 
difference concerns the behavior at the end points: In the renormalon model gi(u) ~ it(l — u) 
for u — > 0, 1, whereas each term in the conformal expansion has the ~ u 2 (l — u) 2 behavior. The 
difference means that the conformal expansion is not converging uniformly at the end points, 
in the renormalon model approximation. As mentioned above, contributions of higher spin 
operators alias the end-point behavior of the renormalon ansatz will be tamed by taking into 
account the anomalous dimensions of higher twist operators, so that what we see in Fig.^has to 
be considered as an upper bound. Phenomenological consequences of this different behavior still 
have to be investigated, but in any case the renormalon approach provides a viable alternative 
to the uses of conformal expansion truncated at the first few orders. 
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3.2.3 Wandzura-Wilczek Contributions and Meson Mass Corrections 

Last but not least, we want to consider meson mass corrections and discuss how they can be 
included in the conformal expansion. To this end we will consider a subset of the distribution 
amplitudes of the p-meson, defined by the matrix elements of nonlocal operators with even 
number of 7-matrices in between the quark and the antiquark operators. This is a close analogue 
to the example of pion twist-three distributions considered above, and the difference is, first, 
that p meson has spin one, and, second, it has a mass which is non-negligible in comparison to 
the mass scale of quark-gluon correlations. Following [33] we define 

(0|d(0) ( r^(a)|p + (P,A)) 



(0|J(0)u(a)|p + (P,A)> 

where P 2 = m 2 , p^ = P^ — n fl m 2 p /(2p + ), p 2 = 0, and the polarization vector eft* of the p- 
meson is decomposed into projections onto the two light-like vectors n M , n M = p^/p+ and the 
orthogonal plane: 

e& = efn^ + e { %p + e%. (180) 

The function (f)±(u) is the leading twist distribution amplitude of the transversely polarized 
p-meson, hu\u) and hi (u) are twist-three distributions for the longitudinally polarized p- 
meson, and h^{u) is the twist-four distribution, for transverse polarization. In addition, there 
exist one three-particle twist-three distribution for longitudinal polarization, and several three- 
particle (and four-particle) twist-four distributions for transversely polarized p mesons. Explicit 
expressions can be found in [3*3*1 [f"""5] , we will not need them in what follows. 

The distributions nP(u) and hi?\u) are very similar to the tensor (fi a (u) and the pseudoscalar 
4> p {u) twist-three distributions of the pion that we have just constructed. In full analogy with 
the latter, hiP (u) and (u) contain a "genuine" twist-three contribution that can be related 
to the existing three-particle distribution using equations of motion and following the scheme 
described above. Note that these relations are between the two-particle and the three-particle 
distributions corresponding to the same (longitudinal) polarization. A peculiar new feature 
is, however, that in addition to this "genuine" twist-three part, hif*(u) and hi? (u) contain 
one more contribution that can be expressed through the distribution amplitude 4>t{u) with 
a lower (leading) twist-two and different (transverse) polarization. Such contributions have 
first been found in [77] for the structure function g 2 (x,Q 2 ) of the deep-inelastic lepton-hadron 
scattering from transversely polarized nucleons, and are generically referred to as Wandzura- 



if. 



(e ( iy - O m ) / due-^M^^) 



+ 2^( e ±y ~ effin^m; / due- ma ^h 3 (u, p) 



la t T 2 (A) 



due- luap+ h { °\u,p), 



(179) 
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Wilczek contributions. We can write 

hf(u) = hf WW {u) + hf qG \u) , hf{u) = h\; )ww (u) + h^ G \u) , (181) 

and will concentrate on the Wandzura- Wilczek parts; the construction of "genuine" twist-three 
contributions is fully analogous to the pion case. 

In order to derive the Wandzura- Wilczek contributions one starts with the same operator 
identities ()165|) (omitting the 75-matrix), takes the appropriate matrix elements and neglects 
contributions of qqG operators. The relations that arise in this way have a unique solution 

hf ww {u) = (2u - 1) ( r d v _ r dv 



h^ WW (u) = 2(u r dv *M+ u fdv+M). (182) 



These are the analogues of the Wandzura- Wilczek contributions to the nucleon structure func- 
tions gx(x, Q 2 ) and h L (x,Q 2 ) [32j. 

Our goal is to understand the structure of these contributions in the context of the conformal 
expansion. Similar as for the pion, we need to go over to auxiliary amplitudes with fixed spin 
projections: 

_(A) ,1 

<0|J(0)7 + 7-«(a)|p + OP,A)> = f p m 2 P —\ du<r iua ^h n {u) , 

P+ Jo 

e (A) ,1 

(0\d(0)^ 1+ u^)\P + (P^)) = f P ™\—\ due- iua *>+h u (u), (183) 

P+ Jo 

which are related to (u) and hlf' (u) by 

h n (u)= + -_JLAi, fc iT („) = + __JL^. (184) 

The conformal expansion of h^(u) and h^(u) is straightforward and is given by 

00 00 
h n (u) = 2uJ2h r !t 3/2 P^ \0 , M«) = 2uJ2hT /2 Pi 1 '°\0 . (185) 

n=0 n=0 

On the other hand 

00 

± H = 6^^01 +2 ^ 2 (O. (186) 

n=0 

In the both cases the superscripts (j = n + 3/2 and j = n + 2) refer to the conformal spin. 
Substituting these expansions in Eq. ()182|) . after some algebra one obtains [SSI 

rrn+3/2 _ 3(77, + 1) +2 „_,_., K + i , _ n .-, I 

ww ~ 2nT3~ 1 ' ww ~ °- 1 1 L ' l n - u > 2 > 4 '--0, 

o-n+3/2 _ 3(n + 1) +1 ) "-■'- _ » + 2 , _ 1 o - > ,| S7) 

-°ww — 2n + 1 ' ww o„ 1 q 'i ' l n ~~ °' °' • ■ ■) ' 



n+3/2 


3(n + l) +1 


WW ~~ 


2n + 1 <?± 


ra+3/2 


3 ( n + 1) -n+2 


WW ~~ 


2n + 3 ^ 
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where H^^/ 2 and li™ 2 are the Wandzura-Wilczek parts of the G-parity even and odd com- 



l WW aLLU - "WW 

binations of the coefficients: 



jjn+3/2 _ 1 

2 



h n+3/2 _ (_ 1 )n^n+3/2 ^ ^n+3/2 = - ^n+3/2 + (_ 1 )^n+3/2 (lgg) 



1 r 



2 



n+3/2 / r , n , n+3/2 



"U ^ V IT 



We see that the conformal spins on the l.h.s. and the r.h.s. of Eq. (|187|) do not coincide which 
calls for an explanation. 

The basic idea is the following j^Sj. The Wandzura-Wilczek relations (|182j) are in fact 
consequence of the Lorentz symmetry. Working out such relations essentially corresponds to 
reexpressing matrix elements of conformal operators between vacuum and the longitudinally 
polarized p-meson in terms of the similar matrix elements involving the transversely polarized 
meson, or vice versa. In our context it is important that the transversely polarized state 
is obtained (in the p-meson rest frame) from the longitudinally polarized state by the spin 
rotation which does not commute with the generators of the collinear conformal group. This 
rotation gives rise to the shift in conformal spin and exactly explains the mismatch appearing 
in Eq. (j!87|) . The moral is that the conformal symmetry is realized for the Wandzura-Wilczek 
contributions as well, but to see this one has to supplement the conformal classification of 
operators by conformal transformation properties of the meson states. 

It is convenient to work in the helicity basis: p-mesons with A = ±1 correspond to the two 
states with transverse polarization, A = stands for the longitudinal polarization. Introducing 
conformal eigenstates L 2 \j, k) = j(J — l)\j, k), L \j, k) = —(j + k)\j, k), cf. ((57}). we can separate 
different terms in Eq. (J186|) as 



tn+2 



oc (n + 2, 0\p(P, A = ±1)) , (189) 



where the proportionality factor is irrelevant for what follows. In the p-meson rest frame the 
A = +1 state is related to the A = state by the spin rotation 

|p(P = m p , A = +1)) oc (M 23 + iM 31 ) \p{P = m p , A = 0)), (190) 

where M23 + iM^i is the appropriate step-up operator of ordinary angular momentum. A 
fast-moving p-meson \p(P, A = ±0)) is then obtained by the Lorentz boost in the appropriate 
direction P p = (P , 0,0, P 3 ): 

\p(P,\))=U(u)\p(P = m p ,\)) (191) 
with U(u) = e - lujM ° 3 , tanh(cu) = P 3 /P - We can thus write 

01 +2 oc(vI/|p(P,A = O)>, (192) 

where 

=U(Lu)(M 23 + iM 31 )U~ 1 (uj)\n + 2,0) . (193) 
Working out the commutation relations one obtains 

|^) = a\n + 3/2, 0) + c 2 \n + 3/2, 1) + c 3 \n + 5/2, 0), (194) 
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where are certain numerical coefficients. Eq. (|194j) shows that the contributions of 0" +2 
to the distribution amplitude of the longitudinal p-meson correspond to the conformal spins 
j = n + 3/2 and j = n + 5/2, which explains the pattern appearing in Eq. (|187j) . This derivation 
is valid also for chiral-even distribution amplitudes, see [33]. 

The structure of meson mass corrections to twist-four distribution amplitudes is much more 
complicated [US] as it turns out that they have several sources. First, mass corrections are 
generated by the contributions of leading twist operators in which case the conditions of sym- 
metry and zero traces imply a certain Lorentz structure that multiplies the relevant reduced 
matrix elements. Such corrections are the direct analogue of the Nachtmann's correction to 
deep inelastic scattering. This result is not complete, however, because in exclusive processes 
one has to take into account higher-twist operators containing total derivatives, and vacuum- 
to-meson matrix elements of such operators reduce to powers of the meson mass times reduced 
matrix elements of leading twist operators. The operators containing the square of the total 
derivative times the leading twist operator are the simplest ones and can be taken into account 
in a systematic way, see [S3]. The major complication arises because of operators that are given 
by (or can be reduced to) the divergence of the leading twist conformal operators considered in 
Sec. 3.1. As well known [7S], conformal operators are divergenless in the conformal-invariant 
theory, so that in QCD their divergence must be proportional to the strong coupling and in- 
volve the gluon field. Indeed, one obtains that this divergence can be expressed in terms of 
contributions of quark-antiquark-gluon operators [3HJ EUl EH] , schematically 



Sandwiching Eq. (|195|) between vacuum and the meson state, one obtains nontrivial relations 
between reduced matrix elements of quark-antiquark-gluon operators on the r.h.s., and the 
meson mass squared times the leading twist matrix elements on the l.h.s. This implies that the 
meson mass effects and the effects of quark-gluon correlations in mesons cannot be separated 
from one another in a meaningful way: The approximation of taking into account the meson 
mass corrections and neglecting contributions of quark-gluon operators appears to be not self- 
consistent, see jHS] for a more detailed discussion. The similar problem is present for photon 
distribution amplitudes |3H] R)T] in which case one would like to separate contributions of higher- 
twist quark-photon operators (contact terms) from the quark-gluon ones, and also for processes 
involving off-forward (generalized) parton distributions if one tries to separate higher-twist 
corrections that are proportional to either the nucleon mass squared, or the four-momentum 
transfer t = (P2 — Pi) 2 [Z^]- Here Pf and P% are the momenta of the initial and the final state 
nucleon, respectively. 

4 Hamiltonian Approach to QCD Evolution Equations 

4.1 Conformal Symmetry and Hermiticity 

We have seen in Sect. 3.1 that the kernel ()142j) of the ER-BL equation for the scale dependence 
of pion distribution amplitudes can be written in terms of the anomalous dimension ()152|) and 




(195) 
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(a) (b) (c) 

Figure 2: Examples of a 'vertex' correction (a), 'exchange' diagram (b) and self-energy insertion 
(c) contributing to the renormalization of quark- ant iquark operators in Feynman gauge. Path- 
ordered gauge factors are shown by the dashed lines. The set of all diagrams includes possible 
permutations. 

the characteristic polynomials of local conformal operators (Gegenbauer polynomials) which 
form a complete basis and are orthogonal with the weight function u(l — u). We have also 
demonstrated in Sect. 2.4.2 that the orthogonality relation follows from conformal symmetry 
and, in particular, is a direct consequence of special conformal invariance. 

In this section we want to elaborate on the possibility to use such representations, first 
suggested in ^U], and rewrite Eq. (|152|) in an abstract operator form such that the evolu- 
tion equation (jl41|) becomes a spectral problem for a one-dimensional Hamiltonian which is 
manifestly conformally invariant: it commutes with the SL(2, M.) generators and is self-adjoint 
(hermitian) with respect to the conformal scalar product defined in (149)) and (|137|) . This re- 
formulation proves to be very convenient for the analysis of more complicated, three-particle 
evolution equations which we consider at the next step. 

To begin with, let us rederive the ER-BL equation in a somewhat different language, as an 
evolution equation for the light-ray operator built from quark and antiquark fields 

Q(a u a 2 ) = ^(ai)7 + [ai, a 2 ] ip(a 2 ) , (196) 

cf. Eq. ([59)1 . This operator serves as a generating function for local quark- ant iquark twist-two 
operators. In order to avoid complications due to mixing with gluon operators, we will tacitly 
assume that the quark fields have different flavor. The renormalization group equation for the 
operator (|196|) can be written as [3S] 

where M is some integral operator acting on the light-cone coordinates of quark fields. Ex- 
panding the both sides of this equation in powers of a± — a 2 one can reproduce the evolution 
equations for local twist-two operators. One of the reasons why we prefer to work with the 
nonlocal operator (jl96|) is that the conformal symmetry becomes manifest in the coordinate 
space. To the one-loop accuracy the evolution kernel H is given by the Feynman diagrams 
shown in Fig. 121 and can be decomposed in three pieces: 

e = 2n^ 2) - 2?4 12) + i . (198) 
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The first contribution, denoted by liX , describes a gluon exchange between the quark (anti- 
quark) and the Wilson line, as displayed in Fig. It is given in the Feynman gauge by 

f 1 du 

[Hi 12) ■Q]{a 1 ,a 2 ) = - / — (l-u){Q{a u 12 ,a 2 ) + Q{a 1 ,a 2 i l )-2Q{a 1 ,a 2 )} , (199) 

Jo u 

where a^ k = atj(l — u) +akU. The second term, H^ 2 \ corresponds to a gluon exchange between 
the quark and the antiquark, see Fig. |2b- It admits the following representation 

[«< ia > • Q](a x , a 2 ) = f [du] Q« 2 \ a u 2 l) , (200) 
Jo 

where the integration measure is defined as [du] = du\du 2 duz8{ui + u 2 + 1*3 — 1). Finally, 
the extra takes into account the self-energy insertions, Fig. |2fc. Notice that the action 

of the evolution kernel EI on the nonlocal quark-antiquark operator Q(ai,a 2 ) reduces to a 
displacement of the quarks along the light-cone in the direction of each other. 

Expanding the nonlocal operator Q(ai,a 2 ) at short distances over the local composite 
operators 

Q(a 1 ,a 2 ) = J2 E T\f\ Qk ^°^ Q fcl , fc2 (0) = #)) 5^7+5+^(0), 

N>0k!+k 2 =N X " 2 ' 

one finds that the evolution kernel (|198|) preserves the total number of derivatives N and, 
therefore, the evolution equation (| 19 7(1 takes the matrix form. This is a consequence of Poincare 
invariance. For given iV the local operators Qki,k 2 (fy m i x under renormalization and their 
mixing matrix has the size (N + 1). As we have seen in Sect. 3.1, the form of this matrix is 
severely constrained by the conformal symmetry so that only the diagonal entries (anomalous 
dimensions) remain undetermined. 

For the time being, we want to stay with the description in terms of nonlocal operators, 
however. According to Eq. ([18)1 . the action of the collinear conformal group on the two-particle 
operator Q(ai,a 2 ) is given by the differential operators 

L a Q(a u a 2 ) = (L M + L 2ia )Q(a u a 2 ) , (201) 

where a = +, — , and L^ a is the differential operator (JTHJ) acting on the light-cone coordinate 
«fc of the fcth quark, k = 1,2. Using Eqs. (j!99j) and ()200|) it is straightforward to verify the 
following relations: 

[H • L a Q](ai, a 2 ) = L a [M- Q](ai, a 2 ) . (202) 

Hence the evolution kernel H, acting on the Hilbert space of nonlocal light-cone operators 
Q(a\,a 2 ) that admit the short-distance expansion, commutes with the generators L + , L_ and 

[H, L+] = [H, L_] = [H, L Q ] = . (203) 
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This is nothing else but the statement of invariance under the projective transformations on 
the light-cone. Eq. (|203|) implies that EI must be a certain function of the two-particle Casimir 
operator 

L\ 2 = L l + \ ( L + L - + L - L +) = ~d ai d a2 a 2 12 , (204) 

where in the last equality we used that the conformal spins of the quark and the antiquark are 
equal to one for the leading twist, j g = jq = 1, cf. ()63|). 

To find the explicit form of the dependence of EI on L\ 2) it is enough to compare their action 
on the homogeneous translation invariant polynomials a™ 2 = ( a i — a 2) ra with n > 0. These 
polynomials diagonalize the operator L\ 2 and are annihilated by L_ 

L\ 2 a n l2 = (n + 2)(n + l)a» , L.a n l2 = . (205) 

As a consequence, they are the highest weights of the irreducible representation of the conformal 
group, cf. Eqs. (jUjJ and (gnj). Substituting Q(ai,a 2 ) -> a™ 2 , we obtain from Eqs. (HUH), 
and PHD 

?4 12 ) a n l2 = 2 [*(n + 2) - *(2)] < 2 , ^ 12) < 2 = l/[(n + 2)(n + 1)] < 2 , (206) 

where \l/(x) = dlnr(x)/dx is the Euler -^—function. This can be cast in an operator form 
introducing the two-particle spin operator J 12 as a formal solution to the operator relation 

L\ 2 = J 12 (J 12 - 1) . (207) 

From Eq. ()205|) it follows that the eigenvalues of J\ 2 are equal to ji 2 = n + 2 and specify the 
total conformal spin of the quark- antiquark state j\ 2 = n + j q + jq, in agreement with the 
general rule for the sum of two conformal SL(2, R) spins: 

[Ji]®b2] = 0bi+j2 + n]. (208) 

n>0 

Using the expressions in (J205|) and (J206)) we obtain 

= 2 [*(J 12 ) - M>(2)] , Hp) = 1/[J 12 (J 12 - 1)] = 1/L 2 2 ( 209) 

and, finally, 

H = 4 [*(J 12 ) - ^(2)] - 2/[J 12 (J 12 - 1)] + 1 . (210) 

Diagonalization of the renormalization group equation fj!97|) requires the solution of the one- 
dimensional Schrodinger equation 

^H(J 12 ) 0(a lt a 2 ) = 7 0(a lt a 2 ) , 7 = ^ 7 (0) + • • • (211) 

47T 47T 

Since Hi 12 '' and Tii 12 ^ are real functions of the two-particle spin operator J\ 2 which is self-adjoint 
on the SL(2,M) representation space, the evolution Hamiltonian EI inherits the same property. 
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The conformal invariance manifests itself through the fact that the Hamiltonian commutes with 
the generators of the conformal group and, as a consequence, in our case it becomes diagonal 
for the states carrying a definite conformal spin. 

Solutions to (|211j) we already know. They are given by homogeneous polynomials (cki — a 2 ) n , 
see (|206J) . which are nothing else but the characteristic polynomials of the conformal operators 
in the adjoint representation, cf. Eq. pUj) 

H(Jia) P n (ai, a 2 ) = (^/C F ) F n (a h a 2 ) , (212) 

and the anomalous dimension is equal to 7„ = Cp H( J\ 2 — > n+2), in agreement with Eq. (|146j) . 

The major advantage of the representation of the ER-BL evolution kernel in the operator 
form (j210j) is that it is covariant, i.e., retains its form if going over to a different representation 
of the conformal group. For example, one may prefer to have an equation for the characteristic 
polynomials of conformal operators in the "direct" rather than adjoint representation. This 
equation can easily be derived applying F n (d ai , d a2 ) to the both sides of the evolution equation 
(|197j) and making use of Eqs. (J5UJ) and One finds that the characteristic polynomial 

satisfies the Schrodinger equation 

H( J 12 ) P n K u 2 ) = ( 7 i 0) /CV) P n («i, u 2 ) , (213) 

where H is given by the same expression as before, Eq. (|210j) . and the only difference is that 
the two-particle spin operator has to be taken in the adjoint representation, Eq. (}3*2*j) 

MJ12 - 1) = Ll + \{L+L- + L„L + ) = -{d Ul - d U2 f Ul u 2 . (214) 

As before L + = L 1+ + L 2t+ , etc. One can verify by explicit calculation that J 12 and, as a 
consequence, the Hamiltonian M(Ji 2 ) are self-adjoint operators on the Hilbert space endowed 
with the scalar product ()137|) with j 1 = j 2 = 1. As a consequence, the eigenfunctions P 

n\U 1; U 2 ) 

form an orthonormal basis and the corresponding eigenvalues take real values. Notice that 
the form of P n (wi, u 2 ) is uniquely fixed by the conformal invariance 

L? 2 P n (ui,u a ) = (n + 2)(7i+l)P ft (Mi,u 2 ), L_ P„( Wl , u 2 ) = , (215) 

cf. Eq. (I2U5|) . where the second, highest weight condition ensures that the polynomials are 
irreducible with respect to multiplication by a power of {u\ + u 2 ). The solutions to Eq. (J215j) 
are given in Eq. (|47j) . 

Using the completeness condition \n)(n\ = 1 with (uiu 2 \n) ~ F n (ui,u 2 ), one can expand 
the Hamiltonian over its eigenstates. This leads to the ER-BL kernel ()152|) : 

00 

HK, u' 2 ; u x , u 2 ) = u x u 2 J2 PnK, u 2 ) {i^/Cp) F n { Ul , u 2 ) . (216) 

n=0 
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4.2 Baryon Distribution Amplitudes 

The Hamiltonian formalism is very powerful for the analysis of three-parton systems. We will 
concentrate on two examples - the nucleon distribution amplitude, 0jv(wj), and the distribution 
amplitude of the A— resonance state with helicity— 3/2, (p 3 ^ 2 (v,i) , both are leading twist. They 
are defined |8(J HH] through the matrix elements of nonlocal three-quark operators of definite 
chirality (gtW = | (1 ± 75) q): 



-1K P+ ^Art(p) [du}e- ip +^ a ^M^) (217) 



and 



Ca, + u}(a 2 ))(A^) uI(as)|A(p)> 



= z/l /2 p+Al f(A^ u J [du] e^^^f^) . (218) 

Here are color indices, C is the charge conjugation matrix, \P(p)) and |A(p)) are the 

proton and A— resonance state described by the spin— 1/2 spinor N(p) and spin— 3/2 Rarita- 
Schwinger field A£(p), respectively. 

The scale dependence of 0^ 2 (uj, fi 2 ) and </>7v(uj, /i 2 ) is governed by the renormalization group 
equations for the three-quark operators with three spinor indices |37I] . 

B 3 £( ai ,a 2i a 3 ) = e^\M}MUM}MW4M)c, (219) 
B^ax, a 2 ,a 3 ) = ^ fc (M T («i))a(^j(« 2 )) fe (#I(a 3 )) c , (220) 

respectively. The superscript 1/2 or 3/2 denotes helicity. The nonlocal operators (|219|) and 
(I220|) have to be interpreted as generating functions for the renormalized local operators 

fel ^2 ^3 

B( ai ,a 2 ,a 3 ) = J2 E ^^^(nD) k ^(nD) k2 m(nD) k ^(0) , (221) 

N>0 k 1 +k 2 +k 3 =N l ' 2 ' 3 ' 

whose reduced matrix elements define moments J[du] u kl u 2 2 u 3 3 (j)(ui) of the distribution ampli- 
tudes. 

All three-quark local operators in (j221|) with the same total number of derivatives N = 
k\ + k 2 + k 3 get mixed by the evolution, so that for given N the size of the mixing matrix is 
N(N + l)/2. As before, conformal symmetry allows one to resolve the mixing with operators 
containing total derivatives. This is done in the standard way, introducing the basis of three- 
particle conformal operators 

O^(0)=¥ N>q (d ai ,d aa ,d a!l )B(a 1 ,a 2 ,a 3 ) _, (222) 



a 4 =0 
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which are enumerated by the subscript q and satisfy the defining equations 



L 2 F N>q ( Ul , u 2 , u 3 ) = (N + 3){N + 2)F N>q ( Ul , u 2 , u 3 ) , L_ F N , q ( Ul , u 2 , u 3 ) = . (223) 

Here L 2 is the three-particle Casimir operator in the adjoint representation. In the present case 
(ji = j 2 = j3 = 1) it can be written as a sum of two-particle Casimir operators L 2 = L\ 2 + 

+ L\ x . The highest weight condition L_ Pat j(? = can be replaced by L Q Pat j(? = (N + 3)Pjv, 9 . 

Note that three-particle conformal operators with different conformal spin N + 3 are or- 
thogonal with respect to the conformal scalar product 

(F N JF NW )= [ [dw] U f- 1 «f- 1 «f- 1 P JVY (u 1 ,« 2) W3)Piv ig («i,W2,W3)~ <W- (224) 
J 

This orthogonality relation can easily be established generalizing the analysis in Sect. 2.4.2 to 
the correlation function of two three-particle conformal operators. 

In difference to the situation that one encounters in the analysis of leading-twist meson 
distribution amplitudes, the constraints (|223|) (or, equivalently, the orthogonality relation (|224j) ) 
do not fix the characteristic polynomials of conformal operators uniquely, cf. Sect. 2.2 and 
Eqs. (I5U|) . (|5ip. Hence we need an extra index q to distinguish between different conformal 
operators with the same conformal spin j — N + ji + j 2 + j 3 = N + 3 alias the same number 
of derivatives N, cf. Eq. (J51|) . Existence of multiple three-particle conformal operators with 
the same spin can easily be understood as the multiplicity of the corresponding irreducible 
component in the sum of three quark spins 

[Ji] ® [J 2 ] ® [Js] = b'l + 32 + Js + N] . (225) 

Af>0 

Applying the summation rule for two conformal spins in Eq. (|208|) . [ji] © \j 2 ] = ^2\jx 2 = 
ji + 32 + n], we find that the total spin j = N + 3 component has a nontrivial multiplicity 
(7V+1) corresponding to possible values of the two-particle spin < j% 2 < N, in agreement with 
our discussion in Sect. 2.2. As the result, for a given N the conformal constraints (J223|) generate 
(N + 1) linear independent solutions for the characteristic polynomials F Ntg (ui,u 2 ,u 3 ). The 
corresponding conformal operators have the same spin N + 3 and, therefore, nothing prevents 
them from mixing with each other. It follows that the mixing matrix is not diagonalized but 
its size is reduced to (N + 1). In other words, the leading twist baryon distribution amplitudes 
have to be written as a double sum 80 



<f>(ui, u 2 , u 3 ; fx 2 ) = 120wiw 2 w 3 0jv, 9 (a*o) ^N, g (ui, U2, «a) ( —, — r ) , (226) 

where ¥^ jQ (ui, u 2 , U3) are certain polynomials that are constrained, but not fixed uniquely by 
Eqs. f|223J) (or, equivalently, by the orthogonality relation (1224(0 and 7^ are the corresponding 
anomalous dimensions. The prefactor 120u\u 2 u 2 in Eq. ()226|) is the asymptotic distribution 
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<g>-----®-----® <£>--- 

a b c 

Figure 3: Examples of a 'vertex' correction (a), 'exchange' diagram (b) and self-energy insertion 
(c) contributing to the renormalization of three-quark operators in Feynman gauge. Path- 
ordered gauge factors are shown by the dashed lines. The set of all diagrams includes possible 
permutations. 

amplitude, cf. Eq. (|159|1 . It is nothing else but the weight function in the conformal scalar 
product and is fixed by the symmetry. On the other hand, jN >q and Pjv,g(wi) have to be found by 
the explicit diagonalization of the (reduced) mixing matrix [HI] , where mixing with the operators 
containing total derivatives can be ignored, i.e., by taking forward matrix elements. Progressing 
to large values of iV one obtains rather complicated matrices which are not symmetric and 
have no obvious structure. For many years the only method to study the evolution of baryon 
distribution amplitudes was based on numerical evaluation of such matrices written in terms 
of contributions of symmetrized Appell's polynomials, see e.g., [S2], which correspond to one 
possible choice (15 1|) of the three-particle conformal operators. It was believed that the role 
of conformal symmetry is exhausted by this representation and the problem has no further 
symmetry; hence, no analytic approach seemed feasible. 

We will demonstrate that this pessimistic conclusion is not warranted. It turns out that all 
important features of the spectrum of anomalous dimensions and eigenfunctions of three-quark 
states can be completely understood in analytic form. The first step is to find an orthonormal 
basis of conformal operators, such that any two operators with different spin iV and/or different 
q are orthogonal with respect to the scalar product ()224j) . (F N)q \F Ntq i) = N5 qq i. Expanding the 
solutions of the evolution equations in this basis, one obtains the reduced mixing matrix 

/ full mixing matrix of 
\ composite operators 

V v 

N(N + l)/2 

which has the same size N + 1 as in the standard approach, but is symmetric. This symmetry 
is a direct consequence of the underlying SL(2,M.) symmetry of the evolution equations and 
explains, e.g., why all anomalous dimensions of three-quark operators are real numbers. Most 
importantly, in this form new (hidden) symmetries can be observed and we will find that the 
evolution equation for the distribution amplitude of the A— resonance is completely integrable 
[Ulinni- Last but not least, for symmetric matrices one can use a simple perturbation theory 
to take into account contributions that are suppressed by a certain parameter, and this proves 
to be crucial in many further applications. 

We begin now with the explicit construction. The nonlocal operators B 3 / 2 f|219J) and B 1 / 2 



hermitian kernels 
in conformal basis 



(227) 



N + l 
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(|22(Jj) do not mix with each other since they belong to different representations of the Lorentz 
group. To one-loop accuracy the corresponding renormalization group equation can be written 

as [TDl EE] 

j/i ^ + P{g) ^} B\a u a 2 , a 3 ) = ^ [H A • £? A ] K a 2 , a 3 ), (228) 

where the evolution kernel Ha depends on the helicity of the three-quark operator and corre- 
sponds to contributions of the Feynman diagrams shown in Fig. El Apart from the different 
color factor these diagrams coincide with those defining the two-particle evolution kernels 7i v 
(PU) and H c (HP), so that 

H A = (1 + l/N c )H x + 3CV/2 , (229) 

where the last term is due to the self-energies, Fig. Efc, the color factor follows from the identity 
e ijktii'tjj' — — £ i'j'k(l + l/N c )/2 and the Hamiltonians H\ are given by the sum of two-particle 
kernels 

n x=3/2 = + n w + ) Hx=i/2 = Hx=3/2 _ H (i2) _ ^ (230) 

where the superscripts are introduced to show that Tiy® and Tie act on the light-cone coor- 
dinates of the quarks l V and '/c'. Note that the difference between the two Hamiltonians for 
A = 3/2 and A = 1/2 in Eqs. ()230|) is due to the fact that the gluon exchange diagram in Fig.Eb 
vanishes unless the two quarks have opposite chirality. 

Since the two-particle kernels 7{y and He k \ by construction, only depend on the corre- 
sponding spin operators the Hamiltonians 7i\ are explicitly SL(2, M) invariant and commute 
with the three-particle generators L a = L l a + L 2a + L 3 a (with a = +, — , 0): 

[H x , L 2 ] = [H x , L + ] = [H x , L_] = [H x , L ] = . (231) 



Finding solution to the renormalization group equation ()228|) is equivalent to solving the 
Schrodinger equation for the characteristic polynomials of multiplicatively renormalizable op- 
erators 

H x { Jia, ^23, J13) • P^,,(mi, w 2 , u 3 ) = E^ q F x N>q ( Ul , u 2 , u 3 ) , (232) 

and the anomalous dimensions are given in terms of the "energies" as 

7^ = (1 + l/N c )E x N>q + 3C F /2 . (233) 

Since the Hamiltonian commutes with the three-particle SL(2, M) generators, we enforce addi- 
tional conditions (J223)) on the eigenfunctions, i.e., require that they are characteristic polynomi- 
als of conformal operators. In addition, since the Hamiltonian is hermitian, all its eigenvalues 
E^ q are real numbers and the eigenfunctions corresponding to different E^ q1 alias different q, 
have to be orthogonal with respect to the conformal scalar product 



56 



The Schrodinder equation (|232|) is written in operator form, and in order to obtain the 
general solution it is convenient to expand the eigenfunctions IPV,g(wi) over an orthonormal 
basis of functions satisfying the same conformal constraints, Eq. (|223jl . A suitable basis can 
be constructed as follows [TTJ |3U1 EZj • We define the set of polynomials Y^ 3 (m 1 , m 2 , u 3 ) with 
n = 0, 1, . . . , N by requiring that in addition to ()223|) they diagonalize the two-particle Casimir 
operator for the (12)-quark pair (j214j) . L\ 2 = J\ 2 (J\2 _ 1) : 

Ll 2 Y { if(u i )= 3l2 Ui2-l)Y { if(u i ), (j 12 = n + 2, < n < N) . (234) 



The solutions to the combined equations ()223|) and ()214|) are given by products of Gegenbauer 
and Jacobi polynomials and are presented for the general case of arbitrary conformal spins 
Ji,j2)j3 in Appendix El For our present purposes it turns out to be convenient to use a 
different normalization, so we choose 

<f = \(N + n + 4)(n + 2) Y^ + ^ =n+2 { Ul ) , (235) 

where the functions Yj- (wj) are given in Eq. (TO) . The funct ions Y^J^(ui) depend on the 
pair of integers N and n which are related to the total conformal spin of the three-quark operator 
ji +J2+J3 + N = 3 + N and the conformal spin of the (12)— pair j 12 = ji +J2 + n — 2 + n with 
ji — H = J3 = 1- 111 what follows we often drop the subscript W if it is clear from the context. 
The particular choice of a quark pair in Eq. (J234|) is of course arbitrary and we might use, e.g., 
L23 to construct a different basis of functions Yn 23 '(«g). Here, the superscript indicates the 
order in which the tensor product of three SL(2) representations has been decomposed into the 
irreducible components. The two set of basis functions Yn 23 ^^ and Y n 12 ^ 3 (ui) are related to 
each other through the Racah 6j— symbols of the SL(2) group (see Appendix B in Ref. [H|). 

One convenient property of the functions Y n 12 ' )3 (uj) is that the action of the Casimir op- 
erators on them turns out to be rather simple. By construction, Y]^ [ui) diagonalize the 
total three-particle Casimir operator L 2 and the Casimir operator for the (12)— quark pair L\ 2l 
whereas the remaining two two-particle Casimir operators turn out to be three-diagonal: 



^ 3 Yi 12)3 (^) = /„ 
^iY n 12 ^) = /„ 



(n + 1) n - L v 11 (n + 2)(n + l) n y lJ (n + 2) 



{n + 1) ( n + 2)(n + l) n v ' (n + 2 



(236) 



where f n is defined in Eq. f|239j) below. This property turns out to be crucial for simplification 
of the evolution equation. 

The most important property, however, is that Y^ 2 J | 3 (-Uj) form an orthogonal basis on the 
space of characteristic polynomials 

(YS 3 |Y^S)~W5nn', (237) 
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the reason being that L\ 2 is a self-adjoint operator with respect to the scalar product (|224|) 
so that its eigenstates with different n have to be mutually orthogonal. As the result, the 
Hamiltonians 7Ya=3/2 and 7Ya=i/2 (I230|) can be written in the conformal basis as 

m kn = (Yl i23) \H x \YS *>)/(\\Y? 3) || ■ ||Yf 3 )||) , (238) 

where [H\] nk = [7~L\\k n are hermitian (symmetric) matrices of the size N + 1 which can be 
expressed in terms of the Racah 6j— symbols of the SL(2, M.) group. 

We look now for the solutions Pjv, 9 of the Schrodinger equation ()232|) . expanding them in 
this basis 

F N , q ( Ulj u 2 ,u 3 ) = ^ C -M Y ^ k 3 (u u u 2j u 3 ) , f k = {k ^l^ ] (N-k+l)(N+kU) , (239) 

where Ck{q) are the expansion coefficients and the factor i h / fk was inserted for later convenience. 
This gives a matrix equation for c n (q) 

N 

E %n ' k c ^)lf* = E k q c n{q)lh , (240) 

k=0 

where it is tacitly assumed that c n (q) also depend on the total spin N and helicity A. The 
parameter q enumerates different solutions to Eq. ()240|) . As was already mentioned, for given 
iV one expects to find N + 1 different eigenstates. Using explicit expressions for the matrices 
[H\] kn is becomes straightforward to solve Eq. (|24U|) and obtain the spectrum of anomalous 
dimensions r y^~^ 2 and 7^~ 3 ^ 2 for fixed N. The results of numerical calculations for < iV < 30 
are shown in Fig. |3J We notice that the both spectra exhibit remarkable regularity. The 
spectrum of H1/2 is very similar to that of H3/2 except the few lowest eigenvalues for each N, 
which are separated from the rest of the spectrum by a gap that remains finite at large N. As 
we will show in the next Section, these properties are a manifestation of a hidden symmetry of 
the Hamiltonian 7^3/2- 

4.3 Complete Integrability 

As was explained in the previous Section, the QCD evolution equations for three-quark op- 
erators have the form of the Schrodinger equations (|232j) describing a three-particle system 
with three degrees of freedom corresponding to quark momentum fractions Mj. In this way, the 
scale dependence of baryon distribution amplitudes in QCD corresponds to a one-dimensional 
quantum mechanical 3-body problem with very peculiar Hamiltonians, (j229J) . ()230|) and ()209|) . 
determined by the underlying QCD dynamics. The conformal symmetry allows us to trade two 
degrees of freedom for two quantum numbers corresponding to the total conformal spin L 2 and 
its projection L , Eq. (J223|) . The remaining degree of freedom is described in the conformal 
basis by the set of coefficients c n introduced in Eqs. ()239|) and ()240|) . Thus, thanks to confor- 
mal invariance, the original 3-body Schrodinger equation ()232|) is reduced to a (complicated) 
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Figure 4: The spectrum of eigenvalues for the helicity-3/2 (left) and helicity-1/2 (right) Hamil- 
tonians Tis/^N) and 7ii/ 2 (iV), respectively. The lines of the largest and the smallest E 3 / 2 are 
indicated on the plot for E\/i by dots for comparison. The solid curves in the left panel show 
three different trajectories for the anomalous dimensions (see text and Fig. EJ). 



one-body problem (|240|) which is not possible to solve analytically for arbitrary N, unless it 
has some additional symmetry. 

It turns out that the Hamiltonian (but not H1/2) possesses such additional 'hidden' 
symmetry. Namely, one can construct the operator 

, — ■ , — ■ 

Q = 2 t L i2' L 2a] = i (dm - d U2 ) (d U2 - d U3 ) (d us - d Ul ) u^Ua , (241) 
that commutes with 7^3/2 and with the SL(2) generators: 

[Q,L a ]=0, [Q,H 3/2 ] = 0. (242) 

Thus, the operator Q is the integral of motion and its eigenvalues q specify uniquely the 
spectrum of the anomalous dimensions of the three-quark operators of helicity A = 3/2. The 
first relation in Eq. ()242j) follows readily from the definition of the charge Q. The second relation 
can be checked using the commutation relations between Q and two-particle Hamiltonians, 
Eqs. (ESDI) and flU} 

[Hi 12 \ Q] = z(Z^ 3 - Ll) , [W? 8 ), Q] = i{Ll x - L\ 2 ) , [H^\Q] = i{Z{ 2 - . (243) 

To prove these operator identities, apply the both sides to the function Y^ 12 -* 3 ^) belonging 
to the conformal basis, use Eqs. ()236|) . and take into account that the operator Q, Eq. (j241D . 
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is represented in the conformal basis by a (N + 1) x (TV + 1) matrix with only two nonzero 
subleading diagonals 



(244) 



Existence of a nontrivial integral of motion Q implies that for A = 3/2 the Schrodinger equation 
()240|) is completely integrable. This allows us to calculate the corresponding spectrum of the 
anomalous dimensions analytically by applying a powerful technique of integrable models. 

Remarkably enough the Hamiltonian H3/2 is well-known in the theory of integrable models 
as a generalization of the celebrated Heisenberg spin magnet describing the nearest-neighbor 
interaction between spin— 1/2 operators (Pauli matrices) with the Hamiltonian 



This model has been studied in 1931 by Bethe as a model of one-dimensional metal and solved 
by the method which has become well-known as the Bethe Ansatz. The Hamiltonian ()245|) 
possesses the set of integrals of motion and the corresponding Schrodinger equation is completely 
integrable for arbitrary number of spins. Later it was recognized that the Heisenberg spin 
magnet can be generalized to spins belonging to higher-dimensional representations of the 
SU(2) and the SL(2,M.) groups [SHI- Although the Hamiltonian becomes much more involved, 
it inherits the property of complete integrability and its spectrum can be studied by the Bethe 
Ansatz technique. 

In application to the QCD evolution equations, one encounters the integrable generalization 
of the Heisenberg magnet for the SL(2,M.) spins. Indeed, an inspection shows that the SL(2) 
generators (|32|) for quarks with conformal spin jk — 1 can be interpreted as the usual spin 
s — — 1 operators. In this way, we may consider the Hamiltonian 7^3/2 as describing the 
system of three interacting spins each acting on its internal space labeled by the coordinates 
Mfc. These spins carry the index of the corresponding particles and form a one-dimensional 
spin chain with three sites which coincides identically with the one-dimensional homogeneous 
Heisenberg spin magnet of noncompact spin s — — 1. This model has been solved using the 
Bethe Ansatz technique and the detailed results can be found in Refs. [84J. In the rest of this 
section we explain the exact solution for the Hamiltonian H3/2 and discuss the main features 
of the spectrum. 

Complete integrability of the Schrodinger equation for A = 3/2 implies that H3/2 is a 
(complicated) function of two and only two mutually commuting operators Q and L 2 . Therefore, 
instead of solving the Schrodinger equation (|232|) directly, one can solve much simpler equations 
(I223j) supplemented by the additional constraint 



QPjv,,(tti) = i (d Ul - d U2 ) (d U2 - d U3 ) (d U3 - d ui ) u x u 2 u z ¥ N>q (ui) = q¥ N ^(ui) (246) 



and find the spectrum of the Hamiltonian 7^3/2 = 7^3/2 (L 2 , Q) by replacing the operators by 
their corresponding eigenvalues. Expanding Pjv,q(wi) over the conformal basis, Eq. ()239|) . and 




(245) 



n 
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making use of Eq. ()244|) , one finds that the equation Q P = q P is equivalent to the following 
three-terms recurrence relation for the coefficients c n , (n = 0, . . . , N): 

qc n = f n (c n+ i + c n _i) , c_i = c N+ i = 0. (247) 

Since the operator Q is self-adjoint with respect to the scalar product (|224|) . its eigenstates 
Piv,q(wj) are orthogonal to each other for different set of quantum numbers. Their norm depends 
on the overall normalization of c n . It is convenient to choose it in such a way that the eigenstates 
Pjv i9 (iti) form an orthonormal set of the states 

(PjV,gFJV',g') = $NN'8qq' • (248) 

The recurrence relations (|247j) represent the system of N+l linear homogeneous equations on 
the coefficients c n . Solution of this system is equivalent to diagonalization of a (N + l) x (N + l) 
matrix with only two subleading diagonals nonzero. The consistency condition for this system 
leads to iV+ 1 quantized values of q, which have the properties of roots of orthogonal polynomi- 
als. That is, they are real, simple and for different N the sets of quantized q are interlaced, see 
Fig. El The set of coefficients c n = c n (N, q) defines uniquely the characteristic polynomial of the 
three-quark operator, Eq. ()239j) . The corresponding exact value of the anomalous dimension is 
given by 

eT = 4 Re £^V„W + 2)-*(2)] gn+jg 

Elo P" + 3) 



The recurrence relations ()247|) combined with the expression for the anomalous dimensions 
and the characteristic polynomials ()239|) allow us to determine the spectrum of local conformal 
three-quark operators of helicity— 3/2, (|222|) . One finds from (|247|) that c n (—q) = (—l) n c n (q) 
and, a consequence, the energy levels corresponding to nonzero values of quantized q are double 
degenerate 

Ef q = Ef_ q . (250) 
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The resulting spectrum of the conserved charge q and the "energy" for < iV < 30 is 
shown in Fig. and Fig. HI respectively. 

As follows from (J218|) . conformal moments of the distribution amplitude of the A-resonance 
with helicity A = 3/2 are related to the (reduced) matrix elements of conformal three-quark 
operators 



[du\ (l>3/2(Ul,U2,U3;iJ 1 )f'N,q(Ul,U 2 ,U 3 ) = (OlOj^/io)!-?) ( ^Jj^j) ■ ( 251 ) 



3/2 :a 



Thanks to the orthogonality condition ()248|) this relation can be inverted to obtain the desired 
expansion of the distribution amplitude envisaged in Eq. (J226|) . 

For fixed N the dominant contribution to the sum in ()226|) comes from the conformal 
operators with the lowest anomalous dimension. For even N such operators are those with 
q = 0. The recurrence relation (|247|) can be easily solved in this case and leads to c 2n = (— l) n c 
and c 2 „+i = 0. Using ()249|) then yields 

E 3 J 2 q=0 = 4^(iV + 3) + A lE - 6 . (252) 

The corresponding curve is shown in Fig. 0] (left) by dots. The corresponding eigenfunction 
(I239J) is a completely symmetric real function of Ui and is given (up to an overall normalization 
factor) by [SU] 



Ul (l - ujcf^l - 2u x ) + u 2 (l - « 2 )cKi(l - 2m 2 ) + u 3 (l - w 3 )CaSi(1 - 2u 3 J?53) 



UiU 2 U 3 PjV, 9 =o(«l, «2, U 3 ) 

= m(i - ui)C%j 

where U\ + u 2 + u 3 = 1. 

The expressions for E^ q and Pjv, 9 for the general case q ^ cannot be written in a closed 
form. A significant simplification occurs for large N, however. In this limit one can solve the 
recurrence relations ()247|) by the WKB method and develop the systematical expansion of the 
spectrum in powers of 1/N. One finds that for iV ^> 1 the energies and conserved charges q 
occupy the bands 

4 In iV < E 3 J 2 q < 6 In iV , 0<g 2 <iV 3 /27. (254) 
The dependence of E^ 2 q on the charge q (the "dispersion curve" ) takes the form 

3 

E 3 J 2 q = 2 In 2 - 6 + 6 7E + 2Re ^i 1 + iX k) + 0{1/N 6 ) (255) 

k=l 

where 7e = 0.57722 ... is the Euler number and are defined by the three roots of the cubic 
equation 2A 3 — A^ — q — 0. 

Possible values of the charge q satisfy the quantization conditions 

giV- 2 lniV-argr(l + iqN' 2 ) +OCL/N) = %(N -21), (256) 

6 
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where £ is a nonnegative integer. Solving (J256)) one obtains the quantized values of the charge 
q = q(N,£). For given £ the values of q form a trajectory which depends analytically on N. 
Each trajectory q(£, N) is mapped into the corresponding trajectory for E^ 2 e as shown in Fig. 0] 
The £— th trajectory starts at iV = £, approaches the 'Fermi surface' (|252|) at N = 2£, gets 
repelled from it and monotonously grows to infinity at large N. The corresponding asymptotic 
expression for the energy reads (with rj = ^J{N + 3)(N + 2)) 



Ef t = 61n *T1 - 3(2? + 11 - g°g + f - 7 - 464P + 696P = 8 ° 2f - 517 +Q( W) (257) 



at large iV > I, and [HI 130] 



= 41n(iV + 3) + 4 7E - 6 + -f0^~ A N - 2£f (258) 

18 In (77 e^ E ) 

in the vicinity of N = 2£. These relations define the asymptotic expansion for the energy levels 
of the Hamiltonian parameterized by the integer £ in the upper and the lower parts of the 
band (|254|) . respectively. The corresponding level spacing is equal to 

5E$ 9 = O (1/ln 2 N) , 5E$ Q ^ N ^ O (1/N) . (259) 

To summarize, the evolution equation for the distribution amplitude is exactly solv- 
able. The physical interpretation of integrability is that we are able to identify a new 'hidden' 
quantum number q, which distinguishes components in the A-resonance with different scale de- 
pendence, Eq. (I226J) . In this case, the coefficients c n in the expansion ()239|) of the eigenfunctions 
of the evolution equation over the complete set of orthogonal conformal polynomials can be 
calculated using a simple three-term recurrence relation (J247|) . The corresponding anomalous 
dimensions are given in terms of c n by Eq. ()249|) . For large iV the spectrum is described to a 
high accuracy by the WKB formulae ()257|) and (|258|) . 

The most interesting result concerns the structure of the component of (f) 3 ^ 2 with the smallest 
anomalous dimension for each N. It represents the leading contribution to the distribution 
amplitude ()226|) in the limit fi 2 — > 00. The corresponding polynomials are known exactly and 
are given in Eq. ()253|) (see Fig. EJ). The physical interpretation of such 'ground states' is most 
transparent in coordinate space. Up to contribution of operators with total derivatives, one can 
represent the three-quark 'ground state' in a concise form as the nonlocal light-cone operator 

m 

if 1 

B^piax, a 2 , a 3 ) = - ^ e ljk / dv M]{a a ) ^]{va a + (1 - v)a b ) M(a b ) . (260) 

^ h 1 9 "i ^ 



a, 6=1, 2, 3 
ajtb 



The Taylor expansion at short distances, a i2 , a 32 — > 0, generates the series of local multiplica- 
tively renormalizable three-quark local operators with the lowest anomalous dimension for each 
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Figure 6: Contributions to the A = 3/2 distribution amplitude 4>a i u i) with lowest anomalous 
dimensions for N = 2 and N = 4. The normalization is arbitrary. 



even TV: 



B { 3 q /2 \a>i,a 2 ,a>3) = ^ 



Ot 12 + «23 + a : 

(iV + 1)! 




N,q=0 



(di,d 2 ,d 3 ) B(an, a 2 ,a 3 ) 



+ ... . 



a l= a 2=°'3 =0 

a 1 +d 2 +d 3 =o 



(261) 



Note the integration in Eq. ()26()j) with unit weight over the position of the quark in the middle 
that goes in between the light-cone positions of the other two quarks, up to permutations. If 
renormalization of the operator is interpreted as interaction, this unit weight can in turn be 
interpreted as the statement that the quark in the middle is effectively 'free': In the 'ground 
state' with the lowest 'energy', the interaction of the quark in the middle with its right and left 
neighbors exactly compensate each other. 

4.4 Further Developments 

Let us now consider the evolution equation for helicity— 1/2 distribution amplitude <f>N(ui). 
We recall that the corresponding Hamiltonian differs from that for (J) 3 / 2 by the additional 
contributions of gluon exchange between the quarks with opposite helicity, see Eq. (I230J) and 



The added terms destroy the complete integrability and hence the Schrodinder equation for 
Ti.i/2 cannot be exactly solved. The numerical calculation for < N < 30 gives the spectrum 
in Fig. EJ As seen from this figure, the spectra of Ti.i/2 and 7^3/2 are very similar in the upper 
part, for larger eigenvalues, and at the same time the two lowest levels of the 7ii/2 Hamiltonian 
appear to be special and 'dive' considerably below the line of lowest eigenvalues of 7^3/2, given 



Fig. El 



Hi/2 = H3/2 + V 



V = -1/4, - l/Ll 



(262) 



by Eq. (J2S2J). 
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The comparison of the spectra in Fig. 0] suggests to consider the operator V in Eq. (J262)) 
as a perturbation. As usual in quantum mechanics, for perturbation theory to be applicable 
the perturbation should be smaller than the level spacing. To verify whether this is the case, 
one has to examine the matrix elements of V over the exact eigenstates of H3/2. The explicit 
calculation shows that the matrix (Pjv ig |V|Pjv l(? ') is strongly peaked at the diagonal q = q' and 
decreases rapidly with \q — q'\. In the lower (q — ► 0) and the upper (q — > N 3 /^/27) part of the 
spectrum it scales for large N as 



(^,n\VKn) q ' Q =° O (1/lnJV) , (F q , tN \V\F qtN ) q ^= /v ^ O (l/N 2 ) , (263) 

respectively. Comparing ()263|) with ()259j) we conclude that the perturbation theory is justified 
for large N for the upper part of the spectrum. A simple calculation then leads to 



\ , 2(21+1) 
V 



+ 0{l/rf) , (264) 



where 77 = ^/(N + 3)(N + 2). 

On the contrary, in the lower part of the spectrum several (of order ~ In N) lowest energy 
eigenstates are affected strongly and their mixing cannot treated perturbatively. It turns out 
that the effective Hamiltonian describing the interactions between the lowest energy levels can 
be brought to the form of (generalized) Kroning-Penney model of a single particle in a periodic 
6— function potential. Using the well-known solution to this model, one can show that the two 

1 /2 

lowest eigenvalues (hence also the anomalous dimensions) E± (N) decouple from the rest of 
the spectrum and in the limit iV — ► 00 are separated from the other eigenvalues by a finite 
constant — the 'mass gap', 

4 /2 (iV)-< /2 9=0 ^-0.3, (265) 

where E^{ 2 q=0 is defined in Eq. (|252|) . The corresponding contributions to the distribution 
amplitude are given for iV — > 00 by 

Ux^aPjf^K^s) = u lU2 u 3 [P^ 1] (2u 3 - l)±P$ 1) {2u 1 - 1)] , (266) 

where P^' are the Jacobi polynomials. Notice that the expression in the square brackets 
does not depend on the momentum fraction of the quark with helicity opposite to that of 
the parent baryon. As a consequence, going over from the conformal three-quark operators 
defined by the characteristic polynomials ()266|) to the generating nonlocal light-cone operator 
B(ax, CH2, oe 3 ), one finds that Eq. (|266|) corresponds, in the same sense as Eq. (|26U|) . to the 
following contribution 



B ± (a 1 ,a 2 ,a 3 ) 



£ ijk 



(M T Mt( a 3)<Ka 2 - «i)± ^J(ai) Wj) MD( a 3)5(a 2 - «s) • (267) 
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Formation of the mass gap in the spectrum of anomalous dimensions is, therefore, naturally 
interpreted as due to binding of the quarks with opposite helicity into scalar diquarks. 

Note that while the expression ()253|) for the eigenfunction is exact, the result ()266|) is only 
valid in the asymptotic In N — ► oo limit. In the coordinate space picture, the restriction to large 
N is translated to the condition that the light-cone separation between the same helicity quarks 
is very large to allow for the formation of a diquark. In the momentum space, the result means 
that at sufficiently large normalization scale /i 2 the quark carrying a very large momentum 
fraction is more often with the same helicity as of the parent baryon. This observation seems to 
be in qualitative agreement with phenomenological models of baryon distribution amplitudes 
derived from QCD sum rules |BT) ISB} IHTj. 

We have demonstrated that the Hamiltonian approach offers a convenient technique for 
finding the scale dependence of the baryon distribution amplitudes. The power of this approach 
is largely due to the fact that it is conformally covariant, with the Hamiltonians being self- 
adjoint operators on the space of characteristic polynomials endowed with the conformal scalar 
product ([224)1 . They only depend on the two-particle conformal spins but the explicit form of 
the dependence is not fixed by the symmetry. Finally, the Hamiltonian 7i 3 /2 possesses a hidden 
symmetry which allows one to solve the corresponding Schrodinger equation exactly and identify 
a new quantum number q that distinguishes different components of the helicity— 3/2 baryon 
distribution amplitude. Although the Hamiltonian H.1/2 does not possess the same symmetry 
due to presence of additional terms, its complete integrability is broken "softly". Its energy 
spectrum has many features common with the spectrum of 7^3/2 and only a few lowest energy 
levels are affected. We would like to stress that integrability of 7i 3 /2 is intrinsically tied to a 
peculiar feature of anomalous dimensions in gauge theories, which involve Euler ^/-functions 
and are rising logarithmically with the spin (dimension) [88J. For large J12 one finds from ()209|) 
that the both Hamiltonians H3/2 and H1/2 have the same universal form 

Hx = 2 In J 12 + 2 In J 23 + 2 In J 31 + 0(1/ Jf 2 , 1/ J 2 3 , . (268) 

The integrability of H3/2 imposes severe restrictions on the form of subleading in l/Jik terms. 

The same approach can be used for studies of the evolution of quark- ant iquark-gluon and 
three-gluon operators. These are interesting, e.g., in connection with the structure function 
92(2, Q 2 ) of polarized deep inelastic scattering, which is attracting increasing interest as it 
provides a direct measurement of three-parton correlations. The scale-dependence of g2(x,Q 2 ) 
can be traced to the renormalization-group equations for the quark-antiquark-gluon operator 

[HEllEni 

S ± (a 1 ,a 2 ,a 3 ) = -^(a^igG x+ {oi 2 ) ± gG ±+ (a 2 )^} ?#(a 3 ) , (269) 

where G^ v = \s^ up \G a ' px t a , G^ v = G^t 0, . In the flavor singlet sector one also has to include 
the three-gluon operator 

5(a u a 2 , a 3 ) = l j-f abc Gl x ( ai )G b +± (a 2 )G c +x (a 3 ) (270) 
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which mixes with the C-even part of the flavor singlet qGq operator (|269|) . 

To begin with, consider the flavor non-singlet sector. In this case the three-gluon operator 
drops out and the evolution of S simplifies dramatically in the limit of large number of colors, 
N c — > oo. The corresponding Hamiltonian reads, e.g., for S + JT] 



H 5+ = N c 



V qG (J 12 ) + U Gq (J 23 ) +0(1/N, 



(271) 



where 



v gG (j) = + 

u Gq (J) = «r[j + i] + *(j-i 



2^(1) 
2V(1) 



3 
4 
3 
4 



(272) 



and it turns out to be equivalent to the Hamiltonian of the completely integrable open Heisen- 
berg magnet [TT]. The corresponding conserved charge is equal to 



{L 



2 r2 
12; ^ 



1 

2* 



23} 



9 L 2 



23 > 
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where {, } stands for an anticommutator. Properties of this Hamiltonian were studied in detail 
in Refs. [SHI EDI EI] and the corrections 0{1/N 2 ) to the spectrum were calculated in Ref. 
The most interesting result is that the lowest energy, albeit the lowest anomalous dimension 
in the spectrum of quark-gluon operators can be exactly found, and the corresponding eigen- 
function coincides to leading logarithmic accuracy with the structure function g 2 itself. Thus, 
a simple evolution equation arises |H5] . 

The study of the flavor singlet sector is much more challenging [H7j. Ignoring quarks for 
a moment, the evolution kernel of three-gluon operators becomes the Hamiltonian of closed 
Heisenberg magnet with three "gluonic" sites of the spin j g = 3/2, whose integrability is 
"softly" broken by additional terms |91j . This can be studied along the same lines as the three- 
quark operators, although the algebra becomes more involved because of higher spins. The most 
interesting part appears to be the mixing between quark-gluon and three-gluon operators which 
for large N can be interpreted as describing the interaction between open and closed Heisenberg 
magnets. It turns out that this mixing has rather peculiar features, which we cannot discuss 
in detail in this review, but the outcome is that the mixing can to a large extent be reduced to 
a few participating levels. Identifying important degrees of freedom (at least for large N) one 
can write down an approximate two-channel evolution equation for the structure function g 2 in 
terms of transverse spin densities |37|. This approximate evolution equation looks as follows. 

To leading logarithmic accuracy, the structure function g 2 (x,Q 2 ) is expressed as 



S2 LL ( 



x,Q 2 



S2 W ( 



X 



Q 2 ) + \Y,<fj A ^Q 



(274) 
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where is the Wandzura-Wilczek contribution, expressible in terms of the twist-two struc- 
ture function g\. The transverse quark distribution Aq^(x,Q 2 ) = Aqr(x,Q 2 ) — AqT(—x,Q 2 ) 
is defined as 

t! 

dv (p, s\ [(a + v) S + (a, v, —a) + (a — v)S~(a, v, —a)] \p, s) = 

= ~is±J dxe 2ixap+ Aq T (x,fi 2 ). (275) 

Here \p, s) is the polarized nucleon state with momentum p^ and polarization vector s^. In 
the leading logarithmic approximation, the gluon contribution only appears due to the QCD 
evolution, in particular 

Q 2 ^9t(x; Q 2 ) = g £ ^-pT(x/y)Ag T (y; Q 2 ) . (276) 
The gluon transverse spin density Ag? is defined in terms of the three-gluon operators |2| 
dv (p, s\ (a + v)0(v, a, — a) + 2a 0(a, v, — a) + (a — v)0(a, — a, v) \p,s) = 

—a 

= 2s ± p + Jdxe 2ixap+ xAg T (x,n 2 ) (277) 
with x having the meaning of the Bjorken variable. The evolution kernels are given by 
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(278) 



Here the first two expressions are accurate up to corrections of order 0(1 — x) for x — > 1 and 
the third one has the accuracy 0((1 — x) 3 ). Since quark and gluon distributions in the nucleon 
strongly decrease for x — > 1, the contribution of large x to the splitting functions is the most 
important one numerically 

4.5 The Regge Limit of QCD 

In this section we shall describe another example in which conformal symmetry plays a crucial 
role — finding the asymptotic behavior of QCD scattering amplitude in the high-energy (Regge) 
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limit. In the BFKL approach [2U El ESI partial waves of the scattering amplitude satisfy a 
Bethe-Salpeter-like equation which involves the effective Hamiltonian 7"^bfkl acting on two- 
dimensional plane of transverse coordinates. It turns out that this Hamiltonian has many 
features in common with the evolution kernels considered in the previous Section. It is invariant 
under conformal transformation on the two-dimensional (Euclidean) transverse plane and is also 
related to integrable Heisenberg magnets. In distinction with the previous case, the conformal 
transformations act on the plane (and not on the line) and the corresponding generators form 
the SL(2, C) algebra (and not the SL(2, M.) algebra). The origin of the SL(2, C) symmetry and 
integrability properties of the Hamiltonian 7Ybfkl remains unclear. However, they allow us to 
solve the Schrodinger equation for this Hamiltonian and obtain the high-energy asymptotics of 
the scattering amplitudes. 

Consider the elastic scattering of two colorless hadronic states in the Regge limit, when 
their center-of-mass energy s is much larger than their masses and momentum transferred in 
the t— channel, s ^> —t, m 2 . For the sake of simplicity and in order for perturbative QCD to 
be applicable, one may choose these states to be heavy onia, i.e. the bound states of heavy 
quark and heavy antiquark with small transverse size ~ 1/m <C 1/Aqcd- To the lowest order 
in the coupling constant the elastic scattering of two onia with the momenta p± ~ y/s] 2n^ and 
P2 — \fsj2n p ' occurs through the exchange by two gluons with momenta k^ and q^ — k^ . Using 
the Sudakov parametrization k^ = otip^ + 0:2^2 + k± the four-dimensional integration over the 
gluon momenta becomes d 4 k = (s/2)daida2d 2 k±. It can be shown that at high energies the 
dominant contribution to the scattering amplitude A(s, t) comes from the momentum region 
such that k 2 pa —k\ and ik — q) 2 pa —{k — q)\ and the integrations over the longitudinal 
momentum fractions a\ and 02 get factorised, see jHSj- As a consequence, the two-gluon 
scattering amplitude can be written in the so-called impact factor representation 

/d 2 k 1 
(27T) 2 k 2 (k - q) 2 $l(/C±; q± ~ k± ^ k± > q± ~ k±) ' (279) 

where the functions $1^ (k±, q± — k±) do not depend on the energy s and describe the coupling 
of two gluons to the onia states with momenta p\ and P2, respectively. Gauge invariance 
implies that these functions vanish if momenta of one of the gluons vanishes $1,2(^1, Q± — 
k±) — >• for k± — > or q± — k± — * 0. The same condition ensures that the amplitude 
(I279J) is infrared finite. The two-gluon exchange amplitude (j279|) is of order A(s,t) ~ s 1 
which corresponds to the constant cross section. Going over to higher orders in perturbation 
theory, one finds that the scattering amplitude receives corrections enhanced by logarithms 
of the energy ~ a 2 (a s ln s) m [l + 0(a s )]. The resummation of such corrections in the leading 
logarithmic approximation (LLA) in energy gives rise to the famous result jMJ E] 

A(s,t)~ iS 1+ ^ ln2 , (280) 

known as the BFKL pomeron. 

The resummation of terms ~ (a s lns) m is done by expanding the scattering amplitude over 
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the partial waves with the complex angular momentum 

/■5+ioo 7 _ 

A(s,t)=tsa 2 s — s"A{u,t), (281) 

where the integration contour goes to the right of all singularities of A(u, t) on the complex 
uo— plane. In this way the expansion over (a s \ns) m is traded for the expansion over (a s /u) m 
with the coefficients that depend on t, and the asymptotical behavior in Eq. (|28U|) arises because 
the resummed expression for A(a s /uj,t) develops a nontrivial singularity for positive u in the 
sum of all orders in a s . 

To leading logarithmic accuracy the partial wave A(a s /u>,t) is given by 

*"' t) = /ffe/ ^ $ i(^'^-^) T -(^'^'^)^(^,^-^), (282) 

where T u} (k±, k' ± ; q±) is the partial wave amplitude for the scattering of two off-shell gluons 
with the virtualities k\, (q± — k±) 2 , k' ± 2 , (q± — k±)' 2 and vacuum quantum numbers in the 
t— channel. Comparing Eq. ()282|) with ()279l) one finds the function T w which, in leading order, 
is given by the product of two gluon propagators: 

T u (k±, k' ± ; q ± ) = - [T^(k ± , k' ± ; q ± ) + 0(a s /u)} , T®{k ± , k' ± ; q ± ) = " ^ ,(283) 

A nontrivial analysis which cannot repeat here, see |il4"| IT2*| IHK] . shows that the resummed 
partial wave T u satisfies a Bethe-Salpeter like equation 

uT u {k ± , k' ± ; q ± ) = T®(k ± , k' ± ; q ± ) + ^ [^bfkl ® T u ] (k ± , k' ± ; q ± ) , (284) 

7T 

where the integral operator Hbfkl describes the interaction between gluons which are propa- 
gating in the t— channel. It acts only on their two-dimensional transverse momenta: 

["Hbfkl ® T u ] (k ± , k' ± ;q±) = 

= h I w=wS [k2±{q " 1)1 + {q ~ k)2±l2± k)2±q2±] Tuj{l±j k ' ±] q±) 

k\ (q- kf_ 



A + (k-l)i ' (g-Oi + ^-Oij^^-^^'^^}- (285) 

Solving Eq. ()284|) by iterations corresponds to the perturbation theory in (a s /u). Instead, in 
order to identify the singularities of the partial wave in the complex uj— plane one has to solve 
Eq. ()284j) exactly. To this end it is convenient to go over from the momentum to the impact 
parameter representation: 

$ a {k, q-k) = J d 2 b lC l 2 b 2 e -*( fc6 0-*((9-fc)6 a )$ a ( 6lj ^ ; (286) 
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where it is implied that all vectors belong to two-dimensional plane of transverse coordinates 
(impact parameters) and we use the same notation for functions in the momentum and the 
coordinate space. 

The rewriting of (|282|) in the impact parameter space reads 



A{u,t) = I 'd 2 b e^ I \{ d 2 b k dX^ih -b ,b 2 -b )TUb 1 ,b 2 ;b' 1 ,b' 2 ))^ 2 (b' 1 ,b' 2 ) 

J fc=l,2 

= J d 2 6oe i{9f,o) (<f(6o)|Tj$(0)), (287) 

where in the second line the notation was introduced for the operator whose kernel is given 
by T UJ (b 1 ,b 2 ;b' 1 ,b' 2 ), defined by flUHJ). As the result, the Bethe-Salpeter equation (gHS) 

takes a simple operator form 

uoT u = T£°) + — Hbfkl T w . (288) 

7T 



Solution of this equation can formally be written as 

T w = L- — Hbfkl) TL 0) , (289) 



TT 



where from one concludes that the singularities of in the cu-plane are determined by the 
eigenvalues of the BFKL operator 

[Hbfkl ■ V>a] (K h) = E a ^ a {b u b 2 ) , (290) 

with a numerating the solutions. The high-energy behavior of the scattering amplitude is 
governed by the right-most singularity of T w , which corresponds to the maximal eigenvalue 
m&x a E a . Eq. ()290|) has the form of a Schrodinger equation for the system of two interacting 
particles on the two-dimensional plane. Such particles can be identified as reggeized gluons and 
the eigenstates ip a (bi, b 2 ) have the meaning of the wave function of the color-singlet compound 
states built from two reggeized gluons. 

The BFKL operator Hbfkl has a number of remarkable properties which allow one to solve 
the Schrodinger equation ()290|) exactly [121 HS3 EH] . First of all, Hbfkl splits into the sum of 
two operators acting on the holomorphic and the antiholomorphic coordinates bi 2 — ► (-21,2, ^1,2); 
Eq. (J23)) . on the transverse plane: 

Hbfkl = H 2 + H 2 , H 2 = d' 1 \n(z 12 ) d zi + \n{z 12 ) d Z2 + \n(d zl d Z2 ) - 2*(1) , (291) 

where z\ 2 = z\ — z 2 and 7i 2 is given by the similar expression in the z— sector. Consider the 
generators of the holomorphic SL(2, C) transformations ()25|) corresponding to the field with 
h = h = 6 

Lk~ = ~dz k , L k = Zf.d Zk , L fci+ = z\d Zk , (292) 

6 We would like to stress that these quantum numbers correspond to a scalar field in a two-dimensional field 
theory and are different from transverse components of a physical gluon field. Because of this the SX(2,C) 
symmetry of Hbfkl does not follow immediately from the conformal symmetry of the QCD Lagrangian. 
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and the corresponding antiholomorphic generators L^,q and Lk,+ given by similar expres- 
sions with Zk replaced by z^, with k = 1,2 enumerating particles. By inspection one finds that 
Hbfkl commutes with all two-particle generators 

[HbfKL, Lx >a + L 2>a ] = [HIbfKL, Lx,a + L 2>a ] = . (293) 

with a = +, — , and, therefore, is invariant under the SL(2, C) transformations of the trans- 
verse plane. This implies that Hbfkl only depends on the two-particle Casimir operators of 
the SL(2, C) group 

h\ 2 = -( Zl - z 2 ) 2 d Zl d Z2 , L 2 12 = -(z 1 -z 2 ) 2 d- Zl d- Z2 , (294) 

and leads to TC 2 = 7i 2 (L\ 2 ) and 7i 2 = 7i 2 (Lf 2 ). As a consequence, solutions of the Schrodinger 
equation ()290|) have to be eigenstates of the Casimir operators 

L\ 2 ip n ,v = h(h - l)Vv ' ^LVv = Hh ~ ■ ( 295 ) 
Here the pair of complex conformal spins 

1 -)- Tl — 1 — 71 

h = — h iv and h = — h iv , (296) 

with a nonnegative integer n and real v specify the irreducible (principal series) representation 
of the SL{2, C) group to which ip ntU belongs to. The solutions to Eqs. ()295|) read 

/ _ \ (l+n)/2+iv , \ (l-n)/2+iu 

^AKh)={^-) (^A , (297) 

where Zjk = Zj — Zk and b = (z , z ) is the collective coordinate, reflecting the invariance of 
Hbfkl under translations. We recall that the solution ()297j) defines the wave function of the 
color-singlet compound state, which is built from two reggeized gluons with the coordinates 
b\ = (zi,zi) and b 2 = (z 2 ,z 2 ). The integer n fixes the two-dimensional Lorentz spin of the 
state, real valued v gives the scaling dimension I = 1 + 2iu, and the two-dimensional vector 6 
sets up the center-of-mass coordinate of the state. 

We conclude that, similar to the two-quark evolution equations discussed in Sect. 4.1, the 
solution to the Schrodinger equation for the compound states of reggeized gluons is uniquely 
fixed by conformal invariance. We would like to stress the following important differences 
between the two cases. The conformal SL(2, R) symmetry of the evolution equations is a 
consequence of the symmetry of the QCD Lagrangian, but the origin of the SL(2, C) symmetry 
of the BFKL operator remains unclear. In the former case, the conformal group acts along 
the line defined by the light-cone direction and conformal spins of operators take positive 
(half)integer values. In the latter case the conformal group acts independently along two 
"orthogonal" z— and z— directions on the two-dimensional Euclidean plane and conformal spins 
of gluonic states take complex values specified in Eq. 
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To find the eigenvalues in Eq. (|29U|) . one substitutes the wave function (|297|) into the 
Schrodinger equation (|29U|) and uses the explicit form of the BFKL kernel (|285|) leading 

to miming 

E n>v = 2^(1) - * ( + iv) - ( - iu) . (298) 



Its maximal value maxE ntU = 4 In 2 corresponds to n = v = 0, or equivalently h = h = 1/2. It 
defines the position of the right-most singularity of the partial wave amplitude and is translated 
into the BFKL behavior of the scattering amplitude in the leading logarithmic approximation, 
Eq. ()28()j) . Due to accumulation of the energy levels E n>u around n = v = 0, the corresponding 
Regge singularity is not a pole but a square-root cut so that ()28()j) is modified by additional 
(a s In s)" 1 / 2 — factor. 

Using ()298|) one can reconstruct the operator form of the BFKL kernel Hbfkl on the rep- 
resentation space of the principal series of the SL(2, C) group 

Hbfkl = \ [H(J 12 ) + H{J 12 )] , H(j) = 2*1/(1) - (j) - (1 - j) , (299) 

where, as before, the two-particle spins are defined as L\ 2 = Jw^Ju — 1) and L\ 2 = Ju^Ju ~ !)• 
Notice that we already encountered similar Hamiltonians in Sect. 4.1 and later found in Sect. 4.3 
that they give rise to complete integrability for the three-particle evolution equations for the 
helicity— 3/2 baryon distribution amplitudes. It turns out that the BFKL kernel ()299|) has the 
same hidden symmetry. In order to see this, we have to go over to the states containing more 
than two particles. In the Regge limit, this amounts to considering color-singlet compound 
states built from three and more reggeized gluons, or briefly N— reggeon states. In this termi- 
nology, the BFKL pomeron is a particular example of the two-reggeon state. In spite of the 
fact that the contribution of iV-reggeon compound states is suppressed by the factor a^~ 2 as 
compared with the leading term, they have to be taken into account at very high energy and 
are expected [HZ] to tame a power rise with energy of the BFKL solution (J28U)) . which is in 
conflict with unitarity constraints. 

The contribution of the iV— reggeon state to the scattering amplitude is given by the ex- 
pression similar to ()287|) in terms of the impact factors $1,2(^1, • • • > &iv) that take into account 
the coupling of iV gluons to the onia states and the kernel T^bi, ... , b N ; b[, . . . , b' N ) describ- 
ing propagation of iV interacting reggeons in the t— channel. The corresponding operator T w 
acts on the two-dimensional coordinates of iV reggeized gluons and satisfies the Bethe-Salpeter 
equation (|289j) . in which corresponds to the product of iV gluon propagators and Hbfkl is 



replaced by the iV— reggeon effective Hamiltonian 

Mn = -^W E m[H(J jk )+H(J jk )], (300) 



l<j<k<N 



where the sum goes over all pairs (j, k) of reggeized gluons. Each term in the sum has the 
color factor given by the product of the color charge of interacting gluons. It is multiplied by 
the pair- wise Hamiltonian ()299|) . depending on the two-particle conformal spins Jjk(Jjk — 1) = 
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—z 2 k d Zj d Zk and Jjk(Jjk~ 1) = ~ z ]k^s j 9z k - The iV-reggeon states satisfy the Bartels-Kwiecinski- 
Praszalowicz equation 

N 

[M N -i> a ](b 1 ,...,b N ) = E a iJj a (b l ,...,b N ), ^2qiJ> a (b 1 ,...,b N ) = 0, (301) 

i=i 

which is a generalization of the Schrodinger equation ()290|) to A^-particle system with pair-wise 
interaction j^Z] . Here the second relation ensures that the total color charge of the state equals 
zero. 

At N = 2 and N = 3 the color factor in (J301j) can be reduced to a number t^t^ = —N c 
and tjtfjl = —N c /2, respectively. For N > 4 several color structures exist and finding the 
general solution to the iV-particle Schrodinger equation (j301|) becomes an extremely difficult 
task. A considerable simplification occurs in the large- N c limit. In this case the relevant 
Feynman diagrams have a topology of cylinder so that the interaction takes place between 
nearest neighbors only, leading to t a -t a k = —8j t k+iN c /2. Thus, in the multi-color limit the 
Hamiltonian (|300|) takes the form 

Mn n c ^oo 1 [H(J kM1 ) + H(J ktk+1 )] , (302) 

l<fc<iV 

where Jn+i,n — Ji,n an d similar for J N+ i <N jHUllIB]- Notice that ()302|) is exact for = 2 and 
N = 3. 

Since the Hamiltonian (J302)) depends on the two-particle spins, it commutes with the total 
conformal SL(2, C) spin of the system L a = J2k=i Lk,a and L a = J2k=i ^k,a- This allows one 
to impose the conformal constraints on its eigenstates 

L 2 ip n ,u = h{h - l)tp n>v , L 2 tp n>u = h(h - > ( 303 ) 

where the conformal spins h and h are defined in Eq. (|296|) . For > 3 these conditions do 
not fix the eigenstates uniquely. Similar to the situation with the evolution equations for three 
quark operators, additional constraints follow from the analysis of symmetry properties of the 
Hamiltonian. 

It turns out that the Hamiltonian ()302|) possesses a set of integrals of motion and 
(with k = 2, . . . , JV) [HI d 

Qk = i k z ™ ■ ■ ■ z h-i3k z 3kh d z h ■ ■ ■ d z Jk _, £>z 3k , (304) 

l<Jl<-<ifc<AT 

Qfc are given by the same expression with z replaced by z. It is easy to verify that Q2 = L 2 
and Q2 — L 2 . The holomorphic charges satisfy the following commutation relations 

[Qk, Qn] = [Qk, L a \ = [Q fc , M N ] = (305) 

and similar relations hold in the antiholomorphic sector. The existence of the integrals of 
motion implies that the Schrodinger equation ()30H) is completely integrable in the multi-color 
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limit. Hence the Hamiltonian depends only on the total set of conserved charges Q2, • • • , Qn 
and their counterparts in the antiholomorphic sector. This allows one to replace the eigenvalue 
problem (j301|) by a simpler one of finding simultaneous eigenstates for the set of conserved 
charges 

Qki>n,v, q (bi, • • • ,M = qk ip n ,u,q(bi, ■ ■ ■ ,M with k — 2, . . . , N , (306) 

supplemented by the dependence of the energy E(q 2 , . . . , g^v) on the integrals of motion. The 
quantization conditions for the charges qk follow from the requirement that ip n ,u,q{bi) is a single- 
valued function on the two-dimensional plane, normalizable with respect to the SL(2, C) in- 
variant scalar product. 

We have seen that the evolution kernel for helicity— 3/2 baryon operator coincides with 
the Hamiltonian of the Heisenberg magnet with the spin operators being the generators of 
the SL(2, R) group. The iV-reggeon Hamiltonian (|3U2|) has the similar hidden symmetry in 
the multi-color limit. Remarkably enough, for arbitrary N this Hamiltonian coincides with 
the Hamiltonian of the Heisenberg magnet with N sites and the spin operators belonging 
to the principal series representation of the SX(2,C) group [Tl|. As follows from Eq. (|292j) . 
L\ = so that the single particle spin equals zero. These features allow one to apply powerful 
methods of integrable models to solve to the spectral problem (J3U2j) . The current situation is 
the following. For N = 3 the solution was found in Refs. [9*81 I99| . For N = 4 the spectrum was 
recently calculated by two groups Refs. [100, HOlj and |102j with somewhat different results. 
The spectrum of 5 < iV < 8 states and the generalization to arbitrary iV were considered in 
Ref. |101j . We refer an interested reader to these papers for further details. 



5 Conformal Symmetry Beyond the Leading Logarithms 

Beyond the leading logarithmic approximation conformal symmetry is certainly broken by the 
trace anomaly of the energy-momentum tensor. In the present section we analyse this breaking 
in detail and establish the following schematic structure of the perturbative series for a generic 
quantity Q 

Q = Q con + ^-AQ , where AQ = power series in a s . (307) 

Here Q con is the result in the formal conformal limit, obtained by setting the /3-function to 
zero by hand. It has full symmetry of a conformally invariant theory. The extra term AQ can 
be perturbatively evaluated in a power series with respect to a s and vanishes to leading order 
(LO). Note that the evaluation of the leading 0(a s ) contribution to AQ requires only little 
effort, since it is sufficient to calculate the Nf proportional terms in (3q = lliV c /3 — 2Nf/3 via 
quark bubble insertions. 

The possibility of the separation of conformally symmetric and /3-proportional contributions 
is by no means trivial. In particular, since the anomalous dimensions are scheme- dependent 
beyond LO, Q con appears to be scheme-dependent as well and one might expect that the 
representation in ()307|) can only be valid in a certain scheme, if at all. We will indeed find that 
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perturbative corrections in general violate conformal Ward identities and destroy the symmetry. 
The most important result is, however, that the n-loop anomalous contribution to the dilatation 
and the n — 1-loop special conformal anomaly are related to each other and they can both be 
removed simultaneously by a finite renormalization. Neglecting the /3-function, the conformal 
covariance can, therefore, be restored in a special, conformal scheme. 

This general result can be used in two different ways. In the conformal scheme, mixing 
matrices involving operators with total derivatives (alias full ER-BL kernels) are completely 
fixed (to all orders in a s ) by the diagonal entries. This implies that, e.g. two-loop off-forward 
anomalous dimensions (or two-loop ER-BL kernel) in an arbitrary scheme, say the modified 
MS (MS), can be restored from their diagonal (forward) parts, making a finite renormalization 
from the conformal to the MS scheme. This transition requires a calculation of the one-loop 
contribution to the special conformal anomaly and is considerably simpler than the full two- 
loop calculation. Another option is to abandon the MS-like schemes completely and do the 
whole analysis in the conformal scheme. We will give examples for both strategies. 

The presentation is organized as follows. In Sect. 15.11 we consider the so-called Crewther 
relation in which case the above-mentioned difficulties are absent. This relation has been tested 
to 0{a^) by explicit calculations and provides a beautiful illustration of the power of the con- 
formal symmetry. In Sect. 15.21 we consider conformal Ward identities in an interacting theory 
and show that the conformal covariance can be restored (up to corrections in (3) in a special 
renormalization scheme — we call it conformal subtraction (CS) scheme — that differs from a 
MS-like scheme by a finite renormalization. In particular the conformal operator product ex- 
pansion, derived in Sect. ESI holds true in the CS scheme beyond the leading order. In Sect. 15 .HI 
the leading order special conformal anomaly is calculated and we demonstrate that using this 
anomaly and imposing the conformal constraints one can derive the NLO expressions for the 
matrices of anomalous dimensions of all twist-two operators in a MS scheme and reconstruct 
the corresponding evolution kernels. Finally, in Sect. 15.41 we consider the application of the 
conformal operator product expansion to light-cone dominated two-photon processes. 



5.1 The Crewther Relation 



The Crewther relation [^Uj (see also |103j ) predates QCD and provides the relation between 
three quantities: the Bjorken sum rule for the structure function of polarized deep inelastic 
scattering 



dx(<fi-g?) (x,Q 2 



9a 



9v 



C Bi (a s (Q 2 )) + 0(l/Q 2 



the isovector part of the ratio of the total cross-sections 

hadrons) 



R(s) 



a(e + e 



a e T e" 



(30£ 



(309) 



at large center-of-mass energy s and the axial anomaly constant, which governs the pion decay 



7T 



77. It was later extended to QCD in Ref. 



In its present form, the so-called 
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Figure 7: The VVA correlator ((HEj) at 
tree level (a). In (b) we display one con- 
tribution out of six to the r.h.s. of the 
CWIs ()315jh where the gray blob de- 
notes the renormalized trace anomaly. 



Ja(x) ( a ) Mv) 



(b) 



generalized Crewther relation reads 



C B j{a s )D(a s ) = 1 + 



[power series in a s ] , 



(310) 



9 



where D = D^s is the flavor non-singlet part of the Adler D function |lU5j . It is related to the 
R— ratio by the dispersion relation 



The Crewther relation ()310|) relates the perturbative expansion of two coefficient functions and 
it is valid up to corrections 0(1/Q 2 ). In comparison with its original form [20 j5 Eq. (|31()j) 
contains an additional term proportional to the ft— function 7 . The Crewther relation can be 
extended to the flavor singlet sector in which case it relates the coefficient function entering 
the Gross-Llewellyn Smith sum rule, Cqls^s); an d the singlet part of the Adler D— function. 
The validity of the Crewther relations in both the flavor singlet and non-singlet sector has been 
verified using the available QCD results for C— and D— functions up to the order C?(af), 



where J a is the electromagnetic current and is the isovector axial- vector current. Neglecting 
quark masses the both currents are conserved in QCD and have zero anomalous dimensions so 
that their scaling dimensions are not altered by the interaction, i = £ can = 3. 

If the conformal symmetry in QCD were exact, the three-point correlation function (|312j) 
would be uniquely fixed up to normalization |lll)| (see Sect. 2.4.2). Moreover, since the scaling 
dimensions of the currents do not depend on a s , C a ^ would be given by the lowest-order triangle 
diagram in Fig. Ufa), up to an overall normalization. Notice, however, that the normalization 

7 In pre-QCD times it was believed that the conformal symmetry is only softly broken, i.e., the trace of the 
energy momentum tensor would provide only power corrections. 




(311) 




C aM (x -y,z-y; a s ) = (0\TJ a (x)Jp(y)J* 



(*)|0> 



(312) 
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of Eq. (|312|) is protected by the U(l) axial anomaly condition 

-C aPl {x -y,z-y;a 8 ) = e a ^ v —— — 5 {A) (x - z)5 {4) (y - z) . (313) 



According to the Adler-Bardeen theorem the r.h.s. of this relation does not receive ra- 

diative corrections and, therefore, in the conformal limit, C aj g 7 coincides with the lowest-order 
result that we denote as A Q/ 3 7 (see Fig. 0(a)) 

C af 3y(x -y,z-y) = A a/37 (x -y,z-y). (314) 

In QCD, Eq. ()314j) is modified because of the non- vanishing (3— function. To understand the 
modifications, note that the currents J a and are transformed under conformal transforma- 
tions as primary fields of spin— 1 and scaling dimension I = 3. Hence C Q/ g 7 obeys the conformal 
Ward identities (fTToT) and (TT251) : 

[V x + V y + V z }C a ^(x-y,z-y) = ^l(0\TJ a (x)J p (y)J^(z) J fa [G a ^G aflu ](w)\0) ,(315) 

[^ + ^ + }C^C a ^(x-y,z-y) = ^(0\TJ a (x) My) J*(z) J dfa2vf[(^G»»](v))\0) . 

Here we introduced a shorthand notation for the differential operators that enter in the l.h.s. 
of ffTTKjl and (fraSjl : 

V x = £ + x-d x and K% = 2x» x ■ d x - x 2 d£ + 2£ x^ - 2x u W u (316) 

with £ = 3 and being the operator of spin rotations for a vector field, Eq. (0). [G a ^ u G aflv } 
denotes the operator renormalized in a MS-like scheme. 

For /3(g) = one can verify that the solution to the conformal Ward identities (J315j) . 
supplemented by the additional condition ()313|) is given by Eq. (]314j) . For (3(g) ^ the 
anomalous terms entering the r.h.s. of Eqs. (|315]) are described by diagrams such as in Fig. 0b). 
They modify the general solution to the system of two inhomogeneous differential equations 
(I315|) by terms proportional to the (3— function 

C a/ 3 7 (x -y,z-y; a,) = A a/3 ^(x -y,z-y) + ^^-A a ^(x -y,z-y; a s ) , (317) 

where A Q ^ 7 = 0(a s ). Since this addition has still to satisfy the anomaly condition (|313J) . 
A a/ 3 7 must have a vanishing derivative with respect to z 7 and also respect the vector current 
conservation. 

On the other hand, the correlation function 1)317]) can be calculated at short distances using 
the operator product expansion. The idea [2H1 11031 1104] is to consider the non-symmetric limit 
with | a; — y\ <C \z — y\, so that one has to expand the product of two electromagnetic currents 
J a (x)Jp(y) as the first step. Choosing y — for simplicity we obtain 

TJ a (x)M0) ~ C* (*)/ + C^(x)J^(0) + ■■■, (318) 
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where the coefficient functions C^{x) and (x) define the perturbative expansion of Cbj(o! s ) 
and D(a s ), Eqs. (|3U8|) and (|311|) . respectively. Inserting this expansion into the correlation 
function (j312j) we find that the contribution of the unit operator in (j318|) vanishes. Making use 
of the operator product expansion for z we obtain 

TJ*(0)J*(z)^C*(z)I + --- , (319) 

so that finally 

C afh (x, z) ~ Cj," (x)C*(z) , x,z^0, x^z. (320) 

Notice that C*L(z) coincides with the isovector part of the correlator (0|TV a (x) Jp(0)\0) = 
C^p(z) and, therefore, it contributes to the flavor non-singlet part of the D— function, Eq. (j310|) . 
Equating the conformal prediction (j317|) and the OPE result ()320|) one obtains nontrivial con- 
straints on the perturbative expansions of the coefficient functions Cj^(z) and C a J i {x). Going 
over to the momentum space one finds |l()4j that Eq. (j317j) and Eq. ()32()|) give rise to the expres- 
sions entering the r.h.s. and the l.h.s. of the Crewther relation ()31()j) . respectively. Repeating 
the similar analysis in the opposite limit z<i with x, z — > one can derive the Crewther 
relation in the flavor singlet sector. Since A a/?7 = 0(a s ), the perturbative series on the r.h.s. 
of Eq. ()310|) starts at order ot 2 s . 

It is worth mentioning that within the so-called Brodsky-Lepage-Mackenzie scale setting 
prescription |112j in which the terms in (3 are absorbed into the redefinition of the scale of 
the running coupling, the generalized Crewther relation (j310|) takes the same form as in the 
conformal invariant theory |113j : 

C(Q)D(Q*) = [1 - &c(Q)) [l + a D (Q*)] = 1 , (321) 

where cxq and &£> are the corresponding effective charges satisfying a geometrical commensurate 
single scale relation. Eq. (j321|) has been verified using the known 0(a^) results, but a general 
proof has not been given so far. 



5.2 Conformally Covariant Renormalization Scheme 

The case of the Crewther relation is special in that all relevant operators are protected by exact 
symmetries and are not renormalized. In the general situation the transformation properties of 
conformal operators can easily be spoiled by renormalization and this is indeed what happens 
in the usual MS scheme. For the sake of simplicity we consider the conformal tower of twist-two 
quark- ant iquark operators in the flavor non-singlet sector 

Q n , = (id + y- n Qy(0), (322) 

where Q^ ,1 (0) is defined in Eq. ()64|) and / > n is the total number of derivatives. For a 
moment let us take the formal conformal limit and set the (3— function to zero. Conformal 
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Ward identities (CWIs) for the Green function made of the renormalized operator [Q n j] (here 
and below [. . .] stands for the renormalized quantity in a MS-like scheme) and arbitrary number 
N of fundamental fields ([O n j]Af/v) = (0| [O n ;]$(xi) . . . $(xat)|0) have the following general form 
(cf. Eqs. (IB~28j) and (JESI): 

N n 



Y, V i{Pnl]X N ) = -J2i £ T^nm + % m ]([O ml ]X N ), (323) 
i=l m=0 
N n 

^JCr ([Q nl ]X N ) = iJ2[ a (n,l)S nm + r nm (l)}(Pmi-i}X N ), (324) 



i=l m=0 



where the gauge-fixing and ghost terms are not displayed, since they do not contribute to 
physical quantities. Here T>i and = n^fC^ denote the differential operators (j316j) acting on 
the primary fields in the monomial Xn, with the scaling dimensions of fields = £| an + 7$. 
The terms ~ 5 nm on the r.h.s. of Eqs. ()323|) . (j324j) originate from the same differential operators 



acting on O n ; which has the Lorentz spin l+l and the canonical dimension £[ an = 2£^ n +l = 3+1. 
The coefficient a(n,l) is equal to 

a(n,l) = -2(l-n)(l-n + 2j n - 1) = -2(1 - n)(n + I + 3) = 2(n - l)(£ c ™ + I) , (325) 

where in the first expression j n = (£^ n + n+ l)/2 is the (canonical) conformal spin, cf. Eq. (fTU|) . 

The remaining entries are the anomalous- dimension matrix 7 = {^ nm } and the so-called 
special-conformal anomaly matrix 7 C (Z) = {7^ m }, cf. Eqs. f|B.30jl and (jB.31j) . They both are 
induced by the trace anomaly and have the standard perturbative expansion 

7(a.) = ^7 (0) + t^7 (1) + 0(al) and f(a s ) = ^T (0) + 0(a 2 s ) . (326) 



We recall that the dilatation Ward identity (|323j) is nothing else but the Callan-Symanzik 
renormalization group equation and what it tells us is that conformal operators in general can 
mix with the operators with the same dimension and lower conformal spin. Due to Poincare 
invariance, the both matrices 7 and j c are in general triangular and, in addition, 7 is indepen- 
dent on the Lorentz spin l + l. At the same time, since the special conformal transformations 
do not commute with the Poincare transformations ([K a ,M^ 7 ] 7^ and [K a , P/?] 7^ 0), both 
a(n, I) and 7 C (0 depend on /. The conformal invariance is broken if the conformal anomaly 
matrices 7 and J C {1) are not diagonal. For instance, if the special conformal anomaly matrix 
7 C (Z) has nonzero off-diagonal terms, then the operator [O nn ] is not annihilated by /C_ and its 
variation is expressed in terms of the operators [O mn _i] with m < n — 1. 

The crucial point is that the two conformal anomalies, 7 and 7 C (Z)> are n °t independent. 
The relation between the two matrices is imposed by the collinear conformal algebra. Using 

N N 

^Pi+iPra] X N ) = i([Om+i] X N ) , ^A^,_ + ([O nJ ] X N ) = -{I + \){[Q nl ] X N ) , (327) 
i=i i=i 
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and applying the identities 



[D, £_] = £_, [JC-,V+] = -2{V + M-+) (328) 
to ([O n /]Ajv) one obtains the desired constraints. They consist of the commutator relation 

[a(l) +7 c (/;a s ),7(a s )] = where a(l) = {a(n, l)S nm } (329) 

and hence 

2(n - m)(n + m + 3)7 nm (a s ) = [7(a s ),7 c (/; a s )] nm , n > m , (330) 
between the anomalous-dimension and special-conformal anomaly matrices, and the relation 

f{l + 1; a a ) - f{k a.) = -2j{a s ) , (331) 



which determines the spin dependence of the special conformal anomaly. Eq. (|33H) ensures 
that the I dependence in the commutator relation (j330p is spurious and cancels out. 

The constraint ()329|) immediately tells us that the anomalous dimension matrix is diagonal 
in LO, since its perturbative expansion gives 

[a(l), 7 {0) ] nm = 2(n-m)(n + m + 3) 7 ffi = , n > m . (332) 

Beyond this approximation one realizes from Eq. (J33U)) that the special conformal anomaly 
induces the off-diagonal entries of the anomalous dimension matrix, starting 0(a 2 s ). 

It is instructive to take a closer look at what is happening to order 0(a s ). To this ac- 
curacy 7 im = jn 5 nm and the dilatation covariance is preserved. In contrast, the covariance 
with respect to the special conformal transformation is lost, i.e., the special conformal anomaly 
matrix 7 C (°) is non-diagonal in a general renormalization scheme. The reason for this is that in 
distinction with one- loop anomalous dimension 7 i°- > , which is induced by UV-divergent parts 
of Feynman diagrams, the special conformal anomaly receives contribution also from UV finite 
parts. Notice, however, that the off-diagonal entries of 7 c< -°' ) (/) can be removed and, as a conse- 
quence, the special conformal covariance can be restored by performing a finite renormalization 

n-2 c (0) 

OS? = W + ^Et^W + «- ( 333 ) 

2tt z — ' aim, n) 

m=0 s ' 

This transformation defines the conformal subtraction (CS) scheme. Indeed, inserting Eq. (J333j) 
in (J324j) and taking into account the identity a(n, I) = a(n, m) + a(m, I) we find that the 
operators O^f do not mix under infinitesimal special conformal transformations and satisfy the 
Ward identity 

N 

]T/Cr (OSSX N ) = i [a(n,l) + ^tSHO] <OS?-i**> + (334) 
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On the other hand, the finite renormalization (|333|) does not affect the dilatation Ward identity 
(|323|) to order ot s : 



N 

a. 



£ A (O c jX N ) = - [tr + ^7i 0) J (O c jX N ) + 0{a% (335) 
i=i 

As we will show below in Eqs. fl356j) and ()359|) . the diagonal matrix elements ^nn\l) are related 
to the anomalous dimensions 7„°' 1 as 7nn (0 = 2(n — l)^^ ■ Therefore, in both the dilatation and 
the special conformal Ward identities, Eqs. (|335|) and ()334|) . the conformal symmetry breaking 
terms, containing one-loop anomalous dimension, can be absorbed through the redefinition of 
the scaling dimension 

C n 4(«,)=C n + 7n(a s ), (336) 
a(n,l) a(n,l) +2(n-l)j n (a s ) = 2(n - l)(£ n (a s ) + I) . 

Thus, the conformal properties of the renormalized operators ()333j) are restored at leading 
order. As in a free theory, the operators ()333)) form the irreducible conformal tower ()4())1 but 
their conformal spin is modified by the anomalous dimension term j n = (£ n (a s ) + n + l)/2. 

Beyond the leading order, the commutator relation ()330|) tells that the special conformal 
anomaly gives rise to off-diagonal entries in the anomalous dimension matrix. Inserting the 
decomposition of the anomalous dimension matrix in the diagonal and off-diagonal parts 

^ = ^+^ND ; ^ ^D = j 7rAm j and ~ND = j 7nm £( n>m )j ; (337) 

into Eq. ()330|) leads to a recurrence relation for 7 ND , which can be solved perturbatively 

7 ND = -^7 D = -£7 D + £ 2 7 D -"--, with gA:=l 6{n ( >m ? W(i),A\ ) (338) 
1 + Q y a(n, m) L J nm J 

by an expansion in powers of 7 c (a s ) JH]- The / dependence appearing here on the r.h.s. is 
spurious. 

The conformally covariant operators O^f are obtained from those in a MS-like scheme by 
the rotation 



n 



E B ™Prm} , (339) 



m=0 



where the matrix B = {B nm } is defined by the requirement that it diagonalizes the anomalous 
dimension matrix 7: 7° = B~ 1 jB. It is given by the following expansion JH| 

B = , 1 — = 1 + £ 7 ND + C ( 7 ND £ 7 ND ) + • ■ ■ with LA := \-0{n > m) A ™ ) (340) 

1 — £7^^* I 7n 7m J 
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Inserting Eq. ([338)1 into Eq. (|34(J|) . one finds after some algebra that the diagonal anomalous 
dimension matrix cancels out, so that the rotation matrix only depends on the special conformal 
anomaly 



B = — ^ = 1 - Jf + J{fjf) , with J A := le(n > m) Anm \ . ( 34 1) 

It remains to be shown that the operator O^f is covariant with respect to the special conformal 
transformations. To this end we apply the transformation ()339|) to the special conformal Ward 
identity (J324jl . Hence O^f appears on the l.h.s. while the r.h.s. reads 

B- 1 [a(l) + 7 C (/)] B = a(l) + B^i^l), B] + B- l f{l)B = {2{n - l)(n + 1 + 3 + j n )5 nm } -(342) 

Here the identity [a(l),B] nm = — {^ c {n)B} nm (see Eq. (A5) of Ref. [21]) has been used. This 
proves that the special conformal anomaly ()342j) in the CS scheme is indeed diagonal and has the 
form ()336j) . The conformal spin of the operators is modified by the anomalous dimensions and 
the operator with the highest conformal weight is annulled by the collinear special conformal 
t r ansf or mat ion . 

We conclude that a CS scheme exists in which covariance of conformal operators holds true 
in the formal conformal limit. Beyond the conformal limit the constraint ([330)1 is modified by a 
(3 proportional term while the condition (j33ip remains unchanged. Besides the commutator re- 
lation ()328)1 and the conformal Ward identities (see (JB.29JI ). the derivation of the full constraint 
requires the renormalization group equation for the trace anomaly. The derivation is tricky and 
all details for the scalar theory in dimensional regularization are given in |18j . The analogous 
calculation for the gauge field theories leads as expected to the absence of any explicit gauge 
dependence in the commutator relation 



a(/)+ 7 c (0 + 2^6(/),7 =0, 6(0 = {b nm (l, 3/2)} . (343) 
The matrix b nm (l, v) appearing here is defined by the expansion of 

C u n (2u - 1) = J2 h nm(h v)C v J2u - 1) (344) 

n=0 



21 - (2u - 1) — 

du 



in terms of Gegenbauer polynomials and reads 

b nm (l, v) = 26(n > m) {(I + m + v)5 nm - [l + (-1)"-™] (m + u)}. (345) 

The solution to the constraint ([343)1 . replacing 7 C (/) by 7 C (0 + 2/3(g)b(l)/g in Eq. ([338)1 . is 
now the off-diagonal part of 7 in the full perturbative theory. Knowing the special conformal 
anomaly at nth order in a s allows us to predict the off-diagonal entries to order n + 1. In NLO 
they read 

(0) _ (o) 119 

tSS (1) = 2 V^ L W ) -^ 6 -) with = ( 346 ) 
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(a) (b) (c) (d) (e) 

Figure 8: The special conformal anomaly at one-loop order in the quark-quark (a, b, c) and 
gluon-quark (d, e) channels. Here the symbol ® and the gray blob denote the operator insertions 
of the conformal operator O n i and J dx2x~A(x), respectively. 

while the evaluation of ^nm remains. 

In the next Section we explain the calculation of the LO special conformal anomaly, give the 
results for both the flavor non-singlet and singlet sectors, and also consider the reconstruction 
of the NLO evolution kernels based on the relation in Eq. f)346j) . 

5.3 Evaluation of Conformal Anomalies and Evolution Kernels 

As the simplest example, we consider the conformal anomalies for flavor non-singlet quark- 
antiquark operators 

[V y + V g ] ([<0yV(2/)V^)} = -(l + 3)([O nl ]tfj(y)^(z)) - i Jd d x(0\[Q nl ]A(xW(y)^(z)) , (347) 

[JC- + JC-] (p nl )il>{y)i>{z)) = ia(nMPm-my)^z))-iJd d x2x- ([O nl ] A(x)1>{y)$(z)), 

(348) 

where (. . .) stands for the vacuum expectation value of the T-product and [O] is the renormal- 
ized operator insertion in a MS-like scheme. For this case it is enough to consider the abelian 
gauge theory, the QCD result can be restored by adding the color factor Cf- It proves to be 
advantageous to choose the Landau gauge £ — > in which the "quark" anomalous dimension 
is zero, 7^ = 0. Within these specifications we have, cf. (|B.14|) 

A(x) = {d-^jG^G^ix)-^^)^ ^(x)-{d-2)^(A tM {x)d 1/ A v (x))+highex orders in a s . 

(349) 

The last, gauge-fixing term in A(x) does not induce UV divergencies and can be omitted. 

For the dilatation Ward identity ()347|) we only have to consider the first term in A(x), 
the insertion of | J d d x {G^ v ) 2 {x) in an internal gluon line, see Fig. |Hfa)-(c). This is simply 
i(g^ v k 2 — k^ky), so that the effective gluon propagator with this insertion is just the usual gluon 
propagator (in Landau gauge) 

Da(k) % (g^k 2 - k^K) D v p{k) = D a/3 (k) with D Q/3 (k) = ^ (g a/3 - k a kp/k 2 ) . (350) 
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Hence we are left with the usual one-loop diagrams and their divergent part is equal to 
— (a g /47re)7n O n j, where "fn^ is the one-loop anomalous dimension of Q n j. In the notations 
of App. B, Zn = In /2. Inserting this result in (|347j) completes the calculation, reproducing 
Eq. (|323|) to the desired accuracy. 

In the case of the special conformal Ward identity ()348|) we have two contributions. The in- 

sertion of i J d x 2x~tp(x)i$ ip(x) is a EOM operator and it can be treated by partial integration 
in the functional integral. For the relevant (divergent) part we have 



d d x2x-{[O nl fy(x)0i>(xW(y)TP(z)) = -Jd d x2x-(( j^i?) +^ x )j^j) VWC*))- 

(351) 

Using the explicit expression for O n ; this can be brought to the form 



8i/j(x) 5ip(x] 

as „,(oA f- a , - a \ln3/2 f 9 x -d y 



x=y=0 



where d x = D x + etc. and the prefactor comes from the renormalization of O n /. This is a local 
operator with the dimension by one unit less than O n i, so it can be expanded in conformal 
operators O m ;_i with m < n. To obtain the coefficients note that the multiplication by x~ 
corresponds in the momentum space to the derivative with respect to the incoming momenta, 
schematically 2i (d p + d q ) [Cn 2 (^^)(p + q) 1 ], and this is the same differential operator that 

appears in the definition of the matrix b(l) in Eq. (J344j) . Consequently, we have the expansion 

/ d " x 2x " (S*<*> + f(x) W>) = 2! ( x + £ b " Mo -"- 1 ■ (352) 



(0) 
m 



where O m /_i on the r.h.s. are not renormalized and possess the UV divergencies — (a s /47re)7 
that have to be subtracted. We find 

/n 
d d x2x-([O nl ]$(x)i3 ^(x)ip(y)ip(z)) = -^E^'^]"-" ([Ond-iMvMz)) (353) 
m=0 

with the 6-matrix given in Eq. ()345|) (for v = 3/2). 

The second part of the special conformal anomaly is induced by the insertion of the operator 
^Jd d x2x~ (G^ u ) 2 (x), see Fig. [SJ (a)-(c). In difference to the case of the dilation anomaly ([35 Up . 
the operator insertion now contains a derivative with respect to the loop momentum. The part 
of the effective gluon propagator that induces a UV divergence can schematically be written as 



2n- dk D aP {k) - D aP {k) 2n- dk ■ (354) 
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In the case of Fig. |S] (a) we can use momentum conservation on the both "quark-gluon-quark" 
vertices to replace these derivatives by those with respect to the external momenta. The result 
for this diagram is thus proportional to the corresponding part of the anomalous dimension 
and the derivatives give rise to the same 6-matrix as above. In the case of Feynman diagrams 
in Figs. Efb) and (c) there is a similar contribution, so that the whole anomalous dimension of 
O n i is restored. However, the derivative with respect to the loop momentum also acts in this 
case directly on the operator insertion O n ; and in this way additional terms arise that will be 
denoted as w nm . More precisely, we define 



n 

^2x-([O ni ]-[Gy(x)^(^W) = ^^{-7 (0) K0+^} <[ami-iMv$(*)>- 

Z r ZTT * — ' I ) nm 

m=0 

(355) 

Inserting Eqs. ()353j) and ()355j) in Eq. ()348j) yields the special conformal anomaly 

7 C (°)(Z) = -b(l)^ + w. (356) 

The u7-term is /-independent and requires an explicit calculation. To this end it is convenient 
to go over from the matrix w nm to the kernel [w(u, v)]+ defined as 

fdu C 3 J 2 (2u - 1) [w(u, v)) + = V WnmCT^V - 1) . (357) 

Jo m=0 

The result for the kernel can be obtained by using non-local light-ray operators instead of O n ; 
(see for instance Eq. (J53Jl ) and reads 

[w(u,v)] + = w{u,v) — 5{u — v) I dtw(t,v) + —S(u — v) / dt (t — v)w(t,v) (358) 

Jo du J 

I \ n nl \ U 2 f U — > 1 — U 

w(u,v) = -C F 9{v - u)— — + < 

v [u — v) z I v — > 1 — V 

where the "+" subtraction is extended in such a way that the second order pole, appearing 
at u = v, is regularized, cf. Eq. (j!42|) . The u>-term is absent in a scalar theory [TH] and 
is responsible for the deviation of the explicit NLO calculations in QCD of the off-forward 
evolution kernels in the MS-scheme from the expectations based on the conformal operator 
product expansion PBl HH1 1114j . Forming moments with respect to Gegenbauer polynomials 
(ED one obtains [H]: 



w nrn = -C F 6(n>m) (l + (-l) n ~ m ) (2m + 3) (359) 

x 



+ -g^P— Ty {A nm - V»(n + 2) + ^(1)} 
[m + l)(m + 2) 



where A nm = if> (2±^±4) - ^ (s=2) + 2ip (n - m) - if) (n + 2) - ^(1). This form of the special 
conformal anomaly matrix and the explicit expression for the [w(u,v)} + kernel are valid for 
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all twist-two quark operators, irrespective of their Dirac structure. Of course, the anomalous 
dimensions for tensor operators differ from those of (axial-)vector ones. 

Let us add that the calculation of conformal anomalies can be systematically extended to 
higher loops. In abelian theory this can be achieved by introducing renormalized operator 
products such as 

n 

[0 A O nl ] = [O a ] [O nl ] - {Z~ A } nm Pmi-i] , (360) 

m=0 

where O a = J d d xx~ G a ^ u (x)G a ^ u (x) and Z A are the corresponding additive renormalization 
constants. In QCD additional counterterms appear due to peculiarity of the coupling constant 
renormalization for non-abelian gauge fields, cf. Eq. (jB.25|) . Within the definition (J360|) the 
quark-quark entry in the special conformal anomaly reads in Landau gauge 



7 



c(0) 



2 [ZM°\ b] + Z [ A im - , with = V } . (361) 



2 



Obviously, the renormalization constant Z A = —2Z^°w + w appears on the r.h.s. of Eq. 
(j355p . 

We are now in a position to discuss the anomalous dimensions in NLO. The special conformal 
anomaly ()356|) inserted in the conformal prediction ()346|) yields together with the forward 
anomalous dimensions the full anomalous dimension matrix to NLO accuracy 

7 (D = 7 D(1) + 7 ND(1) with 7 ND(1) = _ ^(0)^ ^ (^J + 7 (0)) + [.(0)^] j (362) 

where we have introduced the matrices d and g defined as d nm = b nm /a(n,m) and g nm = 
w nm/ a ( n i m ) ■ The corresponding evolution kernel (u, v) is defined by its conformal moments 
as 

„i n 

/ duC^\2u-l)V^\u,v) = --Y,ln 1 icT(2v-l), (363) 

and for the present case it was calculated explicitly in J7|. For completeness, we now present 
the method for the construction of this kernel from the anomalous dimensions (|362j) . This 
method is general and can be used in other cases as well. 

The matrix multiplication in Eq. (|362|) . e.g., 'j^g, corresponds in the momentum fraction 
representation to the convolution of the kernels — 2(V° ® g)(u,v) = —2 J dt V^°\u, t)g(t, v), so 
that all we need to do is to find the kernels g(u, v) and V(u, v) that give rise to the conformal 
moments g nm and [y^ \d] nm , respectively. The precise definition is 

9(u,v) = E E ( 3 } ^ ~ u)uCT{2u - l)g nm CT(2v - 1) (364) 
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and similar for V(u,v), with the replacement g nm — > [7^, rf] nm . 

First, we construct the g(u, v) kernel. Since the Gegenbauer polynomials obey the equation 



^(1 - u)vC%*(2u - 1) = -(n + l)(n + 2)C*/ 2 (2 W - 1) , 



(365) 



the multiplication with a(n, m) /2 = (n+l)(n+2) — (m+l)(m+2) in the conformal moment space 
corresponds to a second order differential operator in the momentum fraction representation. 
Thus, from the definition a(n,m)g nm = w nm it follows that g(u,v) satisfies the inhomogeneous 
second order differential equation 



d 2 d 2 1 

—-—(l — v)vg(u,v) — (l — u)u—-^g(u,v) = -w(u,v) for u^v. 
dv 2 du 1 2 



(366) 



The homogeneous solution of this equation is diagonal with respect to Gegenbauer polynomials 
and is not of interest for our purpose. The particular integral which contains the conformal 
moments g nm for n > m, reads 



g{u,v) 



-C f 



9{v-u) 



v — u 



1 - u 

1 - V 



(367) 



where the "+" subtraction is defined in Eq. 

Second, we construct V(u,v). It turns out that the matrix elements d nm = b nm /a(n,m) 
occur in the expansion of the Gegenbauer polynomials with a shifted index z/ = 3/2 + ein terms 
of those with index v — 3/2 at first order in e: 



C^(2n - 1) = Cl' 2 (2u - 1) - 26 V d nm C 3 J 2 (2u - 1) + 0(e 2 ) , 



m=0 



(368) 



where we are only interested in d nm for n > m. To show this, use the defining differential 
equation for Gegenbauer polynomials 

^(l^) M Cf 2+e (2u-l)-e|-(2u-l)Cf 2+ \2u-l) = -(n+l)(n+2+2e)C^ 2+e (2«-l), (369) 



and the definition (|344j) for the 6-matrix. The desired relation d nm = b nm / 'a(n,m) then follows 
upon the expansion in e. Thus, we look for the kernel V(u,v\e) with eigenfunctions Cn 2+e and 
eigenvalues 7n°' ) for e = 0. This kernel must be symmetric with respect to the weight function 
[(1 — u)u] 1+t , i.e. [(1 — v )v] 1+t V(u, v |e) = [(1 — u)u] 1+e V(v, u\e) and can easily be deduced from 
the LO evolution kernel (|142|) by including an additional factor (u/v) e : 



V{u,v\e) = d 



[v - 11) I ^ I I 



v — u 



1 - u 
1-v 



(370) 



The derivative of V(u, v\e) with respect to e leads to the dotted kernel, which is therefore given 
by a logarithmical modification of V^°\u,v): 



V(u,v) 



0_ 

7k 



V(u, v\e) 



C f 



e=0 



U 



9(v-u)- 1 + 



u 



1 u 
In - + 

v 



(371) 



Collecting everything, we find the following structure of the evolution kernel at NLO 

V(u,v;a s ) = ^V^(u,v) + J ^V^(u,v)+0(al), 
V w (u, v) = D(u, v) + G(u, v) - \v ® V {0) -—V + g® V {0) - V w ® g\ (u, v) 



(372) 



The first two terms D(u,v) and G(u,v) only possess diagonal Gegenbauer moments and are 
distinguished by color factors. They can be restored, at least in principle, through an integral 
transformation of the two-loop DGLAP kernel [22] where, to avoid double counting, one has to 
subtract the diagonal contributions that have already been included in the last term in ()372j) . 
The D-kernel can be represented as a convolution of the two separate diagonal parts of the 
LO kernels (jl42J) . while the F- kernel is related to the crossed-ladder two- loop diagram and is 
rather complicated. It is given by the rational function (u/v)[l + l/(v — u)] which enters the 
LO kernel (|142|) . decorated with Spence functions and double logs. The explanation that the 
crossed-ladder diagram is diagonal is rather simple: it contains no subdivergence and hence 
only a single 1/e pole. Consequently, only tree-level counterterms are required that preserve 
conformal symmetry. The explicit expressions can be found in Ref. |115j . 

The NLO evolution kernel for the tensor quark- ant iquark operator can be constructed along 
the similar lines, see |115|lll6j . It has the same structure ()372|) with g(u, v) defined in Eq. ()367j) . 
but the LO kernel is of course different: 



yT(°)( 



u. V 



C f 



9{v 



u 

u) — 



V v — u 



\u- 


-> 1 


-u) 


\] 


\v- 


-> 1 


-vj 





-?-8{u-v) 



(373) 



The corresponding dotted kernel V T is given by the logarithmical modification of V T<y0 \ in full 
analogy to Eq. (J371j) . It turns out that G T can easily be read off from the known result for 
G(u,v) in the vector case by a simple replacement of the corresponding LO structures, e.g., 
- (l + ~~ > f («^~) • -^ or ^ ne remaining diagonal kernel D T (u, v) one just uses an ansatz, as 
a linear combination of diagonal LO kernels and their convolutions, and finds the corresponding 
coefficients from the known diagonal anomalous dimensions or, alternatively, the corresponding 
DGLAP evolution kernel to NLO. 

The similar procedure allows one to calculate the evolution kernel for gluon tensor twist-two 
operators |117[ lllfi] . where in addition supersymmetric conditions are used to reconstruct the 
corresponding gluonic G(u,v) kernel. Note that the quark and gluon tensor operators belong 
to different spin representations of the Lorentz group and, thus, do not mix with each other 
under renormalization. 
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Figure 9: The renormalization constant Z A at one-loop order. The notations are the same as 
in Fig. El 



5.3.1 The Flavor Singlet Sector 

The technique based on the evaluation of conformal anomalies can be extended to the fla- 
vor singlet sector and indeed it has been employed to calculate the corresponding anomalous 
dimension matrices [221 I2H] and evolution kernels at NLO |115t 11181 H19j . The formalism be- 
comes more cumbersome, since the quark (|fi4|) and the gluon operators (|fi5|) with the same spin 
mix with each other under renormalization. For convenience we introduce the two-dimensional 
vector 



O 



nl 



(id, 



\l—n 




with / > n, 



(374) 



and the 2x2 matrices: 



a 



a 
a 




7 



GC 



7 



GG 



7 



7 



QC 



GQryC GG^yC 



(375) 



Although the superscript of the corresponding Gegenbauer polynomials is different for quark 
and gluon operators in Eq. (|374j) . both of them have the same conformal spin. Hence the 
nonvanishing entries in the matrices a and b are the same in the quark-quark and gluon-gluon 
channels. They are given in Eqs. ()325|) and ()345|) . respectively. The following construction 
applies for the both existing sets of flavor singlet operators which can be distinguished according 
to vector (V) and axial-vector (A) Dirac structure of the quark operators. 

The constraints for the conformal anomalies in the flavor singlet sector have been derived 
in [221 12H]- To achieve the same form as in the non-singlet sector, Eqs. (|331j) and ()343)1 . 



a(/)+7 c (0 + 2^5(0,7 
9 







and 7 C (Z + 1) - 7 C (Z) = -2 7 



(376) 
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we redefined here the original special conformal anomaly matrix (|B.31|) in the gluon-gluon 
channel by a shift GG ^ C — > GG ^y c -f b/3(g)/2g. Consequently, the special conformal anomaly in 
LO reads (in Landau gauge) 



7 C (°) = 



2 [QQ^[i](o) ; h + qq^M(o)- _ 2 ^QG^[i](o) 



2GQzm(o)h 



(377) 



The matrix of renormalization constants enters in the renormalized operator product 
[£> A ][O ni ] where 0~ A = fd d xx~G* u (x)G aiM/ (x), see Eq. (Il3~25l) . It has been calculated at LO 
for both vector and axial- vector cases in [221 12H] • At this order the result in the flavor singlet 
quark-quark channel coincides with the non-singlet result in Eq. ()359|) . Since in the quark- 
gluon mixing channel there is no gluon loop, the operator insertion [O^] cannot couple, and, 
consequently, ^ G Z A is vanishing. It remains to calculate the anomalies in the gluon-quark, cf. 
Fig. Etd,e), and gluon-gluon channels, cf. Fig. The result reads 



c(0) 



7 

abJo) 



b 7 (0) + w 



AB„. 



(378) 

can be found in Refs. [221 H17j . Here 
! w nm are universal, i.e., the same for V and A operators, and are given 
by the conformal moments of the corresponding AB w(u, v) kernels. The quark-quark kernel 
®°w(u,v) = [w(u,v)} + is defined in Eq. ()358[) . ® G w(u,v) = and the remaining two kernels 
read 



where the explicit expressions for y n and 
we only note that AB 



u,v) 



GG„ 



'w(u, v) 



2C F 
-2C A 



6{v 



— U) I u 

V \v 

2 



u 2 9(v 



u 



v 2 {y — u) 2 



l-u 
1-v 



— 9(v - u) + -8{u - v) + 



1 - u 
l-v 



(379) 
(380) 



The constraints (J376|) immediately yield the off-diagonal entries of the anomalous dimension 
matrix at NLO. It has the same structure ()362j) as in the flavor non- singlet sector 



7 (1) = 7° (1) + 7 ND(1) with 7 



ND(1) 



7 



<°> d 



/3ol + 7 {0) ) + 



(381) 



where we introduced the matrices d nm = b nm /a(n,m) and g nm = w nm /a(n,m). These func- 
tions are the same for the vector and axial- vector operators while the corresponding forward 
anomalous dimensions are different (and known to NLO). This result implies the following 
structure of the evolution kernels 



V {1) (u, v) = D(u, v) + G(u, v) - | V ® V (0) - yF + g® V (0) - V (0) ® 



u,v) 



(382) 



where the off-diagonal conformal moments of the kernels V and g coincide with the matrix 
elements of [7^, d] and g. Similar to the flavor non-singlet sector the g kernel can be recon- 
structed from w, Eqs. ()379[) and ()380|) . as the solution of a certain inhomogeneous second order 
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differential equation. The result reads 



{u — > 1 — u 1 
, 383 
D — > 1 — V J 

where the +(— ) sign is to be taken in the diagonal (off-diagonal) entries, respectively. In 
the axial- vector case the dotted kernels are constructed from the LO evolution kernels by the 
logarithmical modification displayed in Eq. (|371|) . while in the vector case there are subtleties 
that induce extra terms |119j . The main challenge is the construction of the diagonal G(u,v) 
kernel. The explicit calculation can be avoided by using six supersymmetric conditions together 
with the two known entries in the quark-quark channel for vector and axial-vector channel, 
respectively. Here it is crucial that the G(u, v) kernel arises from the 1/e pole of crossed-ladder 
diagrams, which preserve both conformal and supersymmetric covariance. This calculation 
includes a number of technical details which will be skipped here. The interested reader can 
find them in Ref. Finally, the D(u,v) kernel has again a simple representation in terms 

of known diagonal kernels and can be fixed by comparison with the known forward anomalous 
dimensions or DGLAP evolution kernels. 

Finally, let us list several consistency checks that have been performed for the conformal 
predictions: 

• In the flavor non-singlet sector the conformal prediction 1)362)1 coincides with the explicit 
NLO calculation of the evolution kernel for the pion distribution [4*4"1 I19j. 

• In the singlet sector the (3q proportional terms in the anomalous dimension matrix ()381j) 
have been checked by the calculation of Feynman with quark bubble insertions |2*2*1 123j. 

• There exist four supersymmetric constraints, arising from the reduction of QCD to su- 
persymmetric M = 1 Yang-Mills theory, for the off-diagonal part (|381|) of the anomalous 
dimensions. They are satisfied in a renormalization scheme that preserves supersymmetry 
|12Uj . As a byproduct it has been shown that two of six constraints for the LO forward 
anomalous dimensions derived in Ref. [TO], are modified beyond the leading order by 
off-diagonal entries. 

• Using superconformal symmetry |121[ 1122) I123j . which is anomalously broken, one is 
able to derive four constraints for the special conformal anomaly. Evaluation of the 
superconformal anomaly in LO shows that these constraints are consistent with the special 
conformal anomaly ()378)1 124] . 

• Rotating the conformal OPE prediction for the hard-scattering amplitudes, valid in the CS 
scheme, to the MS scheme (see next Section) allows one to resum the conformal partial 
waves. In this way one finds a closed expression for this amplitude in the momentum 
fraction representation ()404|) [2U These results coincide with the diagrammatical 
evaluation in Refs. |125| 1126] . 



g(u, v) = 6{v — u) 



-C f 



"In 1-^ 



-G, 
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Let us summarize the results of this Section. The complete set of two-loop anomalous dimen- 
sions and evolution kernels for all twist-two operators has been obtained using the conformal 
approach. A complete consistency check for the special conformal anomaly has been given in 
the considered order. Together with the validity of the conformal operator product expansion, 
shown at NLO, there is a complete understanding how conformal symmetry is implemented 
beyond the LO approximation and its predictive power provides a number of highly nontrivial 
results. 

5.4 Conformal Operator Product Expansion for Two-Photon Processes 

The aim of this Section is to illustrate the application of the conformal operator product 
expansion (COPE) considered in Sect. 2.3 beyond the leading order. The relevant physical 
setup is provided by generic hard exclusive processes including two photons, which are various 
cross-channels of the reaction 7* ((ft) + h(pi) — > 7 ( *)(<?2) + h(p 2 ) in the generalized Bjorken 
kinematics v = p- q — > 00, Q 2 = —q 2 — > 00 with £ = 1/ui — Q 2 /p-q fixed, where p — Pi +P2 and 
Q = (<?i + 92) /2 1451] . Here h(pi) and h(p 2 ) are generic hadronic states with small invariant 
masses (we also allow for the vacuum state) that can be neglected to power accuracy in 1/Q 2 . 
It is convenient to introduce the so-called skewedness parameter rj as the second scaling variable 
rj = A ■ q/p ■ q, where A = p 2 — p\. We assume that —t = —A 2 <C Q 2 and can be neglected as 
well. 

The QCD dynamics in the two-photon processes in the generalized Bjorken limit is described 
by the hadronic tensor 



which receives a dominant contribution from the region x 2 ~ 1/Q 2 - It can be evaluated using 
the operator product expansion of the two currents at short distances x^ — > 0, where to power 
accuracy it is sufficient to retain only leading twist-two operators. We have demonstrated in 
the previous Sections that the conformal covariance of operators can be preserved in QCD, up 
to corrections in /3(g), in the CS scheme. In this scheme, the COPE (J74J) is valid to all orders 
up to corrections proportional to the QCD (3— function. For transverse Lorentz projections of 
the currents — > J± we obtain, therefore 



where Q° s is the renormalized conformal operator with twist t n = 2 + 7 n (a s ) and conformal 
spin j n = n + 2 + ^7 ra (a; s ) in the CS scheme 8 . The hadronic matrix elements, entering Eq. 

8 For simplicity we do not include gluon contributions, and do not display the quark electric charges. 




(384) 



T(u,r),Q 2 ) 



(385) 




+ 0((3(g)), 
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(|385|) . can be parameterized as 



(h(p 2 )\Qn(ux„)\h( Pl )) = e ^)^ +1 ((Q c n s (0))}(v, A 2 ;/ u 2 ) , (386) 

where we used the identity (xA) = f](xp), valid in the Bjorken limit, and introduced the 
notation for the reduced matrix elements. The latter depend on the kinematic invariants 77, A 2 
and on the normalization scale fi 2 . By construction, the conformal operators have autonomous 
scale dependence and satisfy the renormalization group equations 

/^«Qn S (0)» = -7n(a fl M)«Q£ s (0))) . (387) 



We note that beyond the conformal limit both Eqs. ()385|) and (|387|) have to be corrected by 
terms proportional to (3{g) /g ■ 0(a s ). Such terms affect the anomalous dimensions to NLO 
order and the Wilson coefficients to the NNLO accuracy. They vanish in the forward limit. 
Substituting ()386|) into ()385|) and performing the x— integration, one arrives at 



T(^,r ? ,Q 2 ) = ^C n ^ / du 

n=0 W / JO 



2\7nM/2 rl ^+1^(1 _ w )ln+l+7»(0/2 



; (0)))(r/,A 2 ,^ 



[1 - 7]Uj{2u - l)]"+l+7„K)/2 

(388) 

where C n = C n -2 n+1 B(n+l,n+2)/ B(n+2+>y n /2,n+2+>y n /2); C n are the Wilson coefficients in 
the polarized deep inelastic scattering in the flavor non-singlet sector and are known to NNLO 



Cn(a s )=c^ + ^ + J ^cV + 0(a3) with = 1 . (389) 



The integral in (J388j) can be expressed in terms of the 2-^i~ hypergeometric series. It depends 
on the ratio of scaling variables = rjuj and resums the contribution of the conformal tower 
of operators with the highest weight Q^ s (0). Notice that for \r)u\ > 1 the integrand in ()388|) 
has a singularity inside the integration region indicating that the COPE (|388j) is divergent and 
should be defined by the analytic continuation from the \rju)\ < 1 region. It follows that the 
generalized structure function T(uj,r],Q 2 ) is an analytical function on the complex uo— plane 
with a cut that runs between the points £ = ±77. To begin with, let us consider the two limiting 
cases: r] = and 77 = 1. 

The first case, rj = 0, corresponds to the forward scattering A' 1 = 0. The u— integration 
becomes trivial, the ^-variable coincides with the Bjorken variable XBj, and the result for 
T(uj,r] = 0,Q 2 ) reduces to the well-known expansion of the forward Compton amplitude in 
the unphysical region \u\ < 1. In particular, C n (a s ) and 7„ coincide with the Wilson co- 
efficients and the anomalous dimensions in the deep inelastic scattering. Making use of the 
dispersion relations, one can identify the coefficients in front of powers of to in the expansion 
of T(uj,t] = 0,Q 2 ) as the moments of the imaginary part of the Compton amplitude in the 
physical region < XBj = 1/u) < 1- 
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The second case, 77 = 1, is realized, e.g., for the photon-to-pion transition form factor, 
7* + 7* (92) — > 71-0 ip), in which case the matrix elements are evaluated between the vacuum 
and the pion state with momentum p M = A^ 1 = qf + q^. Since the both photons are in the initial 
state, one has to replace q2 — > — qi- The corresponding hard scale is Q 2 = —(qi — q-i) 1 ! /4 = 
— (qf + ?f)/2. Due to parity conservation the hadronic tensor must be of the form 

✓y2 ~2 

T^uj, i] = l, Q 2 ) = ie 2 e fll/oc pq%q$F^(uj, Q) , lo = -±- — § , (390) 

9l + <?2 

where F 17T (uj,Q) defines the form factor and e is the electron charge. In this case it is also 
sufficient to consider transverse components of the currents and we can overtake Eq. (|388|) at 
77 = 1 and conformal operators defined in the axial-vector sector. Taking care of the appropriate 
normalization, the COPE result becomes 

M".G) = ^r E CnH^, Q/n) 3(w * * 2) + 0{i/Q 4 ) , (391) 

3Q 2 j^< 2(2n + 3) 

n even 

where 

r , , n/ v _ 2" +i a n (a s )E(n + l,n + 2) f fjT\ M/2 f l oo n [u(l - «)]" +1+ >' a ^ 2 
O n Ma s , s(n + 2 + 7j2) n + 2 + ^ /2) ^ g2 J ^ d M {i _ ^ u _ 1)]n+1+7n(as)/2 , 

(392) 

f n ~ 132 MeV and the scaling variable |u;| < 1 is called asymmetry parameter. The nonpertur- 
bative parameters n (/x) are related to the reduced matrix elements (7r(p)|Q^ s (0)|0) and define 
the partial wave decomposition of the pion distribution amplitude (|151|) . In the CS scheme 
4> n {v) are multiplicative renormalizable to any order of perturbation theory, up to corrections 
in (3— function. To the NLO accuracy, they obey the renormalization group equation with 
inhomogeneous term proportional to /3q: 



^) (0) , "»_,(!) 
2vr Jn (2vr) 2 Jn 



2 / n •> , n— 2 



m=0 



After the consequent expansion with respect to a s , the partial wave decomposition (|391|) is 
exact in NLO and after rotation (|333|) reproduces the MS result 

To solve ()393)1 one introduces a new set of "renormalization group improved" conformal 
operators that do not mix under renormalization to NLO: 

d n 
H—tfMin) = - 7 nK(//MGu) with <P'M = B^^M , (394) 

m=0 

where 7„(a s (/i)) is given by the two-loop expression defined in ()393|) and the mixing matrix B 
satisfies the differential equation 1128j 

A *-^B nm 0*)= fc), 7 D (a s (/i))l -ft^/A^)} + 0(c£). (395) 

tt/i L Jnm (Z7TJ Z I Jrnn 
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Figure 10: The perturbative QCD prediction for the scaled pion transition form factor 
2Q 2 F.y n (uj, Q) for one virtual and one real photon 



1) using the asymptotic pion dis- 
tribution amplitude (n = 0). The calculations in LO, NLO, and NNLO are shown by the 
dotted, dash-dotted, and solid lines, respectively We use a^ s (/i = Mz) = 0.118 and set the 
scale [i 2 = 2Q 2 = —q\. The data are taken from Refs. |129| I60|. 



Here 7° m (a s ) = S nm j n (a s ). Since the pion transition form factor does not depend on the 
factorization scale /i, the rotated Wilson coefficients defined as 



C' n (u\a s (fi),Q/fi) = C m (u\a a ([jL), Q/n)B mn {n) 

m=n 

satisfy the renormalization group equation 



(396) 



^C' n (u\a s (n),Q/fi) 



2^ 



7i 0) + 



(2tt) 



) 7n 



i 1} + 0(a 



C' n (u\a B (iJ,),Q/(i) . 



(397) 



The perturbative solution of this equation to NLO can be cast in the form of Eq. ()392)1 by 
modifying the Wilson coefficients j2Ej 



CJa, 



<*M a) , w) (2) _ qM_ Pq ln (Q_ 

2tt n (2tt) 2 n 2tt 2 



x 



(0) 



c y 1 4- c v 



2tt 



+ c (o)^(o)_^_ 

u n fn „ (o) 



+ 0(c* 



3 ) 



(39£ 



where the initial condition is chosen in such a way that the conformal covariance holds true for 
fi = Q. Here the derivative in 'j^ acts on C n (u\a s , Qj //) in Eq. ()392|) and amounts to replacing 
a s (fi) —>■ J® d\nfi'a s (fi') in the one-loop expression for the anomalous dimension 7 n (a a ). 

The experimental data are available for the quasi-real photon limit q\ — > which corresponds 
to = 1 and Q 2 = —qf/2 \129[ loT)] and they are reasonably well described, see Fig. [TUJ by 
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the contribution of the lowest conformal partial wave alone. This can easily be found from 
Eqs. (JHU) and flUED 9 



F 77r (u,Q) 



V2f w 
2Q 2 



71 



71^ 



3.583 - 2.25 In 



I* 



+ 0(a 3 s )\, (399) 



where we used Nf = 3. The overall normalization of the LO result is proportional to the sum of 
all conformal partial waves 1 + 02 + 04 + • • • so that the data do not leave much room for higher 
terms, unless there are sign alternating contributions. For now, this is the strongest existing 
evidence that the pion distribution amplitude already at relatively low scales is sufficiently close 
to its asymptotic form. 

It might be possible to measure the pion transition form factor for both virtual photons for 
a certain range of |u;| < 1 at the running facilities CLEO, Barbar and Belle, see e.g. Ref. |13()j . 
This would help to separate contributions of different higher conformal partial waves, and also 
provides a novel QCD test for exclusive processes that is parameter free fl30t 12l>] : 



Q 2 



duF ln (uj, Q) 



V2U 



1 - 



a s (Q) 



3.583 



a 2 s (Q) 



71 



71- 



20.215^^ + 0{af) 



7T" 



where u cnt < 0.4 is required and we set \i = Q. 

Returning to the general case < r] < 1 we recall that the COPE expansion (|388|) diverges 
in the region \r)u\ > 1. This is the similar situation as in deep inelastic scattering where 
the amplitude can only be restored from the moments by analytic continuation through the 
Mellin transform. In the present case we also have to introduce generalized parton distributions 
[JU HB] given by the Fourier transform of light-ray operators (|59j) 



q(x,r], A 2 ,/i 2 



pe^ n - p \h(p 2 )\^(-an) 1+ r^(an)\h( Pl )) 

Z7T 



n-A=rjn-p 



(400) 



where T = {1,75} for vector and axial-vector contributions to the hadronic tensor and the 
variable x has a partonic interpretation of the momentum fraction. Using the definition of the 
conformal operator, Eq. ()64j). one can identify the reduced hadronic matrix elements (j386|) as 
conformal moments of the generalized parton distributions 



(401) 



In the limit of forward scattering rj, A 2 — > this relation coincides with the definition of con- 
ventional parton distributions whereas for the pion transition form factor, 77 = 1, it reproduces 
the definition of the pion distribution amplitude. 



9 Thanks to current conservation, the anomalous dimension 70 vanishes. Thus, C${uj = 1) = 3co/2, where 
Co = Csj enters the Bjorken sum rule Q3l)8fl . 
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The evolution equation for q(x, 77, A 2 , fj, 2 ) follows from the renormalization group equation 
for the light-ray operators |351 1131|. I44[ I45j. In the flavor non-singlet sector it looks like 



2 / a 2 2 



% yf v + x v + y . 
277 ' 277 



-i2N 



(402) 



By construction, at 77 = 1 this equation has to reproduce the ER-BL evolution equation 
for the pion distribution amplitude. This implies that for v = (1 + y)/2, u — (1 + x)/2 with 
< u, v < 1 the kernel V^(m, v) coincides with the ER-BL kernel. Using the well-known result for 
the latter one can restore the expression for the V— kernel entering (J4U2|) |4"H I45j . Substituting 
(|4U1|) into (|388|) allows for the (formal) resummation of all conformal partial waves and results 
into the momentum fraction representation of the hadronic tensor. In this representation the 
hadronic tensor is factorized into a convolution of the hard scattering amplitude 



- = jrC n (u JV \a s ,Q/riC 3 J 2 (-) 



(403) 



and the generalized parton distributions (or generalized distribution amplitudes for the crossed 
processes): 



J-i M 



T(i-,-;«.(A0)T {*--£} 

1 T] T] /i ' 



q(x,r],A 2 ,fx 2 



(404) 



where — (+) is for the (axial-)vector case. Unfortunately, for the CS scheme exact expressions 
for the hard scattering amplitudes ()403|) and the evolution kernels are not known beyond LO. 
To the NLO accuracy they have been evaluated in the MS scheme by the rotation of the 
local COPE to this scheme and reconstructing the corresponding functions in the momentum 
fraction representation by analytic continuation |2U 123 ■ After such a resummation we can 
employ analytical properties of the Compton amplitude to obtain the result in the physical 
region. 



6 Outlook 

To summarize this review, the uses of conformal symmetry in Quantum Chromodynamics in the 
last years were developing along two different lines. The first line of applications is essentially 
algebraic. The parton model, by construction, retains all symmetries of the classical Lagrangian 
and the QCD corrections in the leading logarithmic approximation also do. The simplest and 
the standard application of conformal symmetry to light-cone dominated processes was (and 
remains) the separation of variables corresponding to the transverse and longitudinal degrees 
of freedom. This role is very similar to the partial wave expansion in nonrelativistic quantum 
mechanics. Remember that in the latter case the classification of states in contributions with 
given orbital angular momentum remains to be relevant also in the presence of spin-orbit 
interaction which breaks the rotational symmetry. By the same token, the classification in 
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terms of conformal spin in QCD remains relevant also with higher-order radiative corrections. 
A deeper look in the algebraic structure of the evolution equations in QCD reveals hidden 
symmetries and interesting connections with completely integrable systems, with the conformal 
spins of partons playing the role of spins of a generalized Heisenberg magnet. These connections 
are intriguing and are not really understood. On the theoretical side, they may allow for the 
studies of high-energy QCD using methods inherited from string theory and possible emerging 
of a dual string picture, see [3]. In this way the strong coupling regime of the evolution in 
super symmetric extensions of QCD can be addressed |132t 1133] . On the practical side, these 
developments already led to an almost complete understanding of the spectra of twist-three 
operators, which are relevant for the phenomenology. 

The second line of applications has been using constraints imposed by conformal symmetry, 
and in particular by conformal Ward identities, in order to advance perturbative calculations for 
processes involving off- forward kinematics to higher orders. As the principal result, there exists 
now a complete understanding how conformal symmetry is implemented beyond the leading 
logarithmic approximation and its predictive power provides a number of highly non-trivial 
results at the next-to-leading order. It is shown that the covariance of leading twist conformal 
operators can be ensured in all orders of perturbation theory in a special renormalization 
scheme. This so-called conformal subtraction scheme differs from a MS-like scheme by a finite 
renormalization. The conformal operator product expansion holds true in this scheme, and 
indeed it has been checked in NLO by the explicit calculation. The structure of the conformal 
constraints is such that, generally speaking, n-loop predictions for off-forward amplitudes can 
be obtained from the similar predictions in forward kinematics and the n — 1-loop calculation of 
the special conformal anomaly. The technique is general and can be applied in a broad context. 

There are several other interesting issues that we were not able to cover, in particular 
using conformal symmetry in semiclassical construction of the QCD vacuum, e.g., based on 
instantons. This is, however, another subject. 



Appendices 

A The Conformal Basis 

For a generic three-particle operator a "conformal basis" can be constructed as follows [3*01 13*7] . 
The conformal symmetry allows one to fix the total three-particle conformal spin J = ji+ j% + 

( 12)3 

js + N of a state. We define a set of functions Yj- by requiring that, in addition to a fixed 
J, they also have a definite value of the conformal spin in the given two-particle channel (12), 
for definiteness: j = j\ + ji + n with n = 0, . . . , N. 

Taken together, these two conditions determine the polynomials Yj^ 3 uniquely and yield 
the following expression: 
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Here Pn (x) is the Jacobi-polynomial. The basis functions Yj- (ui) are mutually orthogonal 
with respect to the conformal scalar product: 

[\du] u\ h - 1 Y^\ Ui )Y^)\ Ui ) = Afjj 5jj, S jf , (A.2) 
Jo 

where 

M ■•(./ •./: j-i)\'U -./i -.hi) ■'(•/ ./ •./V.jl'i-/ •./ ./VJ 

33 r(j-j 1 -j 2 +i)r(j+j 1 +j 2 -i)(2i-i) r(j-j-j3+i)r(j+j+j3-i)(2J-i)' s 

(A.3) 

The above construction of the conformal basis involves the obvious ambiguity in what order to 
couple the spins of partons to the total spin J. Choosing a different two-particle channel, one 
obtains a different conformal basis related to the original one through the matrix Q of Racah 
6j-symbols of SL(2, R), e.g., 

i Ji 1)a («*)= E ^A^Yjf 3 ^). (A.4) 

jvti2<j'<J-33 

The properties of the Racah 6j-symbols as well as explicit expressions in terms of the generalized 
hypergeometric series 4^3(1) are summarized in p?7] . 



B Conformal Ward Identities 

To derive the conformal Ward identities (CWIs) we employ the generating functional 

Z(j) = i- J £>$ exp|z[5]+i J d d xJ{x)$(x)X with A/" = J £>$ exp{i[S]} , (B.5) 

where [S] = Jd d xC(x) is the renormalized action and the dimensional regularization is used, 
i.e., d = 4 — 2e. The renormalized connected Green functions result from a functional derivation 
with respect to the sources J and one-particle irreducible Green functions, which are not con- 
sidered here, are obtained from a Legendre transformation. Replacing the integration variable 
$ by the new one where the change in the integral measure is given by the Jacobian 

Det [X>$'/X>$], does not alter the generating functional (|B.5j) . Thus, for an infinitesimal con- 
formal transformation (jUJ) 

= $(x) + e5 G $(x) with <5 G $(x) = g(x, <9)$(x) for G = {P M , M M „, D, K^} , (B.6) 

one finds that the variation of the generating functional (JB.5J) has to vanish to order e: 



% J d d x3(x) {g(x,d)$(x)} + i5 G [S} -i8 G N exp |z[5] +i J d d x J(x)$(x)| =0.(B. 



•7) 
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Here we could safely neglect the Jacobian Det [I + eQ], since it is only a C-number and cancels 
out between numerator and denominator, cf. Eq. (jB.5|) . The variation of the normalization 
SqM = J 5g[S] exp {i[S}} is understood as subtraction of the vacuum expectation value 
of the operator insertion 5g[S] and yields the normal ordering :5g[S]:= 5g[S) — 8qN- For 
simplicity, we do not indicate this operation here and so Eq. (jB.7|) can be cast in the form 

Functional derivatives with respect to J give rise to to the CWIs for connected Green functions. 
Thanks to Poincare invariance, the variations 8p[S] and 8^ [S] vanish identically, while the 
conformal variations Sd[S] and 8^[S] contain the trace anomaly, introduced in Sect. I2.4TT1 In 
what follows we calculate Sd[S] and 8^[S] for QCD. 

Within the covariant gauge fixing the renormalized QCD Lagrangian reads 

c = z# w - - i (^i;) 2 + z w . (b.9) 

Here u and u are the anti-ghost and ghost fields, respectively. The renormalized covariant 
derivatives in the fundamental and adjoint representation are = d fl — in e gXA'j 1 t a and D ab = 
8 dp + fi e gX f acb A c ^ , respectively, and the renormalized field strength tensor is given by G a = 
d^A® — d u A® + fJ, £ gX f abc A b ^A c u . Here \x is a mass parameter which plays the role of the renor- 
malization scale and ensures that the renormalized coupling g is dimensionless. The Z§ factor 
renormalizes the corresponding field: $ unr = while X is needed for the renormalization 

of the coupling: g nnr = fi e gX/ \tZa- We remind that the action is renormalization group invari- 
ant, i.e., the unrenormalized and renormalized actions are equal: S = [S]. In a MS-like scheme 
all renormalization factors are given by Laurent series 

Z = l + ±zW(a,(jA),£(ji)) + ^Z^(a s (fi),^)) + 0(l/e 3 ) for Z = {Z^ Z A , Z u , X} (J3.10) 

where the coefficients Z^(a s (fx), do not explicitly depend on fi (see for instance Ref. |134j ). 
Note that the gauge invariance of the Lagrangian (|B.9)1 is broken by the gauge fixing procedure, 
however, it is invariant under the following renormalized BRST transformations 

-i/i e gXZu u a t a ^5X , 5 BRST ^ = Z w D^ a 5X , 
l -^gXZJ abc uj b uj c 8X , 5 BRST iu a = ^d,A*5X , (B.ll) 

where 5X is a Grassman variable. 

The conformal variations of the renormalized action can be calculated using Eqs. (|92|) and 
(JUS): 

{% } [S]= I ^ {2^ } A(X) With A(X) = Ad(x) ~ k^xW (B - 12) 



5 BRST 4> = - 
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where Ad(x) and the addenda A v K (x) are defined in Eqs. (|94j) and (|98|). As we have mentioned in 
Sect. 12.4.11 A u K (x) does not vanish for the gauge fixing and ghost terms even in four dimensional 
Minkowski space. On the other hand, it is possible to achieve that both conformal variations 
are expressed in the physical sector in terms of the single trace anomaly Ad(x). This is done by 
choosing the scaling dimensions of physical fields in d dimensions to be equal to the canonical 
ones in four dimensions, whereas for the ghost fields we make the following choice: 

£^ = 3/2, £ A = 1, 4 = d-2, and C = 0. (B.13) 

That is, l a is simply the canonical dimension in d dimensions and then the choice £ u = is 
motivated by the BRST transformations ([B.lip . 

Moreover, it turns out that the choice (IB.13|) allows for a straightforward decomposition of 
A(x) in gauge invariant (type A), BRST exact (type B), and EOM operators: 

A(x) = -e {O a (x) + O b (x) + Q a {x) - n^{x)} -(d-2) d»(D Bv (x) . (B.14) 

Corresponding to their subscript, the operator insertions 

v rBRST rBRST 

O a (x) = -f (G%) 2 , O b (x) = —u,"d,A« , O b ,(x) = ~^-^Al , (B.15) 
belong to class A and B and the EOM operators read 



The l.h.s. of the CWI (|B.8jl is finite by definition and so we conclude that the conformal 
variations (jB.12|) are finite operator insertions. This implies that the operators in Eq. (jB.14|) 
generate a 1/e pole, which is annulled by the e prefactor. In the following we renormalize these 
operators. 

The general mixing scheme is the following: class A operators need themselves, class B as 
well as EOM operators as counterterms, while class B operators can only mix with themselves 
or EOM operators |135j . Of course, only operators with the same quantum numbers and the 
same or lower twist can appear as counterterms. The EOM operators are renormalized, since 
they correspond to a partial integration in the generating functional (jB.5|) . The renormalization 
problem at zero momentum transfer can be completely solved by means of differential vertex 
operator insertions, which are generated by the derivation of renormalized Green-functions 
with respect to g and £. [136 . Thus, they are renormalized and this allows us to find the 
integrated and renormalized operator insertions [O a ] and [O b ] |4"9"1 15T)]: 



[A 9 ] = gj-[S\ = [O a ] + [O b ] +f2 A + f2 a , [A*] = ^[S] = 1 -{\0 B ] + ^} , (B.17) 



where operators without argument are defined at zero momentum transfer, for example 12$ 
fd^ 12$ (x). 
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These results together with the general renormalization properties allow us to express the 
conformal anomaly in terms of renormalized operator insertions ^3J 1221 l2Hj : 

A(x) = ^{O a {x)\ + 7a) {{Ob{x)\ + f2 a (x)} - ( 7 a - ^ (2 A (x) 

-( TV) - e)Q^{x) - 2y u Q (x) - (d - 2)d fl [0»{x)) . (B.18) 

Here, 7$ = fi-^ In \fZ~4, denote the anomalous dimensions and /3 e (g;e) = —eg + (3(g) is the 
(3 function in d dimensions. Inserting this result into Eq. (|B.8|) and taking derivatives with 
respect to the sources J we obtain the renormalized CWIs: 

- N B 

VA(4> = -V7*(^7v)--(^[A 9 ]^)-a^[A«]^), (B.19) 

N N 

K (Xn) = - Yl ^ 2 < - -m gi "W - <r(i[A*»]X N ) - 4(i[0%\X N ) . 

i=l i=\ " 

Here we could safely set e to zero, a = —2^ A is the renormalization group coefficient of the 
gauge fixing parameter and in analogy to the differential operator vertex insertions (jB.17|) we 
defined 

[A£] = J d d x 2x, {[O a ] + [O b ] + f2 A + f2 B ) (x) , [A*] = / d d x 2x^~ {[O b ] + Q Q ) (x) .(B.20) 

As the next step we consider the CWIs with one operator insertion O n i. These are derived 
along the same line from the generating functional 

Z(J\O nl ) = ^J V<$> [O nl ($)}expUs}+i J dk3{x)Q(x)\ (B.21) 

and read for the connected Green functions 

([O^] (5 G X N )) = -((5 G [O nl ]) X N ) - ([O nl ]i (8 G [S]) X N ) . (B.22) 

For simplicity we assume that the composite operators are closed under renormalization, i.e., 
we neglect possible counterterms of class B and EOM operators: 

[O nl ] = ^Z nm (a s ,e)O m * with Z(a s , e) = i + -Z [1 \a s ) + 0(1/ e 2 ) . (B.23) 

771=0 



The l.h.s. of the CWIs (IB. 22)) is finite by definition and so its r.h.s., too. The conformal 
variations 5c[O n ;] of the renormalized operator insertion are the same as at tree level, except 
for special conformal transformation, which yields a divergent expression: 

n 

S D [Om] = (I + 3) [0„,] , 6 K [O^] = -iV Iza^Z' 1 } [O m ,_!] . (B.24) 

771=0 
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The divergencies in Za{V)Z~ x have to be annulled by UV- divergencies that arise from the 
operator product [O n i] 8^ [S] . In the case of dilatation we conclude that [O^^ofS'] is finite. 
However, also this operator product will induce new UV- divergencies that are multiplied by 
e, see discussion above. The UV- divergencies of these operator products are caused by the 
singularities in [O n j]A(;r) that are concentrated in x = 0, see also Eq. (jB.12|) . 

The details of the renormalization procedure [O n ;]A(x) can be found for Abelian gauge 
theory in Ref. ^3J CHI an d for QCD in Ref. [23]- Here we only present the gauge-invariant 
contribution of [O a {x)\, cf. (|B15jl : 



n . n 

i[0 A (x)][O nl ] = i[0 A {x)O nl ] - ^ D \x) V \Z A \ [O ml ] - l -d + ^ D \x) V jzj) [O^-x] 

m=0 m=0 

The renormalization is done by an additive subtraction, where counterterms with higher deriva- 
tives, denoted by the ellipses, will not contribute to the CWIs. A similar formula exists for 
the renormalization of the product [Ob + ^o](^)[O n j], which is already finite in Landau gauge. 
These ansatze are motivated by considering the differential vertex insertions which also relate 
the renormalization matrices Z A and Zb to the anomalous dimensions. The term containing 
the variation of the composite operators with respect to the gauge field ensures the correct 
renormalization group equation for the composite operators The products of conformal 
and EOM operators can easily be treated by partial integration in the functional integral (IB. 211) . 
Instead of the MS prescription, we define 

(i[O nl n$(x)]X n ) = —([O n i&(x)———X n ) , (B.26) 

o<P{x) 

where the unrenormalized operator product is related to the renormalized one as 

i[o nl ]n*(x) = i[o nl n^(x)) - $(a0— ^[On,] . (B.27) 

Employing Eqs. (jB.22|) and (jB.24|) - (jB.27|) the CWIs can be written after some algebra as 

N n 



J2^([Om}X N ) = -Y,{(l + 3)i + 7} nm ([O m iXN}--(i[O nl A 9 }X N ) (B.28) 

i=l m=0 " 

N 



N 



8 

£7 ([Oni]X N ) = i W) + 7 c } nm ([O m i-i]X N ) - ^[O^Ai]^) (B.29) 

i=l m=0 9 



N 

-a(i[O nl Ai)X N ) + 2(i[O nl A; xt ]X N ) - J2l^7([Oni)X N ) 

i=0 
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where A cxt = J d d x2n-xd ll O^. The finiteness of the anomalous matrices 

7 = lim \--Z A -(-- J A ) Z b + 2( 7 ^ - e)ZP Q Z- 1 ) , ( B.30) 



9 \9 



= 2^P Q - 2 [ZW P Q ] + ££ ] + ^ , where P Q = (\ ° j , 
= lim (z |a(Z) - 2( T ^ - e)P Q &l Z" 1 - ^Z^ - - lA ) Z~ - a(l) 



9 \9 

= -2^P Q b + 2 [ZW P Q b] + ZW- + Z§~ , (B.31) 

follows from the fact that all other terms in the CWIs (jB.28|) and (jB.29|) are finite by definition. 
In Landau gauge they are calculable from the counterterms in Eq. (jB.25|) alone, since Zb and 
Z^ are vanishing. Let us add that the definition ()B.30j) coincides with the usual one for the 
anomalous dimensions matrix: 



7 = 7z + 2 7A P G + 2 7v ,P Q with 7z = - {^Jj^j Z' 1 = 9^~ g Z [1] and P G = K ^(B. 



32) 
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